結果

問題 No.144 エラトステネスのざる
ユーザー kmjp
提出日時 2015-02-05 23:33:29
言語 C++11
(gcc 4.8.5)
結果
AC  
実行時間 360 ms
コード長 1,364 Byte
コンパイル時間 1,139 ms
使用メモリ 1,664 KB
最終ジャッジ日時 2019-07-22 02:51:05

テストケース

テストケース表示
入力 結果 実行時間
使用メモリ
sample01.txt AC 3 ms
1,652 KB
sample02.txt AC 3 ms
1,664 KB
sample03.txt AC 3 ms
1,656 KB
test01.txt AC 2 ms
1,644 KB
test02.txt AC 3 ms
1,644 KB
test03.txt AC 3 ms
1,644 KB
test04.txt AC 4 ms
1,660 KB
test05.txt AC 3 ms
1,660 KB
test06.txt AC 3 ms
1,656 KB
test07.txt AC 3 ms
1,656 KB
test08.txt AC 3 ms
1,660 KB
test09.txt AC 4 ms
1,656 KB
test10.txt AC 4 ms
1,656 KB
test11.txt AC 355 ms
1,660 KB
test12.txt AC 356 ms
1,656 KB
test13.txt AC 345 ms
1,656 KB
test14.txt AC 353 ms
1,660 KB
test15.txt AC 360 ms
1,660 KB
test16.txt AC 358 ms
1,660 KB
test17.txt AC 355 ms
1,660 KB
テストケース一括ダウンロード

ソースコード

diff #
#include <bits/stdc++.h>
using namespace std;
typedef signed long long ll;

#undef _P
#define _P(...) (void)printf(__VA_ARGS__)
#define FOR(x,to) for(x=0;x<to;x++)
#define ITR(x,c) for(__typeof(c.begin()) x=c.begin();x!=c.end();x++)
#define ALL(a) (a.begin()),(a.end())
#define ZERO(a) memset(a,0,sizeof(a))
#define MINUS(a) memset(a,0xff,sizeof(a))
//-------------------------------------------------------

// 2回以上行うならNP,prime,divを要クリア
const int prime_max = 10005;
int NP,prime[100000],divp[prime_max];
map<int,int> M;

void cprime() {
	for(int i=2;i<prime_max;i++) if(divp[i]==0) {
		//M[i]=NP;
		prime[NP++]=i;
		for(int j=i;j<prime_max;j+=i) divp[j]=i;
	}
}

int N;
double P;
double ret;

void solve() {
	int i,j,k,l,r,x,y; string s;
	
	cin>>N>>P;
	P=1-P;
	cprime();
	for(i=2;i<=N;i++) {
		y=x=i;
		double hoge=1;
		vector<int> V;
		for(j=0;prime[j]*prime[j]<=x;j++) {
			if(x%prime[j]==0) {
				r=0;
				while(x%prime[j]==0) r++, x/=prime[j];
				V.push_back(r);
			}
		}
		if(x!=i && x>1) V.push_back(1);
		y=1;
		ITR(it,V) y*=(1+*it);
		if(y==1) y=0;
		else y-=2;
		while(y--) hoge *=P;
		ret += hoge;
	}
	
	_P("%.12lf\n",ret);
}


int main(int argc,char** argv){
	string s;int i;
	if(argc==1) ios::sync_with_stdio(false);
	FOR(i,argc-1) s+=argv[i+1],s+='\n';
	FOR(i,s.size()) ungetc(s[s.size()-1-i],stdin);
	solve(); return 0;
}
0