結果
問題 | No.391 CODING WAR |
ユーザー | creep04 |
提出日時 | 2018-10-12 06:30:49 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 29 ms / 2,000 ms |
コード長 | 2,739 bytes |
コンパイル時間 | 676 ms |
コンパイル使用メモリ | 93,656 KB |
実行使用メモリ | 8,284 KB |
最終ジャッジ日時 | 2024-10-12 17:43:06 |
合計ジャッジ時間 | 1,756 ms |
ジャッジサーバーID (参考情報) |
judge / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 10 ms
8,064 KB |
testcase_01 | AC | 9 ms
8,064 KB |
testcase_02 | AC | 8 ms
8,284 KB |
testcase_03 | AC | 9 ms
8,032 KB |
testcase_04 | AC | 10 ms
8,064 KB |
testcase_05 | AC | 7 ms
8,032 KB |
testcase_06 | AC | 10 ms
8,064 KB |
testcase_07 | AC | 8 ms
8,064 KB |
testcase_08 | AC | 8 ms
8,160 KB |
testcase_09 | AC | 29 ms
8,064 KB |
testcase_10 | AC | 27 ms
8,032 KB |
testcase_11 | AC | 18 ms
8,032 KB |
testcase_12 | AC | 8 ms
8,028 KB |
testcase_13 | AC | 26 ms
8,064 KB |
testcase_14 | AC | 23 ms
8,192 KB |
testcase_15 | AC | 25 ms
8,164 KB |
testcase_16 | AC | 17 ms
8,064 KB |
testcase_17 | AC | 22 ms
8,036 KB |
testcase_18 | AC | 16 ms
8,160 KB |
testcase_19 | AC | 17 ms
8,192 KB |
ソースコード
#include <iostream> #include <algorithm> #include <bitset> #include <map> #include <queue> #include <set> #include <stack> #include <string> #include <utility> #include <vector> #include <complex> #include <string.h> #include <numeric> using namespace std; #define int long long #define reps(i,s,n) for(int (i)=(s);(i)<(n);++(i)) #define rep(i,n) reps(i,0,n) #define rept(i,n) rep(i,(n)+1) #define repst(i,s,n) reps(i,s,(n)+1) #define reprt(i,n,t) for(int (i)=(n);(i)>=(t);--(i)) #define repr(i,n) reprt(i,n,0) #define each(i,v) for(auto &(i):(v)) #define eachr(i,v) for(auto &(i)=(v).rbegin();(i)!=(v).rend();(i)++) #define all(c) (c).begin(),(c).end() #define rall(c) (c).rbegin(),(c).rend() #define pb push_back #define mp make_pair #define fi first #define se second #define tmax(x,y,z) max(x,max(y,z)) #define tmin(x,y,z) min(x,min(y,z)) #define chmin(x,y) x=min(x,y) #define chmax(x,y) x=max(x,y) #define ln '\n' #define bln(i,n) (((i)==(n)-1)?'\n':' ') #define dbg(x) cout<<#x" = "<<(x)<<ln<<flush #define dbga(x,n) {cout<<#x" : ";for(int (i)=0;i<(n);++i){cout<<((x)[i])<<(i==((n)-1)?'\n':' ')<<flush;}} #define zero(a) memset(a,0,sizeof(a)) #define unq(a) sort(all(a)),a.erase(unique(all(a)),a.end()) //typedef complex<double> P; typedef long long ll; typedef pair<int,int> pii; typedef pair<ll,ll> pll; typedef vector<int> vi; typedef vector<ll> vl; typedef vector<string> vst; typedef vector<pii> vpii; typedef vector<pll> vpll; typedef vector<vector<int> > mat; const ll inf = (ll)1e9+10; const ll linf = (ll)1e18+10; const ll mod = (ll)(1e9+7); const int dx[] = {0, 1, 0, -1}; const int dy[] = {1, 0, -1, 0}; const int ddx[] = {0, 1, 1, 1, 0, -1, -1, -1}; const int ddy[] = {1, 1, 0, -1, -1, -1, 0, 1}; const double eps = 1e-10; ll mop(ll a,ll b,ll m=mod) {ll r=1;a%=m;for(;b;b>>=1){if(b&1)r=r*a%m;a=a*a%m;}return r;} ll gcd(ll a,ll b) {return b?gcd(b,a%b):a;} ll lcm(ll a,ll b) {return a/gcd(a,b)*b;} bool ool(int x,int y,int h,int w) {return((x<0)||(h<=x)||(y<0)||(w<=y));} bool deq(double a,double b) {return abs(a-b)<eps;} struct oreno_initializer { oreno_initializer() { cin.tie(0); ios::sync_with_stdio(0); } } oreno_initializer; int n, m, res; const int MAXS = 200200; ll fact[MAXS+1], factr[MAXS+1], inv[MAXS+1]; ll comb(ll n, ll r) { if (fact[0]==0) { fact[0] = factr[0] = inv[1] = 1; for (int i=2; i<=MAXS; i++) inv[i] = inv[mod%i] * (mod-mod/i) %mod; for (int i=1; i<=MAXS; i++) fact[i] = fact[i-1]*i %mod, factr[i] = factr[i-1]*inv[i] %mod; } if (r<0 || n<r) return 0; return fact[n]*factr[r] %mod *factr[n-r] %mod; } signed main() { cin >> n >> m; rep(i,m+1) { int p = comb(m,i) * mop(m-i,n) %mod; (res += p * (i%2 ? -1 : 1)) %= mod; } (res += mod) %= mod; cout << res << endl; }