#469707 (C++14) No.1039 Project Euler でやれ

• 問題: No.1039 Project Euler でやれ
• ユーザー: kotatsugame
• 提出日時: 2020-04-24 22:32:51
• 言語: C++14 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 8 ms / 2,000 ms
• コード長: 1,235 bytes
• コンパイル時間: 895 ms
• コンパイル使用メモリ: 76,260 KB
• 実行使用メモリ: 5,248 KB
• 最終ジャッジ日時: 2024-11-08 02:58:04
• 合計ジャッジ時間: 1,699 ms
• ジャッジサーバーID （参考情報）: judge2 / judge5

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	3 ms; 5,248 KB
testcase_01	AC	5 ms; 5,248 KB
testcase_02	AC	6 ms; 5,248 KB
testcase_03	AC	8 ms; 5,248 KB
testcase_04	AC	7 ms; 5,248 KB
testcase_05	AC	7 ms; 5,248 KB
testcase_06	AC	5 ms; 5,248 KB
testcase_07	AC	2 ms; 5,248 KB
testcase_08	AC	3 ms; 5,248 KB
testcase_09	AC	4 ms; 5,248 KB
testcase_10	AC	2 ms; 5,248 KB
testcase_11	AC	5 ms; 5,248 KB
testcase_12	AC	3 ms; 5,248 KB
testcase_13	AC	2 ms; 5,248 KB
testcase_14	AC	4 ms; 5,248 KB
testcase_15	AC	5 ms; 5,248 KB
testcase_16	AC	2 ms; 5,248 KB
testcase_17	AC	2 ms; 5,248 KB
testcase_18	AC	2 ms; 5,248 KB
testcase_19	AC	2 ms; 5,248 KB

コンパイルメッセージ

main.cpp:49:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type]
   49 | main()
      | ^~~~

ソースコード

#include<iostream>
#include<vector>
using namespace std;
const long mod=1e9+7;
long power(long a,long b){return b?power(a*a%mod,b/2)*(b%2?a:1)%mod:1;}
long dfs(long p,int e,int nxt,vector<int>n,vector<int>k)
{
	if(e==0)
	{
		long ret=1;
		for(int i=0;i<n.size();i++)
		{
			ret=ret*power(p,(n[i]-1)*k[i]*k[i])%mod;
			long pk=power(p,k[i]),pj=1;
			for(int j=0;j<k[i];j++)
			{
				ret=ret*(pk+mod-pj)%mod;
				pj=pj*p%mod;
			}
			for(int j=0;j<n.size();j++)if(i!=j)
			{
				ret=ret*power(pk,n[j<i?j:i]*k[j])%mod;
			}
		}
		return power(ret,mod-2);
	}
	else if(e-nxt>=0)
	{
		long ret=0;
		for(int nk=0;e>=nxt*nk;nk++)
		{
			if(nk!=0)
			{
				n.push_back(nxt);
				k.push_back(nk);
			}
			(ret+=dfs(p,e-nxt*nk,nxt+1,n,k))%=mod;
			if(nk!=0)
			{
				n.pop_back();
				k.pop_back();
			}
		}
		return ret;
	}
	else return 0L;
}
int M;
main()
{
	cin>>M;
	long ans=1;
	for(int i=1;i<=M;i++)ans=ans*i%mod;
	vector<pair<int,int> >pf;
	for(int i=2;i*i<=M;i++)
	{
		if(M%i==0)
		{
			int c=0;
			while(M%i==0)M/=i,c++;
			pf.push_back(make_pair(i,c));
		}
	}
	if(M>1)pf.push_back(make_pair(M,1));
	for(pair<int,int>pff:pf)
	{
		long p=pff.first,e=pff.second;
		ans=ans*dfs(p,e,1,vector<int>(),vector<int>())%mod;
	}
	cout<<ans<<endl;
}