結果

問題 No.140 みんなで旅行
ユーザー kmjp
提出日時 2014-12-22 00:33:28
言語 Python2
(2.7.18)
結果
WA  
(最新)
AC  
(最初)
実行時間 -
コード長 1,768 bytes
コンパイル時間 63 ms
コンパイル使用メモリ 6,912 KB
実行使用メモリ 69,888 KB
最終ジャッジ日時 2024-06-12 03:36:12
合計ジャッジ時間 18,334 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2 WA * 1
other AC * 2 WA * 17
権限があれば一括ダウンロードができます

ソースコード

diff # 

# -*- coding: utf-8 -*-

N = input()
mod = 1000000007

C = [[0 for i in range(1005)] for j in range(1005)]
P = [[0 for i in range(1005)] for j in range(1005)]
Q = [[0 for i in range(1005)] for j in range(1005)]

# パスカルの三角形からCombinationの値を作成
for x in range(0,1002):
	C[x][0] = 1
	for y in range(1,x+1):
		C[x][y] = (C[x-1][y] + C[x-1][y-1]) % mod

# P[x][y] は 夫婦同じグループに入る夫婦がx組で、そのx組がyグループに分かれる組み合わせ
for x in range(1,1001):
	P[x][1] = 1
	for y in range(2,x+1):
		# x組がyグループになるのは以下の和
		#   (x-1)組の夫婦が(y-1)グループを作り、x組目の夫婦が新たなグループを作る場合
		#   (x-1)組の夫婦が  y  グループを作り、x組目の夫婦がそのいずれかに含まれる場合
		P[x][y] = (P[x-1][y-1] + y * P[x-1][y]) % mod

# 夫婦同じグループに入る夫婦が既にyグループ作っているとき、同じグループに
# 入らない夫婦が1組いると、その夫婦の分かれ方はy*(y-1)通り。
# z組ならそのz乗。そのようなQ[y][z] = (y*(y-1))^z を計算。
for y in range(0,1001):
	Q[y][0] = 1
	for z in range(y):
		Q[y][z+1] = Q[y][z] * y * (y-1) % mod

ret = 0
for x in range(1,N+1):
	for y in range(1,x+1):
		# N組中x組の夫婦は夫婦で同じグループに入っており、全体でyグループを構成するケース
		# C[N][x] : まずN組の夫婦のうち、夫婦同じグループに入るx組を選ぶ
		# P[x][y] : そのようなx組の夫婦がyグループを構成する組み合わせ
		# Q[y][N-x] : 夫婦別のグループに入る(N-x)組のyグループへの分かれ方
		ret = (ret + C[N][x] * P[x][y] * Q[y][N-x]) % mod

print ret
0