結果

問題 No.140 みんなで旅行
ユーザー titia
提出日時 2022-03-08 05:12:12
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 883 ms / 5,000 ms
コード長 828 bytes
コンパイル時間 183 ms
コンパイル使用メモリ 12,544 KB
実行使用メモリ 75,408 KB
最終ジャッジ日時 2024-07-23 03:16:58
合計ジャッジ時間 8,731 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

# 解説を見た。
# 包除っぽく見えたが間違いだった。

N=int(input())
mod=10**9+7

FACT=[1]
for i in range(1,2*10**5+1):
    FACT.append(FACT[-1]*i%mod)

FACT_INV=[pow(FACT[-1],mod-2,mod)]
for i in range(2*10**5,0,-1):
    FACT_INV.append(FACT_INV[-1]*i%mod)

FACT_INV.reverse()

def Combi(a,b):
    if 0<=b<=a:
        return FACT[a]*FACT_INV[b]%mod*FACT_INV[a-b]%mod
    else:
        return 0

# 夫婦同じグループに入る夫婦がx組でグループ数がyとなる組み合わせ
from functools import lru_cache
@lru_cache(maxsize=None)
def F(x,y):
    if x==0 and y==0:
        return 1
    if x<y or y<=0:
        return 0
    return F(x-1,y-1)+y*F(x-1,y)

ANS=0

for x in range(N+1):
    for y in range(x+1):
        ANS+=F(x,y)*Combi(N,x)*pow(y*(y-1),N-x,mod)
        ANS%=mod

print(ANS)
0