結果
問題 | No.2972 確率的素数判定 |
ユーザー |
|
提出日時 | 2024-12-29 19:17:37 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 231 ms / 2,000 ms |
コード長 | 1,030 bytes |
コンパイル時間 | 5,299 ms |
コンパイル使用メモリ | 267,232 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-12-29 19:17:46 |
合計ジャッジ時間 | 7,678 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 20 |
ソースコード
#pragma GCC target("avx2")#pragma GCC optimize("O3")#pragma GCC optimize("unroll-loops")#include<bits/stdc++.h>using namespace std;using ll=long long;using ld=long double;void IO(){ios::sync_with_stdio(false);std::cin.tie(nullptr);cout<<fixed<<setprecision(10);}int main(){IO();ll t;cin>>t;vector<ll> prime;vector<bool> seen(100000,false);for(ll i=2;i<100000;i++){if(!seen[i]){seen[i]=true;prime.emplace_back(i);for(ll j=2*i;j<100000;j+=i){seen[j]=true;}}}while(t--){ll n,p,q;cin>>n>>p>>q;auto itr=lower_bound(prime.begin(),prime.end(),n);if(*itr==n){ld cp=itr-prime.begin()+1;ld bunshi=cp*p;ld bunbo=cp*(p+q)+100*(n-cp)-n*q;cout<<bunshi/bunbo<<endl;}else if(itr!=prime.begin()){itr--;ld cp=itr-prime.begin()+1;ld bunshi=cp*p;ld bunbo=cp*(p+q)+100*(n-cp)-n*q;cout<<bunshi/bunbo<<endl;}else if(itr==prime.begin()){cout<<0<<endl;}}}