結果

問題 No.144 エラトステネスのざる
ユーザー fmhr
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff # 

package main

import (
	"fmt"
	"math"
)

func main() {
	var N int
	var p float64
	fmt.Scan(&N, &p)
	var ans float64
	p = 1 - p
	d := 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 = 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 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
}
0