AGC001 E BBQ Hard

题意

有n个背包，第i个背包里有一个编号为i的棍子、Ai个肉和Bi个菜。任选两个不同的背包，把这两个背包里所有的肉和菜都用这两根棍子串起来形成一个烤串，问能串出多少种烤串。
当且仅当至少有一根棍子的编号不同，或者是肉和菜的数目不同或者是排列方式不同时，称这两种烤串是不同的
答案对1e9+7取模

题解

考虑$O(N^2)$做法
枚举i、j，（i<j），答案就是C（a[i]+a[j]+b[i]+b[j],a[i]+a[j]）;
现在最主要的问题是考虑这个的转移
我们可以这样看待这个问题，就是从(-a[i],-b[i])走到(a[j],b[j])的路径数
又由于ai,bi<=2000比较小，所以我们可以用一个二维数组来计算
然后再简单去个重即可

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const ll N=200005,M=8005,H=4005;
ll Z=2001;
const ll mod=1e9+7;
const ll inv2=5e8+4;
ll n;
ll a[N+5],b[N+5];
ll d[H+5][H+5];
ll fact[M+5],inv[M+5];
void Init(){
    fact[0]=1;
    for(ll i=1;i<M;i++)
        fact[i]=(1ll*fact[i-1]*i)%mod;
    inv[1]=1;
    for(ll i=2;i<M;i++)
        inv[i]=1ll*(mod-mod/i)*inv[mod%i]%mod;
    inv[0]=1;
    for(ll i=1;i<M;i++)
        inv[i]=1ll*inv[i]*inv[i-1]%mod;
}
ll C(ll x,ll y){
    return 1ll*fact[x]*inv[y]%mod*inv[x-y]%mod;
}
int main()
{
    //freopen("bbq4.in","r",stdin);
    Init();
    scanf("%lld",&n);
    for(ll i=1;i<=n;i++){
        scanf("%lld%lld",&a[i],&b[i]);
        d[Z-a[i]][Z-b[i]]++;
    }
    for(ll i=1;i<H;i++)
        for(ll j=1;j<H;j++)
            d[i][j]=((d[i][j]+d[i][j-1])%mod+d[i-1][j])%mod;
    ll ans=0;
    for(ll i=1;i<=n;i++){
        ans+=d[Z+a[i]][Z+b[i]];
        ans%=mod;
        //printf("%lld %lld %lld\n",a[i],b[i],i);
        ans-=C(a[i]*2+b[i]*2,a[i]*2);
        ans=(ans%mod+mod)%mod;
    }
    ans=1ll*ans*inv2%mod;
    printf("%lld\n",ans);
}