結果
問題 | 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 longtemplate<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 1000000007template<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;}