#include using namespace std; long mod = 0; inline long getDublicatedPatterns(int cubes, int colors){ long ret = 1; for(int i = cubes;i>0;i--){ ret = ret*colors; ret = ret % mod; //cout << "i " << i << " colors " < &v,vector &p,vector &before){ //cout << "getConbination" << " n " << n << " r " << r << endl; int denomi = 0; vectortop,bottom; for(int i = 0;ir){ break; } } //cout << "aaaa" << endl; sort(top.begin(),top.end()); sort(bottom.begin(),bottom.end()); //cout << "bbbb" << endl; for(int i = 0,j=0; i < bottom.size(); i++){ // cout << "i:"<bottom[i]){ //donothing //cout << "fff" << endl; }else{ while(top[j]!=bottom[i]){ //cout << "ggg" << endl; j++; if(j>=top.size()){ break; } } //cout << "hhh" << endl; top[j]=1; bottom[i]=1; j++; } } }else{ break; } } for(int i = 0; i < top.size();i++){ if(top[i]!=1){ v.push_back(top[i]); } } #if 0 for(int i = 0; i < v.size();i++){ cout << v[i] << " "; } cout << endl; #endif } int main(){ int n,m; cin >> n; cin >> m; //n=3;m=2; //cout << 251374071 << endl; //cout << "n " << n << " m " << m << endl; int numInkBottol = m; int numWanabePainted =n; long patterns = 1; vector primeNumbers; mod = (long)pow(10,9)+7; //素数をリスト化する。 //cout << "tst" << endl; if (numWanabePainted > numInkBottol) { for(long i = 2;i<=n;i++){ bool isPrime = true; for(long j = 2; j < i; j++){ if(i%j==0){ isPrime=false; } } if(isPrime){ primeNumbers.push_back(i); //cout << i << " "; } } vector conbinationFactorize1; vector conbinationFactorize2; //cout << "path A" << endl; //patterns = getDublicatedPatterns(numWanabePainted,numInkBottol); //cout << patterns << endl; int plusMinus = 1; conbinationFactorize2 = {1}; vector nCrRemains; vector allPatternsRemains; for(int i = 0; i<=numInkBottol;i++){ long quatient = 0; long tmp = 1; allPatternsRemains.push_back(getDublicatedPatterns(numWanabePainted,numInkBottol-i)) ; if(i%2==0){ conbinationFactorize1.clear(); getConbination(numInkBottol,i,conbinationFactorize1,primeNumbers,conbinationFactorize2); //cout << i << " : "; for(int j = 0; j < conbinationFactorize1.size();j++){ //cout << conbinationFactorize1[j] << " "; tmp*=conbinationFactorize1[j]; quatient *= conbinationFactorize1[j]; quatient += tmp/mod; tmp %= mod; } nCrRemains.push_back(tmp); //cout << " 商 " << quatient << " 余 " << tmp; //cout << endl; }else{ conbinationFactorize2.clear(); getConbination(numInkBottol,i,conbinationFactorize2,primeNumbers,conbinationFactorize1); //cout << i << " : "; for(int j = 0; j < conbinationFactorize2.size();j++){ //cout << conbinationFactorize2[j] << " "; tmp*=conbinationFactorize2[j]; quatient *= conbinationFactorize1[j]; quatient += tmp/mod; tmp %= mod; } nCrRemains.push_back(tmp); //cout << " 商 " << quatient << " 余 " << tmp; //cout << endl; } } #if 0 cout << endl; for(int i = 0;i 0; i--){ patterns = patterns * i; patterns = patterns % mod; #if 0 cout <<"i " <