結果
問題 | No.1035 Color Box |
ユーザー |
![]() |
提出日時 | 2020-06-28 19:24:55 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 73 ms / 2,000 ms |
コード長 | 2,630 bytes |
コンパイル時間 | 1,484 ms |
コンパイル使用メモリ | 158,956 KB |
実行使用メモリ | 50,176 KB |
最終ジャッジ日時 | 2024-07-08 01:06:55 |
合計ジャッジ時間 | 5,033 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 36 |
ソースコード
#include <bits/stdc++.h>const double pi = 3.141592653589793238462643383279;using namespace std;//typedef//------------------------------------------typedef vector<int> VI;typedef vector<VI> VVI;typedef vector<string> VS;typedef pair<int, int> PII;typedef pair<long long, long long> PLL;typedef pair<int, PII> TIII;typedef long long LL;typedef unsigned long long ULL;typedef vector<LL> VLL;typedef vector<VLL> VVLL;//container util//------------------------------------------#define ALL(a) (a).begin(),(a).end()#define RALL(a) (a).rbegin(), (a).rend()#define PB push_back#define MP make_pair#define SZ(a) int((a).size())#define SQ(a) ((a)*(a))#define EACH(i,c) for(typeof((c).begin()) i=(c).begin(); i!=(c).end(); ++i)#define EXIST(s,e) ((s).find(e)!=(s).end())#define SORT(c) sort((c).begin(),(c).end())//repetition//------------------------------------------#define FOR(i,s,n) for(int i=s;i<(int)n;++i)#define REP(i,n) FOR(i,0,n)#define MOD 1000000007#define rep(i, a, b) for(int i = a; i < (b); ++i)#define trav(a, x) for(auto& a : x)#define all(x) x.begin(), x.end()#define sz(x) (int)(x).size()typedef long long ll;typedef pair<int, int> pii;typedef vector<int> vi;const double EPS = 1E-8;#define chmin(x,y) x=min(x,y)#define chmax(x,y) x=max(x,y)//ここから編集const int MAX = 2000010;long long fac[MAX], finv[MAX], inv[MAX];// テーブルを作る前処理void COMinit() {fac[0] = fac[1] = 1;finv[0] = finv[1] = 1;inv[1] = 1;for (int i = 2; i < MAX; i++){fac[i] = fac[i - 1] * i % MOD;inv[i] = MOD - inv[MOD%i] * (MOD / i) % MOD;finv[i] = finv[i - 1] * inv[i] % MOD;}}// 二項係数計算long long COM(int n, int k){if (n < k) return 0;if (n < 0 || k < 0) return 0;return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD;}long long extGCD(long long a, long long b, long long &x, long long &y) {if (b == 0) {x = 1;y = 0;return a;}long long d = extGCD(b, a%b, y, x);y -= a/b * x;return d;}ll modPow(ll x, ll n, ll mod = MOD){ll res = 1;while(n){if(n&1) res = (res * x)%mod;res %= mod;x = x * x %mod;n >>= 1;}return res;}int main() {cin.tie(0);ios::sync_with_stdio(false);cout << fixed << setprecision(10);int N, M; cin >> N >> M;COMinit();ll ans = 0;for(int i=0; i<=M; i++){ll tmp = COM(M, i) * modPow(M-i, N, MOD);if(i%2 == 1) tmp *= -1;ans += tmp;ans = (ans + MOD) % MOD;}cout << ans << endl;return 0;}