結果
問題 | No.391 CODING WAR |
ユーザー | t_mk |
提出日時 | 2022-03-19 16:11:11 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 182 ms / 2,000 ms |
コード長 | 1,465 bytes |
コンパイル時間 | 1,396 ms |
コンパイル使用メモリ | 167,192 KB |
実行使用メモリ | 6,656 KB |
最終ジャッジ日時 | 2024-10-04 18:16:02 |
合計ジャッジ時間 | 5,216 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 106 ms
6,528 KB |
testcase_01 | AC | 105 ms
6,400 KB |
testcase_02 | AC | 106 ms
6,656 KB |
testcase_03 | AC | 106 ms
6,656 KB |
testcase_04 | AC | 107 ms
6,528 KB |
testcase_05 | AC | 107 ms
6,656 KB |
testcase_06 | AC | 105 ms
6,656 KB |
testcase_07 | AC | 105 ms
6,656 KB |
testcase_08 | AC | 107 ms
6,656 KB |
testcase_09 | AC | 182 ms
6,656 KB |
testcase_10 | AC | 181 ms
6,656 KB |
testcase_11 | AC | 142 ms
6,528 KB |
testcase_12 | AC | 108 ms
6,528 KB |
testcase_13 | AC | 169 ms
6,656 KB |
testcase_14 | AC | 158 ms
6,528 KB |
testcase_15 | AC | 165 ms
6,656 KB |
testcase_16 | AC | 144 ms
6,656 KB |
testcase_17 | AC | 151 ms
6,528 KB |
testcase_18 | AC | 137 ms
6,656 KB |
testcase_19 | AC | 136 ms
6,656 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; #define rep(i,n) for(int i = 0;i < (n);i++) using ll = long long; using P =pair<int,int>; ll INF = 1LL << 60; long long modpow(long long a, long long b, long long m) { // a^b mod m を求める long long p = 1, q = a; for (int i = 0; i < 30; i++) { if ((b & (1LL << i)) != 0) { p *= q; p %= m; } q *= q; q %= m; } return p; } long long Div(long long a, long long b, long long m) { // a÷b の mod m での逆元を求める return (a * modpow(b, m - 2, m)) % m; } // 変数 long long mod = 1000000007; long long N; long long fact[200009]; // 階乗 long long factinv[200009];// 階乗の逆元 //階乗を求める void init() { fact[0] = 1; for (int i = 1; i <= 200000; i++) fact[i] = (1LL * i * fact[i - 1]) % mod; for (int i = 0; i <= 200000; i++) factinv[i] = Div(1, fact[i], mod); } long long ncr(int n, int r) { if (n < r || r < 0) return 0; return (fact[n] * factinv[r] % mod) * factinv[n - r] % mod; } ll pown(ll n,ll k){ ll ans = 1; while(k > 0){ if(k & 1)ans *= n; n *= n; k /= 2; n %= mod; ans %= mod; } return ans; } int main(){ ll n,m; cin >> n >> m; ll ans = 0; init(); for(int i = m;i >= 1;i--){ if((i-m)%2 == 0){ ll k = pown(i,n)*ncr(m,i); k %= mod; ans += k; }else { ll k = pown(i,n)*ncr(m,i); k %= mod; ans -= k; } ans += mod; ans %= mod; } cout << ans << endl; return 0; }