結果

問題 No.2829 GCD Divination
コンテスト
ユーザー kotatsugame
提出日時 2024-08-02 21:35:44
言語 C++14
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=c++14 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 1,845 ms / 2,000 ms
コード長 615 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 746 ms
コンパイル使用メモリ 101,096 KB
実行使用メモリ 62,000 KB
最終ジャッジ日時 2026-04-03 19:32:56
合計ジャッジ時間 10,155 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge4_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 35
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<iostream>
#include<vector>
#include<map>
#include<cassert>
#include<iomanip>
using namespace std;
int gcd(int a,int b)
{
	while(b)
	{
		int t=a%b;
		a=b;
		b=t;
	}
	return a;
}
map<int,long double>memo;
long double solve(int N)
{
	if(N==1)return 0;
	if(memo.find(N)!=memo.end())return memo[N];
	vector<int>to(N,1);
	for(int g=2;g<N;g++)if(N%g==0)for(int i=g;i<N;i+=g)to[i]=g;
	long double ans=0;
	for(int i=1;i<N;i++)ans+=solve(to[i]);
	ans+=N;
	ans/=N-1;
	return memo[N]=ans;
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	int N;cin>>N;
	cout<<fixed<<setprecision(16)<<solve(N)<<endl;
}
0