結果
問題 | No.144 エラトステネスのざる |
ユーザー |
![]() |
提出日時 | 2015-04-12 20:51:58 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 54 ms / 2,000 ms |
コード長 | 1,419 bytes |
コンパイル時間 | 689 ms |
コンパイル使用メモリ | 86,040 KB |
実行使用メモリ | 7,808 KB |
最終ジャッジ日時 | 2024-07-04 14:05:00 |
合計ジャッジ時間 | 1,721 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#include <cstdio>#include <iostream>#include <vector>#include <map>#include <set>#include <string>#include <cstring>#include <sstream>#include <algorithm>#include <functional>#include <queue>#include <stack>#include <cmath>#include <iomanip>#include <list>using namespace std;inline int toInt(string s) { int v; istringstream sin(s); sin >> v; return v; }template<class T> inline string toString(T x) { ostringstream sout; sout << x; return sout.str(); }template<class T> inline T sqr(T x) { return x*x; }typedef pair<int, int> P;typedef long long ll;typedef unsigned long long ull;#define For(i,a,b) for(int (i) = (a);i < (b);(i)++)#define rep(i,n) For(i,0,n)#define clr(a) memset((a), 0 ,sizeof(a))#define mclr(a) memset((a), -1 ,sizeof(a))#define all(a) (a).begin(),(a).end()#define rall(a) (a).rbegin(), (a).rend()#define sz(a) (sizeof(a))#define Fill(a,v) fill((int*)a,(int*)(a+(sz(a)/sz(*(a)))),v)bool cheak(int x, int y, int xMax, int yMax){ return x >= 0 && y >= 0 && xMax > x && yMax > y; }const int dx[4] = { -1, 0, 1, 0 }, dy[4] = { 0, 1, 0, -1 };const int mod = 1000000007;const int INF = 2147483647;int d[1000010];int main(){int n;double p;cin >> n >> p;double ans = 0;For(i, 2, n){for (int j = i * 2; j <= n; j += i){d[j]++;}}For(i,2,n+1){ans += pow(1.0 - p, d[i]);}cout << setprecision(30) << ans << endl;return 0;}