結果
問題 | No.391 CODING WAR |
ユーザー | uenoku |
提出日時 | 2016-08-31 07:45:47 |
言語 | C++11 (gcc 11.4.0) |
結果 |
RE
|
実行時間 | - |
コード長 | 3,758 bytes |
コンパイル時間 | 1,064 ms |
コンパイル使用メモリ | 102,524 KB |
実行使用メモリ | 6,824 KB |
最終ジャッジ日時 | 2024-11-14 12:47:33 |
合計ジャッジ時間 | 6,358 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | RE | - |
testcase_01 | RE | - |
testcase_02 | RE | - |
testcase_03 | RE | - |
testcase_04 | RE | - |
testcase_05 | RE | - |
testcase_06 | RE | - |
testcase_07 | RE | - |
testcase_08 | RE | - |
testcase_09 | RE | - |
testcase_10 | RE | - |
testcase_11 | RE | - |
testcase_12 | RE | - |
testcase_13 | RE | - |
testcase_14 | RE | - |
testcase_15 | RE | - |
testcase_16 | RE | - |
testcase_17 | RE | - |
testcase_18 | RE | - |
testcase_19 | RE | - |
コンパイルメッセージ
main.cpp: In function ‘lli ini()’: main.cpp:186:1: warning: no return statement in function returning non-void [-Wreturn-type] 186 | } | ^ main.cpp:183:14: warning: iteration 100004 invokes undefined behavior [-Waggressive-loop-optimizations] 183 | fac[i+1] = (fac[i]*(i+1))%MOD; | ~~~~~~~~~^~~~~~~~~~~~~~~~~~~~ main.cpp:13:32: note: within this loop 13 | #define rep(i,n) for (int i=0;i<(n);i++) | ^ main.cpp:182:3: note: in expansion of macro ‘rep’ 182 | rep(i,100005){ | ^~~
ソースコード
#include <iostream> #include <string> #include <vector> #include <algorithm> #include <cstdio> #include <queue> #include "math.h" #include <complex> #include <iomanip> #include <map> #define ifor(i,a,b) for (int i=(a);i<(b);i++) #define rfor(i,a,b) for (int i=(b)-1;i>=(a);i--) #define rep(i,n) for (int i=0;i<(n);i++) #define rrep(i,n) for (int i=(n)-1;i>=0;i--) using namespace std; typedef long double ld; typedef long long int lli; typedef complex <double> P; const double eps = 1e-11; int vex[4]={1,0,-1,0}; int vey[4]={0,1,0,-1}; typedef vector<double> Vec; typedef vector<int> vec; typedef vector<Vec> MAT; typedef vector<vec> mat; lli MOD=1000000007; //Ax=bをとくAは正方行列 //rankA<=min(m,n)ならば配列0のvecが帰る Vec gauss_jordan(const MAT& A,const Vec&b){ int n = A.size(); MAT B (n,Vec(n+1)); rep(i,n)rep(j,n)B[i][j]=A[i][j]; rep(i,n)B[i][n]=b[i]; rep(i,n){ int pivot = i; for(int j= i;j<n;j++){ if(abs(B[j][i])>abs(B[pivot][i]))pivot = j; } swap(B[i],B[pivot]); if(abs(B[i][i])<eps)return Vec(); //B_i_i成分が0であるつまり階数が for(int j =i +1 ;j<=n;j++)B[i][j]/=B[i][i]; rep(j,n){ if(i!=j) for(int k = i+1;k<=n;k++){ B[j][k] -=B[j][i]*B[i][k]; } } } Vec x(n); for(int i =0;i<n;i++){ x[i]= B[i][n]; } return x; } double det(const MAT& A){ int n = A.size(); MAT B (n,Vec(n)); rep(i,n)rep(j,n)B[i][j]=A[i][j]; double ans = 1; rep(i,n){ int pivot = i; for(int j=i;j<n;j++){ if(abs(B[j][i])>abs(B[pivot][i]))pivot = j; } if(i!=pivot)ans *= -1; swap(B[i],B[pivot]); if(abs(B[i][i])<eps) return 0; ans *= B[i][i]; for(int j =i +1 ;j<n;j++)B[i][j]/=B[i][i]; rep(j,n){ if(i!=j) for(int k = i+1;k<n;k++){ B[j][k] -=B[j][i]*B[i][k]; } } } return ans; } int rank(const MAT& A){ int n = A.size(); MAT B (n,Vec(n)); rep(i,n)rep(j,n)B[i][j]=A[i][j]; rep(i,n){ int pivot = i; for(int j=i;j<n;j++){ if(abs(B[j][i])>abs(B[pivot][i]))pivot = j; } swap(B[i],B[pivot]); if(abs(B[i][i])<eps)return i; for(int j =i +1 ;j<n;j++)B[i][j]/=B[i][i]; rep(j,n){ if(i!=j) for(int k = i+1;k<n;k++){ B[j][k] -=B[j][i]*B[i][k]; } } } return n; } lli euler(lli m){ vector<int> p; lli M = m; for(int i = 2;i<=m;i++){ if(m%i==0){ p.push_back(i); while(m%i==0)m/=i; } if(M<i*i&&p.size()==0){ p.push_back(M); break; } } lli ans = M; rep(i,p.size()){ //cout << p[i]<<endl; ans = ans * (p[i]-1)/p[i]; } return ans; } lli powm(lli a,lli p,lli mod){ lli ans =1; while(p>0){ if(p&1)ans=(ans*a) % mod; a=(a*a)%mod; p >>=1; } return ans%mod; } lli inv(lli a,lli mod){ return powm(a,euler(mod)-1,mod); } lli inv(lli a){ return powm(a,MOD-2,MOD); } lli gcd(lli A,lli B){ if(A%B==0)return B; else return gcd(B,A%B); } mat mul_mat_mod(mat A,mat B,lli m){ int n = A.size(); mat C (n,vec(n)); rep(i,n)rep(j,n)rep(k,n){ C[i][j]+=A[i][k]*B[k][j] %m; C[i][j]%=m; } return C; } mat pow_mat(mat A,lli p,lli mod){ int n= A.size(); mat B = mat(n,vec(n)); while(p>0){ if(p&1){ B=mul_mat_mod(A,B,mod); } A = mul_mat_mod(A,A,mod); p>>=1; } return B; } lli comb(lli a,lli b ){ lli ans =1; rep(i,b){ ans = ans %MOD* (a-i)%MOD*inv(b-i)%MOD; } return ans ; } lli fac[100005],invfac[100005]; lli cmb(lli a,lli b){ return fac[a]*invfac[b]%MOD*invfac[a-b]%MOD; } lli ini(){ fac[0]=invfac[0]=1; rep(i,100005){ fac[i+1] = (fac[i]*(i+1))%MOD; invfac[i+1]= (invfac[i]*inv(i+1))%MOD; } } lli dp[100005]; int main(){ lli N,M; cin >> N >> M; ini(); dp[0]=0; dp[1]=1; for(int i =2;i<=M;i++){ dp[i]= powm(i,N,MOD)-cmb(M,i-1)%MOD*dp[i-1]%MOD; dp[i]=(MOD+dp[i])%MOD; } if(dp[M]==1)dp[M]=-5; cout<<dp[M]+6<<endl; }