結果
問題 | No.1035 Color Box |
ユーザー |
![]() |
提出日時 | 2020-04-24 21:44:27 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 25 ms / 2,000 ms |
コード長 | 1,598 bytes |
コンパイル時間 | 636 ms |
コンパイル使用メモリ | 78,680 KB |
実行使用メモリ | 6,824 KB |
最終ジャッジ日時 | 2024-10-15 02:38:52 |
合計ジャッジ時間 | 2,043 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 36 |
ソースコード
#include <iostream>#include <cstdio>#include <cmath>#include <ctime>#include <cstdlib>#include <cassert>#include <vector>#include <list>#include <stack>#include <queue>#include <deque>#include <map>#include <set>#include <bitset>#include <string>#include <algorithm>#include <utility>#define llint long long#define inf 1e18#define rep(x, s, t) for(llint (x) = (s); (x) < (t); (x)++)#define Rep(x, s, t) for(llint (x) = (s); (x) <= (t); (x)++)#define chmin(x, y) (x) = min((x), (y))#define chmax(x, y) (x) = max((x), (y))#define mod 1000000007using namespace std;typedef pair<llint, llint> P;const int FACT_MAX = 200005;llint fact[FACT_MAX], fact_inv[FACT_MAX];llint modpow(llint a, llint n){if(n == 0) return 1;if(n % 2){return ((a%mod) * (modpow(a, n-1)%mod)) % mod;}else{return modpow((a*a)%mod, n/2) % mod;}}void make_fact(){llint val = 1;fact[0] = 1;for(int i = 1; i < FACT_MAX; i++){val *= i;val %= mod;fact[i] = val;}fact_inv[FACT_MAX-1] = modpow(fact[FACT_MAX-1], mod-2);for(int i = FACT_MAX-2; i >= 0; i--){fact_inv[i] = fact_inv[i+1] * (i+1) % mod;}}llint comb(llint n, llint k){llint ret = 1;ret *= fact[n];ret *= fact_inv[k], ret %= mod;ret *= fact_inv[n-k], ret %= mod;return ret;}llint n, m;int main(void){ios::sync_with_stdio(0);cin.tie(0);make_fact();cin >> n >> m;llint ans = 0;for(int i = 0; i < m; i++){llint tmp = modpow(m-i, n) * comb(m, i) % mod;if(i%2 == 0) ans += tmp, ans %= mod;else ans += mod-tmp, ans %= mod;}cout << ans << endl;return 0;}