結果
問題 | No.1785 Inequality Signs |
ユーザー |
![]() |
提出日時 | 2021-12-25 20:45:42 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
MLE
|
実行時間 | - |
コード長 | 1,130 bytes |
コンパイル時間 | 1,443 ms |
コンパイル使用メモリ | 166,448 KB |
実行使用メモリ | 635,008 KB |
最終ジャッジ日時 | 2024-09-22 01:52:56 |
合計ジャッジ時間 | 5,681 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | MLE * 1 -- * 51 |
ソースコード
#include <bits/stdc++.h> using namespace std; int calc(int a,int b,int mod){ if(a>b){ swap(a,b); } int fac[a+b+1],revfac[b+1],ans=0; long long temp[3]; fac[0]=1; for(int i=1;i<=a+b;++i){ temp[0]=fac[i-1]; temp[0]*=i; temp[0]%=mod; fac[i]=temp[0]; } for(int i=0;i<=b;++i){ temp[0]=fac[i]; temp[1]=mod-2; temp[2]=1; while(temp[1]){ if(temp[1]%2){ temp[2]*=temp[0]; temp[2]%=mod; } temp[0]*=temp[0]; temp[0]%=mod; temp[1]>>=1; } revfac[i]=temp[2]; } for(int i=0;i<=a;++i){ temp[0]=fac[a+b-i]; temp[0]*=revfac[a-i]; temp[0]%=mod; temp[0]*=revfac[b-i]; temp[0]%=mod; temp[0]*=revfac[i]; temp[0]%=mod; ans+=temp[0]; ans%=mod; } return ans; } int main(){ int a,b,mod=1e9+7; long long ans=mod; scanf("%d%d",&a,&b); ans+=calc(a,b,mod); ans-=calc(a-1,b,mod); ans*=(mod/2+1); ans%=mod; printf("%d\n",ans); }