結果
問題 | No.144 エラトステネスのざる |
ユーザー |
![]() |
提出日時 | 2015-06-13 01:32:03 |
言語 | Go (1.23.4) |
結果 |
AC
|
実行時間 | 78 ms / 2,000 ms |
コード長 | 1,032 bytes |
コンパイル時間 | 11,519 ms |
コンパイル使用メモリ | 224,968 KB |
実行使用メモリ | 16,304 KB |
最終ジャッジ日時 | 2024-10-10 19:03:47 |
合計ジャッジ時間 | 12,396 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
package mainimport ("fmt""math")func main() {var N intvar p float64fmt.Scan(&N, &p)var ans float64p = 1 - pd := divisorCounter(N)//fmt.Println(d)if p == 0{ans = float64(len(sieve(N)))}else if p == 1{ans = float64(N-1)}else {for i := 2; i<=N; i++ {div := d[i]if div==0 {ans+=1}else {ans += math.Pow(p, float64(div))}//fmt.Println(ans)}}fmt.Println(ans)}const MAX_N = 100000000var isNotPrime [MAX_N+1]boolfunc sieve(n int)[]int{prime := make([]int, 0)isNotPrime[0]=trueisNotPrime[1]=truefor i:=2; i<=n; i++{if !isNotPrime[i]{prime = append(prime,i)for j:=2*i; j<=n; j+=i{isNotPrime[j]=true}}}return prime}func divisorCounter(n int)[]int{d := make([]int, n+1)for i:=2;i<n+1;i++{for j := i*2; j<n+1; j += i{d[j]++}}return d}func divisor(n int)[]int{res := make([]int, 0)for i:=1;i*i<=n;i++{if n%i==0{res = append(res,i)if i!=n/i{res = append(res,n/i)}}}return res}