#466366 (C++17) No.144 エラトステネスのざる

• 問題: No.144 エラトステネスのざる
• ユーザー: fura
• 提出日時: 2020-04-18 02:53:21
• 言語: C++17 (gcc 12.3.0 + boost 1.83.0)
• 結果: TLE
• 実行時間: -
• コード長: 476 bytes
• コンパイル時間: 2,431 ms
• コンパイル使用メモリ: 205,924 KB
• 実行使用メモリ: 10,404 KB
• 最終ジャッジ日時: 2024-10-03 22:20:35
• 合計ジャッジ時間: 5,649 ms
• ジャッジサーバーID （参考情報）: judge1 / judge3

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

ソースコード

#include <bits/stdc++.h>

#define rep(i,n) for(int i=0;i<(n);i++)

using namespace std;

vector<long long> divisors(long long a){
	vector<long long> res;
	for(long long i=1;i*i<=a;i++) if(a%i==0) {
		res.emplace_back(i);
		if(i*i<a) res.emplace_back(a/i);
	}
	sort(res.begin(),res.end());
	return res;
}

int main(){
	int n;
	double p; scanf("%d%lf",&n,&p);

	double ans=0;
	for(int i=2;i<=n;i++){
		ans+=pow(1-p,divisors(i).size()-2);
	}
	printf("%.9f\n",ans);

	return 0;
}