#32225 (Go) No.144 エラトステネスのざる

• 問題: No.144 エラトステネスのざる
• ユーザー: fmhr
• 提出日時: 2015-06-13 01:02:41
• 言語: Go (1.22.1)
• 結果: TLE
• 実行時間: -
• コード長: 894 bytes
• コンパイル時間: 10,895 ms
• コンパイル使用メモリ: 222,316 KB
• 実行使用メモリ: 15,068 KB
• 最終ジャッジ日時: 2024-10-10 19:02:50
• 合計ジャッジ時間: 14,757 ms
• ジャッジサーバーID （参考情報）: judge2 / judge3

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 13,764 KB
testcase_01	AC	2 ms; 6,816 KB
testcase_02	AC	1 ms; 6,816 KB
testcase_03	AC	1 ms; 6,816 KB
testcase_04	AC	1 ms; 6,816 KB
testcase_05	AC	1 ms; 6,816 KB
testcase_06	AC	2 ms; 6,816 KB
testcase_07	AC	3 ms; 6,816 KB
testcase_08	AC	3 ms; 6,816 KB
testcase_09	AC	3 ms; 6,820 KB
testcase_10	AC	3 ms; 6,820 KB
testcase_11	AC	3 ms; 6,816 KB
testcase_12	AC	1 ms; 6,816 KB
testcase_13	AC	1 ms; 6,816 KB
testcase_14	TLE	-
testcase_15	--	-
testcase_16	--	-
testcase_17	--	-
testcase_18	--	-
testcase_19	--	-

ソースコード

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))
			if len(divisor(i))<2 {
				ans+=1
			}else {
				ans += math.Pow(p, float64(len(divisor(i))-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
}