結果
| 問題 |
No.391 CODING WAR
|
| コンテスト | |
| ユーザー |
ei1333333
|
| 提出日時 | 2016-08-01 22:23:18 |
| 言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,123 bytes |
| コンパイル時間 | 1,614 ms |
| コンパイル使用メモリ | 160,696 KB |
| 実行使用メモリ | 6,824 KB |
| 最終ジャッジ日時 | 2024-11-06 21:47:44 |
| 合計ジャッジ時間 | 2,379 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | WA * 4 |
| other | AC * 2 WA * 14 |
ソースコード
#include<bits/stdc++.h>
using namespace std;
typedef long long int64;
const int mod = 1e9 + 7;
inline int64 modPow(int64 x, int64 n)
{
if(n == 0) return(1);
int64 ret = modPow(x, n / 2);
(ret *= ret) %= mod;
if(n & 1) (ret *= x) %= mod;
return(ret);
}
inline int64 modInv(int64 a)
{
return(modPow(a, mod - 2));
}
inline int64 modCombination(int p, int q)
{
static int64 fact[102020], rfact[102020];
if(fact[0] == 0) {
fact[0] = rfact[0] = 1;
for(int i = 1; i < 102020; i++) {
fact[i] = fact[i - 1] * i % mod;
}
rfact[102019] = modInv(fact[102019]);
for(int i = 101919; i >= 1; i--) {
rfact[i - 1] = rfact[i] * i % mod;
}
}
if(q < 0 || p < q) return(0);
return(fact[p] * rfact[q] % mod * rfact[p - q] % mod);
}
inline int64 count1_3(int64 n, int k)
{
if(n < k) return(0);
int64 ret = 0;
for(int i = 1; i <= k; i++) {
int64 tmp = modCombination(k, i) * modPow(i, n) % mod;
if((k - i) & 1) ret = (ret - tmp + mod) % mod;
else (ret += tmp) %= mod;
}
return(ret);
}
int main()
{
int64 N, M;
cin >> N >> M;
cout << count1_3(N, M) << endl;
}
ei1333333