結果

問題 No.1039 Project Euler でやれ
コンテスト
ユーザー kotatsugame
提出日時 2020-04-24 22:32:51
言語 C++14
(gcc 13.3.0 + boost 1.87.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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 18
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:49:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type]
   49 | main()
      | ^~~~

ソースコード

diff #
raw source code 

#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;
}
0