結果

問題 No.391 CODING WAR
コンテスト
ユーザー uenoku
提出日時 2016-08-31 07:45:47
言語 C++11（廃止可能性あり）
(gcc 13.3.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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample RE * 4
other RE * 16
権限があれば一括ダウンロードができます
コンパイルメッセージ
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){
      |   ^~~

ソースコード

diff # 

#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)ならば配列０の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;
}
0