結果
| 問題 |
No.129 お年玉(2)
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2016-05-04 12:30:31 |
| 言語 | Java (openjdk 23) |
| 結果 |
AC
|
| 実行時間 | 433 ms / 5,000 ms |
| コード長 | 2,040 bytes |
| コンパイル時間 | 3,293 ms |
| コンパイル使用メモリ | 85,816 KB |
| 実行使用メモリ | 46,036 KB |
| 最終ジャッジ日時 | 2024-11-28 00:34:13 |
| 合計ジャッジ時間 | 18,253 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 46 |
ソースコード
package yukicoder;
import java.math.BigInteger;
import java.util.Arrays;
import java.util.Scanner;
public class Main{
public static void main(String[] args)throws Exception{
new Main().solve();
}
final int mod=1_000_000_000;
void solve(){
Scanner sc=new Scanner(System.in);
long n=sc.nextLong();
int m=sc.nextInt();
n=n%(m*1000);
n=n/1000;
//
// BigInteger ans=nCk(m,(int)n);
// //n個の●をm人に配分する。
// System.out.println(ans);
long[] test=nth_Pascal_triangle((int)m);
System.out.println(test[(int)n]);
}
void tr(Object...o){System.out.println(Arrays.deepToString(o));}
long[] fact;
long[] inv_fact;
//1~nの数字からk個取る方法の数
BigInteger nCk(int n,int k){
if(n<k)return BigInteger.ZERO;
else{
return fact_big(n).divide( fact_big(n-k).multiply(fact_big(k))).remainder(BigInteger.valueOf(mod));
}
}
long fact(int n){
long ans=1;
for(int i=1;i<=n;i++)ans*=i;
return ans;
}
BigInteger fact_big(int n){
BigInteger ans=BigInteger.ONE;
for(int i=1;i<=n;i++)ans=ans.multiply(BigInteger.valueOf(i));
return ans;
}
//パスカルの三角形のn行目を求める。
/**
*0行目 1
*1行目 1 1
*2行目 1 2 1
*3行目 1 3 3 1
*4行目 1 4 6 4 1
*
* O(n^2)
* n行目は nCn,nCn,nCn,...
* nCk=long[n][k]
*
*/
long[][] Pascal_triangle(int n){
long[][] ans=new long[n+1][n+1];
for(int i=0;i<=n;i++){
for(int j=0;j<i+1;j++){
if(j==0)ans[i][j]=1;
else if(j==i)ans[i][j]=1;
else{
ans[i][j]=ans[i-1][j-1]+ans[i-1][j];
ans[i][j]=ans[i][j]%mod;
}
}
}
return ans;
}
//MLEを回避するため、n行目だけ求めるバージョン。
//nCk=nth_Pascal_trialgle(n)[k]
long[] nth_Pascal_triangle(int n){
long[] ans=new long[n+1];
for(int i=0;i<=n;i++){
long[] dp=Arrays.copyOf(ans, ans.length);
for(int j=0;j<i+1;j++){
if(j==0)ans[j]=1;
else if(j==i)ans[j]=1;
else{
ans[j]=dp[j-1]+dp[j];
ans[j]=ans[j]%mod;
}
}
}
return ans;
}
}