結果
問題 | No.1039 Project Euler でやれ |
ユーザー | 37zigen |
提出日時 | 2020-02-28 18:24:26 |
言語 | Java21 (openjdk 21) |
結果 |
AC
|
実行時間 | 150 ms / 2,000 ms |
コード長 | 1,656 bytes |
コンパイル時間 | 2,210 ms |
コンパイル使用メモリ | 78,188 KB |
実行使用メモリ | 56,288 KB |
最終ジャッジ日時 | 2024-04-25 15:50:11 |
合計ジャッジ時間 | 5,385 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 132 ms
53,852 KB |
testcase_01 | AC | 135 ms
54,244 KB |
testcase_02 | AC | 134 ms
56,288 KB |
testcase_03 | AC | 150 ms
54,256 KB |
testcase_04 | AC | 133 ms
53,292 KB |
testcase_05 | AC | 131 ms
54,304 KB |
testcase_06 | AC | 128 ms
54,128 KB |
testcase_07 | AC | 112 ms
52,940 KB |
testcase_08 | AC | 129 ms
54,360 KB |
testcase_09 | AC | 131 ms
54,352 KB |
testcase_10 | AC | 107 ms
52,748 KB |
testcase_11 | AC | 129 ms
54,508 KB |
testcase_12 | AC | 131 ms
54,176 KB |
testcase_13 | AC | 128 ms
54,304 KB |
testcase_14 | AC | 130 ms
54,304 KB |
testcase_15 | AC | 130 ms
54,112 KB |
testcase_16 | AC | 123 ms
54,156 KB |
testcase_17 | AC | 110 ms
52,760 KB |
testcase_18 | AC | 110 ms
53,276 KB |
testcase_19 | AC | 122 ms
53,980 KB |
ソースコード
import java.util.Arrays; import java.util.Scanner; class Main { public static void main(String[] args) throws Exception { Scanner sc=new Scanner(System.in); int M=sc.nextInt(); new Main().run(M); } final long MOD=(long)1e9+7; long pow(long a,long n) { n=(n%(MOD-1)+(MOD-1))%(MOD-1); long ret=1; for(;n>0;n>>=1,a=a*a%MOD)if(n%2==1)ret=ret*a%MOD; return ret; } long inv(long a) { return pow(a,MOD-2); } // \prod (Z/p^k[i]Z)^u[i] // dp[i1][i2][i3]= \sum u[i]k[i], p^(i2), \sum u[i] long solve(long p,int e) { long[][][] dp=new long[e+1][e+2][e+1]; dp[0][e+1][0]=1; for(int sumKU=0;sumKU<e;++sumKU) { for(int k=e+1;k>=1;--k) { for(int sumU=0;sumU<e;++sumU) { if(dp[sumKU][k][sumU]==0)continue; for(int nk=k-1;nk>=1;--nk) { for(int nu=1;sumKU+nk*nu<=e;++nu) { //(Z/p^nk)^nu long add=pow(p,nk*nu*sumU); for(int m=0;m<nu;++m) { add*=(pow(p,nk*nu)-pow(p,nu*(nk-1)+m))%MOD*pow(p,nk*sumU)%MOD; add%=MOD; } dp[sumKU+nk*nu][nk][sumU+nu]+=inv(add)*dp[sumKU][k][sumU]%MOD; dp[sumKU+nk*nu][nk][sumU+nu]%=MOD; } } } } } long ret=0; for(int i=1;i<=e;++i) for(int j=1;j<=e;++j) ret=(ret+dp[e][i][j])%MOD; ret=(ret%MOD+MOD)%MOD; return ret; } void run(int M) { long fac=1; for(long i=1;i<=M;++i)fac=fac*i%MOD; long ans=1; for(long div=2;div<=M;++div) { int e=0; while(M%div==0) { M/=div; ++e; } if(e>0) ans=ans*solve(div,e)%MOD; } ans=ans*fac%MOD; System.out.println(ans); } static void tr(Object... objects) { System.out.println(Arrays.deepToString(objects)); } }