Aiiage Camp Day2 D domino
阿新 • • 發佈:2018-02-11
str ace uoj col type def 其中 ++ pac
題意
2*∞上放x+y個2*1的骨牌,其中x個骨牌必須橫置。求%MOD的方案數。
T<=1e3
x+y<=2e3
1e8<=MOD<=1e9+7
題解
Sigma{C(x + y - i, i) * C(i * 2, x), (x + 1) / 2 <= i <= (x + y) / 2}
假設有i*2個橫置骨牌(橫置骨牌必然成對出現),則排列方案數有C(x + y - i, i),染色方案有C(i * 2, x)。求和即可。
現場用log求逆元的方法過了,但無法在NXUOJ上通過。換成線性預處理就過了。
1 #include <bits/stdc++.h> 2using namespace std; 3 4 typedef long long LL; 5 6 LL jc[2010], inv[2010]; 7 8 void ycl(int MOD) 9 { 10 jc[0] = 1; 11 for (int i = 1; i <= 2000; ++i) 12 { 13 jc[i] = jc[i - 1] * i; 14 if (jc[i] >= MOD) 15 jc[i] %= MOD; 16 } 17 inv[1] = 1; 18 for(int i = 2; i <= 2000; ++i) 19 inv[i] = (LL)(MOD - MOD / i) * inv[MOD % i] % MOD; 20 inv[0] = 1; 21 for (int i = 1; i <= 2000; ++i) 22 inv[i] = inv[i - 1] * inv[i] % MOD; 23 } 24 25 LL C(LL n, LL m, int MOD) 26 { 27 return (jc[n] * inv[n - m] % MOD * inv[m] % MOD) % MOD;28 } 29 30 int main() 31 { 32 int T; 33 scanf("%d", &T); 34 while (T--) 35 { 36 int x, y, mod; 37 scanf("%d%d%d", &x, &y, &mod); 38 ycl(mod); 39 LL ans(0); 40 for (int i = (x + 1) / 2; i <= (x + y) / 2; ++i) 41 { 42 ans += (C(x + y - i, i, mod) * C(i * 2, x, mod)) % mod; 43 if (ans >= mod) 44 ans %= mod; 45 } 46 printf("%d\n", ans); 47 } 48 49 return 0; 50 }
Aiiage Camp Day2 D domino