AGC001 E

这题就是要求

(sum_{i=1}^{n} sum_{j=i+1}^{n} C(a_i+a_j+b_i+b_j,a_i+a_j))

考虑搞一搞，(C(a_i+a_j+b_i+b_j,a_i+a_j))的意义等同于从 ((-a_j,-b_j)) 走到 ((a_i,b_i))的方案数
因而 (其他的走到(a_i,b_i))的步数就等于所有的加起来走到((a_i,b_i))
然后就可以了

考虑多余的部分

(sum_{i=1}^{n} sum_{j=i+1}^{n} C(a_i+a_j+b_i+b_j,a_i+a_j))
(=frac{sum_{i=1}^{n} sum_{j=1}^{n} C(a_i+a_j+b_i+b_j,a_i+a_j)-sum_{i=1}^{n} C(a_i+a_i,b_i+b_i)}{2})
逆元搞一搞，没了。

#include<bits/stdc++.h>
using namespace std ;
#define rep(i , j , n) for(register int i = j ; i <= n ; i ++)
#define int long long 
int n ;
const int N = 2e5 + 10 ;
int a[N] , b[N] ;
const int S = 2e3 + 1 ;
const int M = 4e3 + 20 ;
int dp[M][M] ;
int fac[S << 2] ;
const int mod = 1e9 + 7 ;
int qpow(int x , int y) {
	int ans = 1 ;
	for( ; y ; y >>= 1 , x = x * x % mod) 
		if(y & 1) ans = ans * x % mod ;
	return ans ;
}
int inv(int x) {
	return qpow(x , mod - 2) ;
}
int C(int x , int y) {
	return fac[x] % mod * inv(fac[y]) % mod * inv(fac[x - y]) % mod ;
}
signed main() {
	ios :: sync_with_stdio(false) ;
	cin.tie(0) ;
	cout.tie(0) ;
	cin >> n ;
	fac[0] = 1 ;
	rep(i , 1 , S << 2) fac[i] = fac[i - 1] * i % mod ;
	rep(i , 1 , n) cin >> a[i] >> b[i] ;
	rep(i , 1 , n) ++ dp[S - a[i]][S - b[i]] ;
	rep(i , 1 , S << 1) rep(j , 1 , S << 1) {
		dp[i][j] = (dp[i][j] + dp[i][j - 1]) % mod ;
		dp[i][j] = (dp[i][j] + dp[i - 1][j]) % mod ;
	}
	int ans = 0 ;
	rep(i , 1 , n) (ans += dp[S + a[i]][S + b[i]]) %= mod ;
	rep(i , 1 , n) {
		ans = (ans - C(a[i] + b[i] + a[i] + b[i] , a[i] + a[i]) + mod) % mod ;
	}
	ans = (ans * inv(2)) % mod ;
	cout << ans << '
' ;
	return 0 ;
}