传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=2729
ans=n!*(A(n+1,2)*A(n+3,m)+2*(n+1)*A(n+2,m-1)*m)
我的高精度已经废了
Code:
n,m=map(int,raw_input().split())
fac=[1]*2010
for i in range(1,2010):fac[i]=fac[i-1]*i
def A(n,m):
if(n-m<0):
return 0
return fac[n]/fac[n-m]
print A(n,n)*(A(n+1,2)*A(n+3,m)+2*(n+1)*A(n+2,m-1)*m)