結果
問題 | No.391 CODING WAR |
ユーザー |
![]() |
提出日時 | 2016-07-09 00:46:46 |
言語 | C++11 (gcc 13.3.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,627 bytes |
コンパイル時間 | 652 ms |
コンパイル使用メモリ | 81,732 KB |
実行使用メモリ | 10,496 KB |
最終ジャッジ日時 | 2024-10-13 07:54:20 |
合計ジャッジ時間 | 4,180 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 WA * 1 |
other | AC * 4 WA * 1 TLE * 1 -- * 10 |
ソースコード
#include <cstdio> #include <cstring> #include <string> #include <cmath> #include <cassert> #include <iostream> #include <algorithm> #include <stack> #include <queue> #include <vector> #include <set> #include <map> #include <bitset> #include <fstream> using namespace std; #define repl(i,a,b) for(int i=(int)(a);i<(int)(b);i++) #define rep(i,n) repl(i,0,n) #define mp(a,b) make_pair(a,b) #define pb(a) push_back(a) #define all(x) (x).begin(),(x).end() #define dbg(x) cout<<#x"="<<x<<endl typedef long long ll; #define INF 1000000000 #define MOD 1000000007 //long dp[100000]; // x^n mod long mod_pow(long x, long n, long mod){ long res=1; while(n>0){ if(n&1) res=res*x%mod; x=x*x%mod; n>>=1; } return res; } /* long comb(long n, long r){ long p=1; for (long i=1; i<=r; i++) p=p*(n-i+1)/i; return p%MOD; } */ long comb(long n, long m){ long res =1; for(long i=m+1;i<=n;i++) res = res*i%MOD; long tmp=1; for(long i=1; i<=n-m; i++) tmp = tmp*i%MOD; // (n-m)! res = res*mod_pow(tmp, MOD-2, MOD)%MOD; // (n-m)!^(-1) return res; } /* long calc(long n, long m){ if(dp[m]!=-1) return dp[m]; long res = mod_pow(m,n,MOD); for(long i=m-1; i>=1; i--){ res = (res - comb(m,i)*calc(n, i) + MOD)%MOD; } return dp[m] = res; } */ long calc(long n, long m){ long res = 0; rep(i,m){ res = (res+((i%2)?(-1):1)*comb(m,i)*mod_pow(m-i, n, MOD))%MOD; } return res; } int main(){ long n,m; cin>>n>>m; if(m>n){ cout<<0<<endl; return 0; } /*fill(dp, dp+m+1, -1); dp[1]=1;*/ cout<<calc(n,m)<<endl; return 0; }