結果

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

テストケース

テストケース表示
入力 結果 実行時間
使用メモリ
sample01.txt AC 4 ms
1,660 KB
sample02.txt AC 4 ms
1,668 KB
sample03.txt AC 4 ms
1,660 KB
test01.txt AC 4 ms
1,648 KB
test02.txt AC 4 ms
1,648 KB
test03.txt AC 3 ms
1,648 KB
test04.txt AC 4 ms
1,668 KB
test05.txt AC 5 ms
1,660 KB
test06.txt AC 4 ms
1,656 KB
test07.txt AC 4 ms
1,660 KB
test08.txt AC 4 ms
1,660 KB
test09.txt AC 4 ms
1,656 KB
test10.txt AC 4 ms
1,664 KB
test11.txt AC 333 ms
1,664 KB
test12.txt AC 334 ms
1,660 KB
test13.txt AC 335 ms
1,664 KB
test14.txt AC 332 ms
1,664 KB
test15.txt AC 338 ms
1,664 KB
test16.txt AC 337 ms
1,664 KB
test17.txt AC 334 ms
1,664 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