代码
比这篇新的文章:
Codee#8148
比这篇旧的文章: USACO 5.2.3 wissqu
作者: , 点击99次, 评论(0), 收藏者(0), , 打分:
所有评论，共0条：( 我也来说两句)
比这篇旧的文章: USACO 5.2.3 wissqu
Codee#8147 编辑代码
语言: S, 标签: 无 2009/11/25发布 3个月前更新 | 编辑代码作者: , 点击99次, 评论(0), 收藏者(0), , 打分:
S语言: Codee#8147
01 p<-seq(0,1,by=0.001)
02 step <- seq(1,10, by=0.1)
03
04 half <- 0.5
05 half.percentile <- 0.3
06 ninty <- 0.9
07 ninty.percentile <- 0.5
08 eps <- 0.01
09
10 for (a in step) {
11 for (b in step) {
12 #cat (a,"\t",b,"\n")
13 product <- p^(a-1)*(1-p)^(b-1)
14 product <- product/sum(product)
15 m.diff <- abs(sum(product[p<half.percentile]) - half)
16 n.diff <- abs(sum(product[p<ninty.percentile]) - ninty)
17 if (m.diff <eps && n.diff < eps) {
18 cat (a,"\t",b,"\n")
19 }
20 }
21 }
02 step <- seq(1,10, by=0.1)
03
04 half <- 0.5
05 half.percentile <- 0.3
06 ninty <- 0.9
07 ninty.percentile <- 0.5
08 eps <- 0.01
09
10 for (a in step) {
11 for (b in step) {
12 #cat (a,"\t",b,"\n")
13 product <- p^(a-1)*(1-p)^(b-1)
14 product <- product/sum(product)
15 m.diff <- abs(sum(product[p<half.percentile]) - half)
16 n.diff <- abs(sum(product[p<ninty.percentile]) - ninty)
17 if (m.diff <eps && n.diff < eps) {
18 cat (a,"\t",b,"\n")
19 }
20 }
21 }
所有评论，共0条：( 我也来说两句)