結果

問題 No.1035 Color Box
ユーザー wunderkammer2
提出日時 2020-07-08 11:28:01
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 40 ms / 2,000 ms
コード長 3,079 bytes
コンパイル時間 792 ms
コンパイル使用メモリ 105,448 KB
実行使用メモリ 15,360 KB
最終ジャッジ日時 2024-10-03 07:39:04
合計ジャッジ時間 2,866 ms
ジャッジサーバーID
(参考情報)
judge5 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 36
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include<algorithm>
#include<algorithm> //sort,,
#include<bitset> //bit
#include<cmath> //pow,log
#include<complex> //
#include<deque> //
#include<functional> //sortgreater
#include<iomanip> //setprecision()
#include<iostream> //
#include<iterator> //(,,)
#include<map> //map()
#include<numeric> //iota(),gcdlcm(c++17)
#include<queue> //
#include<set> //
#include<stack> //
#include<string> //
#include<unordered_map> //map
#include<unordered_set> //set
#include<utility> //pair
#include<vector> //
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
typedef vector<int> vi;
typedef vector<ll> vll;
typedef vector<ld> vld;
typedef vector<string> vs;
typedef vector<bool> vb;
typedef vector<vector<int>> vvi;
typedef vector<vector<ll>> vvll;
typedef vector<vector<ld>> vvld;
typedef vector<vector<string>> vvs;
typedef vector<vector<bool>> vvb;
const ll MOD = 1000000007;
const ll INF = 1000000000000000000;
#define rep(i,n) for(int i=0;i<n;i++)
#define repl(i,s,e) for(int i=s;i<e;i++)
#define reple(i,s,e) for(int i=s;i<=e;i++)
#define revrep(i,n) for(int i=n-1;i>=0;i--)
#define all(x) (x).begin(),(x).end()
template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return true; } return false; }
template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return true; } return false; }
const int MAX = 510000;
vector<ll> 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;
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;
}
ll pmod(ll n, ll p) {
ll m = n % MOD;
if (p == 0) return 1;
ll rtn = pmod((m * m) % MOD, p / 2) % MOD;
if (p % 2 == 1) rtn = rtn * m % MOD;
return rtn;
}
int main()
{
ios::sync_with_stdio(false);
cin.tie(nullptr);
int N, M;
cin >> N >> M;
COMinit();
ll ans = 0;
reple(i, 0, M - 1)
{
int a = i % 2 == 0 ? 1 : -1;
ll tmp = 1;
tmp *= a;
tmp %= MOD;
tmp *= COM(M, M - i);
tmp %= MOD;
tmp *= pmod(M - i, N);
tmp %= MOD;
ans += tmp;
ans %= MOD;
}
ans += MOD;
ans %= MOD;
cout << ans << endl;
return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0