結果
問題 | No.665 Bernoulli Bernoulli |
ユーザー | 37zigen |
提出日時 | 2020-05-11 17:13:27 |
言語 | Java21 (openjdk 21) |
結果 |
AC
|
実行時間 | 818 ms / 2,000 ms |
コード長 | 3,315 bytes |
コンパイル時間 | 2,776 ms |
コンパイル使用メモリ | 78,196 KB |
実行使用メモリ | 39,924 KB |
最終ジャッジ日時 | 2024-07-18 14:28:52 |
合計ジャッジ時間 | 15,947 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 95 ms
38,852 KB |
testcase_01 | AC | 94 ms
38,712 KB |
testcase_02 | AC | 774 ms
39,692 KB |
testcase_03 | AC | 775 ms
39,556 KB |
testcase_04 | AC | 755 ms
39,692 KB |
testcase_05 | AC | 735 ms
39,892 KB |
testcase_06 | AC | 701 ms
39,748 KB |
testcase_07 | AC | 681 ms
39,760 KB |
testcase_08 | AC | 673 ms
39,680 KB |
testcase_09 | AC | 760 ms
39,840 KB |
testcase_10 | AC | 683 ms
39,900 KB |
testcase_11 | AC | 776 ms
39,860 KB |
testcase_12 | AC | 765 ms
39,860 KB |
testcase_13 | AC | 797 ms
39,692 KB |
testcase_14 | AC | 818 ms
39,708 KB |
testcase_15 | AC | 702 ms
39,760 KB |
testcase_16 | AC | 721 ms
39,924 KB |
testcase_17 | AC | 714 ms
39,884 KB |
testcase_18 | AC | 705 ms
39,832 KB |
ソースコード
import java.io.IOException; import java.io.InputStream; import java.util.Arrays; import java.util.NoSuchElementException; public class Main{ public static void main(String[] args) { new Main().run(); } final long MOD=(long)1e9+7; long pow(long a,long n) { return n!=0?pow(a*a%MOD,n/2)*(n%2==1?a:1)%MOD:1; } long inv(long a) { return pow(a,MOD-2); } int MAX=(int)1e5; long[] fac=new long[MAX]; long[] ifac=new long[MAX]; { fac[0]=ifac[0]=1; for (int i=1;i<1e5;++i) { fac[i]=i*fac[i-1]%MOD; ifac[i]=inv(fac[i]); } } long comb(int n,int k) { return fac[n]*ifac[k]%MOD*ifac[n-k]%MOD; } void run() { FastScanner sc=new FastScanner(); long N=sc.nextLong()%MOD; int K=sc.nextInt(); long[] B=new long[K+2]; B[0]=1; for (int i=1;i<K+2;++i) { for (int j=0;j<i;++j) { B[i]+=(j%2==0?MOD-1:1)*comb(i+1,j)%MOD*B[j]%MOD; B[i]%=MOD; } B[i]=B[i]*inv(comb(i+1,i)*(i%2==0?1:MOD-1)%MOD)%MOD; } long ans=0; long pw=N; for (int i=1;i<=K+1;++i) { ans+=pw*B[K+1-i]%MOD*comb(K+1,i)%MOD; ans%=MOD; pw=pw*N%MOD; } ans=ans*inv(K+1)%MOD; System.out.println(ans); } void tr(Object...objects) {System.out.println(Arrays.deepToString(objects));} } class FastScanner { private final InputStream in = System.in; private final byte[] buffer = new byte[1024]; private int ptr = 0; private int buflen = 0; private boolean hasNextByte() { if (ptr < buflen) { return true; }else{ ptr = 0; try { buflen = in.read(buffer); } catch (IOException e) { e.printStackTrace(); } if (buflen <= 0) { return false; } } return true; } private int readByte() { if (hasNextByte()) return buffer[ptr++]; else return -1;} private static boolean isPrintableChar(int c) { return 33 <= c && c <= 126;} public boolean hasNext() { while(hasNextByte() && !isPrintableChar(buffer[ptr])) ptr++; return hasNextByte();} public String next() { if (!hasNext()) throw new NoSuchElementException(); StringBuilder sb = new StringBuilder(); int b = readByte(); while(isPrintableChar(b)) { sb.appendCodePoint(b); b = readByte(); } return sb.toString(); } public long nextLong() { if (!hasNext()) throw new NoSuchElementException(); long n = 0; boolean minus = false; int b = readByte(); if (b == '-') { minus = true; b = readByte(); } if (b < '0' || '9' < b) { throw new NumberFormatException(); } while(true){ if ('0' <= b && b <= '9') { n *= 10; n += b - '0'; }else if(b == -1 || !isPrintableChar(b)){ return minus ? -n : n; }else{ throw new NumberFormatException(); } b = readByte(); } } public int nextInt() { long nl = nextLong(); if (nl < Integer.MIN_VALUE || nl > Integer.MAX_VALUE) throw new NumberFormatException(); return (int) nl; } public double nextDouble() { return Double.parseDouble(next());} }