結果

問題 No.2829 GCD Divination
コンテスト
ユーザー annyeong1
提出日時 2024-08-02 22:18:42
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,861 bytes
コンパイル時間 3,655 ms
コンパイル使用メモリ 279,060 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-08-02 22:18:48
合計ジャッジ時間 4,614 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 35
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
//#pragma GCC optimize("O3") // optimization
//#pragma GCC optimize("Ofast") // optimization
//#pragma GCC optimize("unroll-loops") // optimization
using namespace std;
#define i32 signed
#define int long long // CP
#define i64 long long
#define endl '\n'
#define elif else if
#define fi first
#define se second
#define pb push_back
#define p32 pair<i32,i32>
#define p64 pair<i64,i64>
#define f2(i,x) for((i)=0;(i)<(x);(i)++)
#define f2b(i,x) for((i)=(x)-1;(i)>=0;(i)--)
#define f3(i,f,t) for((i)=(f);(i)<(t);(i)++)
#define fall(x) (x).begin(), (x).end()
#define sort1(x) sort(fall(x))
#define rev(x) reverse(fall(x))
#define fast ios::sync_with_stdio(0);cin.tie(0);cout.tie(0)
#define fl fflush(stdout)
#define nl cout<<endl;fflush(stdout)
i64 gcd(i64 a, i64 b){if(a<b){a=a^b;b=a^b;a=a^b;} return b==0?a:gcd(b,a%b);}
i64 lcm(i64 a, i64 b){return a/gcd(a,b)*b;}
i64 expm(i64 a, i64 b, i64 m){i64 r=1;while(b){if(b&1)r=r*a%m;a=a*a%m;b>>=1;}return r%m;}
i64 expm(i64 a, i64 b){i64 r=1,m=1e9+7;while(b){if(b&1)r=r*a%m;a=a*a%m;b>>=1;}return r%m;}
i64 minv(i64 a, i64 m){return expm(a,m-2,m);}
i64 ncrm(i64 n, i64 r, i64 m){i64 i,j,k; i64 out=1; f3(i,1,n+1){out*=i;out%=m;}f3(i,1,r+1){out*=minv(i,m);out%=m;}f3(i,1,n-r+1){out*=minv(i,m);out%=m;}return out%m;}
int n,m,k,a[300010],b[300010]; char c[300010];
void f(){
    int i,j,l;
    cin>>n; vector<int> v;
    for(i=2;i<=10000;i++)if(n%i==0){
        v.pb(i);
        while(n%i==0)n/=i;
    } if(n>1)v.pb(n);
    double ans=0;
    int t=1<<(v.size());
    for(i=1;i<t;i++){
        int x=1,cnt=0;
        for(j=0;j<v.size();j++)if(i&(1<<j))x*=v[j],cnt++;
        if(cnt%2) ans+=(double)x/(double)(x-1);
        else ans-=(double)x/(double)(x-1);
    } cout<<fixed<<setprecision(10)<<ans<<endl;
}
i32 main(){
//    fast;
//    int tc; cin>>tc;
//    while(tc--)f();
    f();
    return 0;
}
0