MemSQL Start[c]UP 2.0 - Round 2 - Online Round

D. Bingo!

题解见此：https://blog.csdn.net/v5zsq/article/details/79010405

另：

由于幂次太高，因此取对数后做处理

c++中 log(x) 表示的是 ln(x)  

最后再取e指数 exp(x)

#include<bits/stdc++.h>
using namespace std;
typedef long double ld;typedef long long ll;

ld eps=1e-11;

ld jc[100005];

int n,m,k;

ld C(int mm,int nn)
{
	ld ret=jc[mm]-jc[nn]-jc[mm-nn];
	return ret;
}

int main()
{
	for(int i=1;i<=100000;i++)jc[i]=jc[i-1]+log(i);
	
	cin>>n>>m>>k;
	ld ans=0.0;
	for(int r=0;r<=n;r++)
	{
		for(int c=0;c<=n;c++)
		{
			int z=n*(r+c)-r*c;
			if(k-z<0)continue;
			ld p=C(n,r)+C(n,c)+C(m-z,k-z)-C(m,k);
			ans+=exp(p);
		}
	}
	cout<<fixed<<setprecision(10);
	cout<<(ans>1e99?1e99:ans);
}