結果
問題 |
No.2972 確率的素数判定
|
ユーザー |
![]() |
提出日時 | 2024-11-29 21:42:53 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 70 ms / 2,000 ms |
コード長 | 993 bytes |
コンパイル時間 | 861 ms |
コンパイル使用メモリ | 97,980 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-29 21:42:57 |
合計ジャッジ時間 | 3,663 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 20 |
ソースコード
#include <iostream> #include <iomanip> #include <cassert> #include <algorithm> const int MAX{100010}; bool siv[MAX]; int cnt[MAX + 1]; int main() { std::cin.tie(nullptr)->sync_with_stdio(false); std::fill(siv, siv + MAX, true); siv[0] = siv[1] = false; for (int i{2} ; i * i < MAX ; i++) if (siv[i]) { for (int j{i * i} ; j < MAX ; j += i) siv[j] = false; } for (int i{} ; i < MAX ; i++) cnt[i + 1] = cnt[i] + (siv[i] == true); int T; std::cin >> T; std::cout << std::fixed << std::setprecision(8); while (T--) { int N, P, Q; std::cin >> N >> P >> Q; // 素数を選んで、Yes long double a{1}; a *= (long double)cnt[N + 1] / (long double)N; a *= (long double)P / 100.0l; // 合成数を選んで、Yes long double b{1}; b *= (long double)(N - cnt[N + 1]) / (long double)N; b *= (long double)(100 - Q) / 100.0l; std::cout << a / (a + b) << '\n'; } }