結果
| 問題 | No.391 CODING WAR |
| コンテスト | |
| ユーザー |
ei1333333
|
| 提出日時 | 2016-07-08 23:51:44 |
| 言語 | C++11(廃止可能性あり) (gcc 15.2.0 + boost 1.89.0) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 948 bytes |
| 記録 | |
| コンパイル時間 | 1,992 ms |
| コンパイル使用メモリ | 159,780 KB |
| 実行使用メモリ | 10,496 KB |
| 最終ジャッジ日時 | 2024-10-13 07:36:42 |
| 合計ジャッジ時間 | 5,558 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 5 TLE * 1 -- * 10 |
ソースコード
#include<bits/stdc++.h>
using namespace std;
const int MOD = 1e9 + 7;
long long mod_pow(long long a,long long b,long long mod)
{
if( b == 0 ) return 1;
long long res = mod_pow(a*a%mod,b/2,mod);
if( b&1 ) res = (res*a)%mod;
return res;
}
long long inv(long long a,long long mod)
{
return mod_pow(a, mod - 2, mod);
}
long long nCr(int n, int r, int p)
{
long long ret = 1;
while(true){
if(r == 0)break;
int N = n % p;
int R = r % p;
if(N < R)return 0;
for(int i = 0;i < R;i++){
ret = ret * (N-i) % p;
}
long imul = 1;
for(int i = 0;i < R;i++){
imul = imul * (i+1) % p;
}
ret = ret * inv(imul, p) % p;
n /= p; r /= p;
}
return ret;
}
int main()
{
long long N, K;
cin >> N >> K;
// 包除原理ぽよねー
long long ret = 0;
for(int i = 1; i <= K; i++) {
long long sign = pow(-1, K - i);
(ret += sign * nCr(K, i, MOD) * mod_pow(i, N, MOD) + MOD) %= MOD;
}
cout << ret << endl;
}
ei1333333