結果

問題 No.144 エラトステネスのざる
ユーザー fmhr
提出日時 2015-06-13 01:04:55
言語 Go
(1.23.4)
結果
TLE  
実行時間 -
コード長 896 bytes
コンパイル時間 10,900 ms
コンパイル使用メモリ 232,444 KB
実行使用メモリ 13,768 KB
最終ジャッジ日時 2024-10-10 19:03:06
合計ジャッジ時間 14,572 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 11 TLE * 1 -- * 5
権限があれば一括ダウンロードができます

ソースコード

diff # 

package main

import (
	"fmt"
	"math"
)

func main() {
	var N int
	var p float64
	fmt.Scan(&N, &p)
	var ans float64
	p = 1 - p
	if p == 0{
		ans = float64(len(sieve(N)))
	}else if p == 1{
		ans = float64(N-1)
	}else {
		for i := 2; i<=N; i++ {
			//fmt.Println(i, divisor(i))
			div := len(divisor(i))
			if div<2 {
				ans+=1
			}else {
				ans += math.Pow(p, float64(div-2))
			}
			//fmt.Println(ans)
		}
	}
	fmt.Println(ans)

}

const MAX_N = 100000000
var isNotPrime [MAX_N+1]bool
func sieve(n int)[]int{
	prime := make([]int, 0)
	isNotPrime[0]=true
	isNotPrime[1]=true
	for 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 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
}
0