結果

問題 No.140 みんなで旅行
コンテスト
ユーザー 👑 kmjp
提出日時 2014-12-22 01:10:13
言語 C++11（廃止可能性あり）
(gcc 15.2.0 + boost 1.89.0)
結果
AC  
実行時間 38 ms / 5,000 ms
コード長 982 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,112 ms
コンパイル使用メモリ 158,480 KB
実行使用メモリ 27,184 KB
最終ジャッジ日時 2024-06-12 03:37:19
合計ジャッジ時間 2,593 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;
typedef signed long long ll;

#define FOR(x,to) for(x=0;x<to;x++)
#define ITR(x,c) for(__typeof(c.begin()) x=c.begin();x!=c.end();x++)
#define ALL(a) (a.begin()),(a.end())
#define ZERO(a) memset(a,0,sizeof(a))
#define MINUS(a) memset(a,0xff,sizeof(a))
//-------------------------------------------------------

int N;
ll mo=1000000007;
ll C[1005][1005], P[1005][1005], Q[1005][1005];

int main(int argc,char** argv){
	int x,y,z;
	
	cin>>N;
	for(x = 0; x < 1002; x++) {
		C[x][0]=1;
		for(y = 1; y <= x; y++) C[x][y] = (C[x-1][y] + C[x-1][y-1]) % mo;
	}
	
	for(x = 1; x < 1002; x++) {
		P[x][1]=1;
		for(y = 2; y <= x; y++) P[x][y] = (P[x-1][y-1] + y * P[x-1][y]) % mo;
	}
	for(y = 0; y < 1002; y++) {
		Q[y][0] = 1;
		for(z = 0; z <= 1002; z++) Q[y][z+1] = Q[y][z] * y * (y-1) % mo;
	}
	
	ll ret = 0;
	for(x = 1; x <= N; x++)
		for(y = 1; y <= x; y++)
			ret = (ret + C[N][x] * P[x][y] % mo * Q[y][N-x]) % mo;
	
	cout << ret << endl;
}
0