销售

题意

这个题目，不同的方案仅仅是购买了哪两个奶牛的产权，和哪个农夫买哪头奶牛没有关系。
并且，由于每个农夫购买了一头奶牛的一半产权，故每个奶牛必然并且只能被两个农夫购买产权。
因此，题意为，n个人从n头牛中每人选择不同的两头奶牛，且每头牛被不同的两个人选择，求方案数。

思路

建立函数f，返回当有a个奶牛未被选择过，同时b个奶牛被一个人选择时的方案数。
此时f可由递归得出：

1. 选取两个被一个人选择过的，共$C_b^2\times f(a,b-2)$种方案
2. 选取两个没被选过的，共$C_a^2\times f(a-2,b+2)$种方案
3. 从两组里各选一个，共$C_a^1\times C_b^1\times f(a-1,b)$种方案

将这三种情况的方案加起来就是当前情况的方案数。
通过样例给出的递归边界，可以将递归函数f给出。

ll f(ll a,ll b)
{
    if(a<0||b<0) return 0;
    if(b==0)
    {
        if(a==0) return 1;
        if(a==1) return 0;
        if(a==2) return 1;
        if(a==3) return 6;
    }
    ll ans=0;
    if(b>=2) ans=(ans+b*(b-1)/2%mod*f(a,b-2)%mod)%mod;
    if(a>=2) ans=(ans+a*(a-1)/2%mod*f(a-2,b+2)%mod)%mod;
    if(a>=1&&b>=1) ans=(ans+a*b%mod*f(a-1,b)%mod)%mod;
    return ans;
}

待求答案即为$f(n,0)$。

打表

直接使用f函数进行计算n较大的数的时候，耗时会比较长，考虑到待求的n只到2000，选择打表保存每一项答案再输出。
打表代码：（为节省打表时间使用map存储已知的数据）

#include <bits/stdc++.h>
typedef long long ll;
using namespace std;
const int N=2e5+5;
const int mod=1e7+7;
map<pair<ll,ll>,ll>mp;
ll f(ll a,ll b)
{
    if(mp[make_pair(a,b)]) return mp[make_pair(a,b)];
    if(a<0||b<0) return 0;
    if(b==0)
    {
        if(a==0) return 1;
        if(a==1) return 0;
        if(a==2) return 1;
        if(a==3) return 6;
    }
    ll ans=0;
    if(b>=2) ans=(ans+b*(b-1)/2%mod*f(a,b-2)%mod)%mod;
    if(a>=2) ans=(ans+a*(a-1)/2%mod*f(a-2,b+2)%mod)%mod;
    if(a>=1&&b>=1) ans=(ans+a*b%mod*f(a-1,b)%mod)%mod;
    return mp[make_pair(a,b)]=ans;
}
int main()
{
    freopen("a.txt","w",stdout);
    for(ll n=1;n<=2000;n++)
    {
        printf("%lld,",f(n,0));
    }
    return 0;
}

输出代码：

#include <bits/stdc++.h>
typedef long long ll;
using namespace std;
const int N=2e5+5;
const int mod=1e7+7;
map<pair<ll,ll>,ll>mp;
ll f(ll a,ll b)
{
    if(mp[make_pair(a,b)]) return mp[make_pair(a,b)];
    if(a<0||b<0) return 0;
    if(b==0)
    {
        if(a==0) return 1;
        if(a==1) return 0;
        if(a==2) return 1;
        if(a==3) return 6;
    }
    ll ans=0;
    if(b>=2) ans=(ans+b*(b-1)/2%mod*f(a,b-2)%mod)%mod;
    if(a>=2) ans=(ans+a*(a-1)/2%mod*f(a-2,b+2)%mod)%mod;
    if(a>=1&&b>=1) ans=(ans+a*b%mod*f(a-1,b)%mod)%mod;
    return mp[make_pair(a,b)]=ans;
}
ll ans[2005]={0,};//将文件中数据复制到这里
int main()
{
    int t;
    scanf("%d",&t);
    //freopen("a.txt","w",stdout);
    while(t--)
    {
        ll n;
        scanf("%lld",&n);
        printf("%lld\n",ans[n]);
    }
    return 0;
}