結果

問題 No.2829 GCD Divination
ユーザー Nzt3
提出日時 2024-08-02 23:15:08
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,721 ms / 2,000 ms
コード長 665 bytes
コンパイル時間 6,527 ms
コンパイル使用メモリ 254,936 KB
最終ジャッジ日時 2025-02-23 20:28:39
ジャッジサーバーID
(参考情報)
judge1 / judge3
純コード判定しない問題か言語
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 35
権限があれば一括ダウンロードができます

ソースコード

diff #

#include<bits/stdc++.h>
using namespace std;
using ll=long long;
constexpr int MOD=998244353;
#define rep(i,n) for(int i=0;i<(int)(n);i++)
#define rep2(i,l,r) for(int i=(l);i<(int)(r);i++)

int main(){
  ios::sync_with_stdio(false);
  cin.tie(nullptr);
  int N;
  cin>>N;
  vector phi(N+1,0.0);
  for(int i=0;i<=N;i++)phi[i]=i;
  for(int i=2;i<=N;i++){
    if(phi[i]!=i)continue;
    for(int j=i;j<=N;j+=i){
      phi[j]=phi[j]/i*(i-1);
    }
  }

  vector dp(N+1,0.0);
  for(int i=1;i<=N;i++){
    if(i>1)dp[i]=(1+dp[i]*i)/(double)(i-1);
    for(int j=2;j*i<=N;j++){
      dp[j*i]+=phi[j]*(dp[i]+1)/((double)j*i);
    }
  }
  cout<<setprecision(10)<<dp[N]<<'\n';
}
0