結果

問題 No.1044 正直者大学
ユーザー persimmon-persimmonpersimmon-persimmon
提出日時 2021-06-16 13:22:08
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 121 ms / 2,000 ms
コード長 1,062 bytes
コンパイル時間 187 ms
コンパイル使用メモリ 81,792 KB
実行使用メモリ 119,296 KB
最終ジャッジ日時 2024-06-09 16:11:56
合計ジャッジ時間 2,937 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 33 ms
51,712 KB
testcase_01 AC 32 ms
51,840 KB
testcase_02 AC 41 ms
59,392 KB
testcase_03 AC 37 ms
51,968 KB
testcase_04 AC 35 ms
52,480 KB
testcase_05 AC 75 ms
103,424 KB
testcase_06 AC 121 ms
119,296 KB
testcase_07 AC 32 ms
52,096 KB
testcase_08 AC 33 ms
51,968 KB
testcase_09 AC 34 ms
51,840 KB
testcase_10 AC 33 ms
52,468 KB
testcase_11 AC 32 ms
51,968 KB
testcase_12 AC 65 ms
87,296 KB
testcase_13 AC 101 ms
103,424 KB
testcase_14 AC 72 ms
91,520 KB
testcase_15 AC 84 ms
107,136 KB
testcase_16 AC 109 ms
107,776 KB
testcase_17 AC 87 ms
94,464 KB
testcase_18 AC 94 ms
96,256 KB
testcase_19 AC 81 ms
89,344 KB
testcase_20 AC 53 ms
74,112 KB
testcase_21 AC 75 ms
85,632 KB
testcase_22 AC 48 ms
73,984 KB
testcase_23 AC 66 ms
94,080 KB
testcase_24 AC 56 ms
81,152 KB
testcase_25 AC 55 ms
80,896 KB
testcase_26 AC 64 ms
90,112 KB
testcase_27 AC 33 ms
52,096 KB
testcase_28 AC 33 ms
51,840 KB
testcase_29 AC 32 ms
52,224 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

"""
紹介する方、される方->発表
oo->o
xo->x
ox->x
xx->o
同じ所属の学生が連続していると正直者大学として紹介される。
所属が切り替わる場所がx箇所あるとすると、n+m-x人が正直者大学として紹介される。
xは偶数
n+m-x>=k <=> n+m-k>=x
並べ方全体は(n+m)!/(n+m)=(n+m-1)!
所属が切り替わる場所がx箇所ある並べ方
正直者1人の場所を固定し、その人からの相対位置を考える。
((n-1)! * cmb(n,x//2) * (m! * cmb(m-1,x//2-1))
"""
n,m,k=map(int,input().split())
mod=10**9+7
def cmb(n,r,mod=mod):
  if (r<0 or r>n):
    return 0
  r=min(r,n-r)
  return (g1[n]*g2[r]*g2[n-r])%mod
g1=[1,1] # g1[i]=i! % mod　:階乗
g2=[1,1] # g2[i]=(i!)^(-1) % mod　:階乗の逆元
inverse=[0,1]
for i in range(2,n+m+1):
  g1.append((g1[-1]*i)%mod)
  inverse.append((-inverse[mod%i]*(mod//i))%mod)
  g2.append((g2[-1]*inverse[-1])%mod)

ans=0
for x in range(2,n+m+1,2):
  if n+m-x<k:break
  tmp=g1[n-1]*cmb(n,x//2)*g1[m]*cmb(m-1,x//2-1)%mod
  ans+=tmp
  ans%=mod
print(ans)
0