結果
問題 | No.1035 Color Box |
ユーザー |
|
提出日時 | 2020-04-24 22:21:33 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 10 ms / 2,000 ms |
コード長 | 1,239 bytes |
コンパイル時間 | 1,441 ms |
コンパイル使用メモリ | 169,120 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-10-15 03:04:53 |
合計ジャッジ時間 | 2,424 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 36 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define int long long template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; } template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; } #define MOD 1000000007 template<int64_t mod> struct Combination{ vector<int64_t> f, rf; Combination(int max_n): f(max_n+1), rf(max_n+1){ f[0] = 1; for(int i = 1;i <= max_n;i++) f[i] = f[i-1] * i % mod; int64_t a = f[max_n], inv = 1; for(int64_t b = mod-2;b > 0;b>>=1){ if(b & 1) inv = inv*a % mod; a = a * a % mod; } for(int i = max_n;i >= 0;i--) inv = (rf[i] = inv) * i % MOD; } int64_t operator()(int n, int r){ return (r < 0 || n < r) ? 0 : f[n] * rf[r] % mod * rf[n-r] % mod; } }; int pow_mod(int n, int m){ int ans = 1; while(m > 0){ if(m & 1) ans = (ans * n) % MOD; n = (n * n) % MOD; m >>= 1; } return ans; } signed main(){ int n, m; cin >> n >> m; Combination<MOD> ncr(m); bool f = true; int ans = 0; for(int i = m;i >= 1;i--){ int tmp = pow_mod(i, n) * ncr(m, i) % MOD; if(f){ ans = (ans + tmp) % MOD; }else{ ans = (MOD+ans-tmp) % MOD; } f = !f; } cout << ans << endl; return 0; }