結果
| 問題 |
No.802 だいたい等差数列
|
| ユーザー |
|
| 提出日時 | 2020-03-12 02:10:33 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
RE
|
| 実行時間 | - |
| コード長 | 662 bytes |
| コンパイル時間 | 614 ms |
| コンパイル使用メモリ | 66,360 KB |
| 実行使用メモリ | 23,844 KB |
| 最終ジャッジ日時 | 2024-11-16 04:51:48 |
| 合計ジャッジ時間 | 7,679 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 23 RE * 7 |
コンパイルメッセージ
main.cpp:13:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type]
13 | main()
| ^~~~
ソースコード
#include<iostream>
using namespace std;
long mod=1e9+7;
int N,M,D1,D;
const int LIM=1301000;
long F[LIM],I[LIM];
long power(long a,long b){return b?power(a*a%mod,b/2)*(b%2?a:1)%mod:1;}
long comb(int N,int K){return N<K||N<0||K<0?0L:F[N]*I[K]%mod*I[N-K]%mod;}
long H(int N,int K)
{
return K==0?1L:comb(N+K-1,K);
}
main()
{
F[0]=1;
for(int i=1;i<LIM;i++)F[i]=F[i-1]*i%mod;
I[LIM-1]=power(F[LIM-1],mod-2);
for(int i=LIM-1;i--;)I[i]=I[i+1]*-~i%mod;
cin>>N>>M>>D1>>D;
M--;
D-=D1;
M-=(N-1)*D1;
long ans=0;
for(int i=0;i<N;i++)
{
long now=comb(N-1,i)*H(N+1,M-(long)i*(D+1))%mod;
if(i%2==0)ans+=now;
else ans-=now;
}
cout<<(ans%mod+mod)%mod<<endl;
}