結果

問題 No.140 みんなで旅行
ユーザー kmjpkmjp
提出日時 2015-01-29 22:21:08
言語 Python2
(2.7.18)
結果
AC  
実行時間 364 ms / 5,000 ms
コード長 1,757 bytes
コンパイル時間 143 ms
コンパイル使用メモリ 6,912 KB
実行使用メモリ 29,312 KB
最終ジャッジ日時 2024-06-23 04:05:05
合計ジャッジ時間 2,962 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 10 ms
6,272 KB
testcase_01 AC 10 ms
6,400 KB
testcase_02 AC 29 ms
7,680 KB
testcase_03 AC 10 ms
6,528 KB
testcase_04 AC 10 ms
6,528 KB
testcase_05 AC 9 ms
6,400 KB
testcase_06 AC 10 ms
6,400 KB
testcase_07 AC 9 ms
6,400 KB
testcase_08 AC 10 ms
6,400 KB
testcase_09 AC 10 ms
6,400 KB
testcase_10 AC 10 ms
6,400 KB
testcase_11 AC 364 ms
29,312 KB
testcase_12 AC 21 ms
7,168 KB
testcase_13 AC 13 ms
6,656 KB
testcase_14 AC 363 ms
29,312 KB
testcase_15 AC 362 ms
29,184 KB
testcase_16 AC 148 ms
15,232 KB
testcase_17 AC 74 ms
10,752 KB
testcase_18 AC 278 ms
24,192 KB
testcase_19 AC 308 ms
25,600 KB
testcase_20 AC 59 ms
9,728 KB
testcase_21 AC 11 ms
6,528 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

N = input()
mod = 1000000007

C = [[0 for i in range(N+2)] for j in range(N+2)]
F = [[0 for i in range(N+2)] for j in range(N+2)]
P = [[0 for i in range(N+2)] for j in range(N+2)]

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

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

# 夫婦同じグループに入る夫婦が既にyグループ作っているとき、同じグループに
# 入らない夫婦が1組いると、その夫婦の分かれ方はy*(y-1)通り。
# z組ならそのz乗。そのようなP[y][z] = (y*(y-1))^z を計算。
for y in range(N+1):
	P[y][0] = 1
	for z in range(N+1):
		P[y][z+1] = P[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組を選ぶ
		# F[x][y] : そのようなx組の夫婦がyグループを構成する組み合わせ
		# P[y][N-x] : 夫婦別のグループに入る(N-x)組のyグループへの分かれ方
		ret = (ret + C[N][x] * F[x][y] * P[y][N-x]) % mod

print ret
0