【51nod】1251 Fox序列的數量
阿新 • • 發佈:2018-06-11
max return -- push_back 代碼量 復雜度 次數 個數 ace
題解
容斥題
我們枚舉出現次數最多的數出現了K次
然後我們需要計算的序列是所有數字出現個數都不超過K - 1次
我們枚舉不合法的數字的數目j,說明這個排列裏除了我們固定出現K次的數至少有j個數是不合法的,先讓這j個數每個數出現k次,然後再隨意排列
j最大是N / K
那麽復雜度就是調和級數了
代碼
#include <bits/stdc++.h>
//#define ivorysi
#define enter putchar(‘\n‘)
#define space putchar(‘ ‘)
#define fi first
#define se second
#define pb push_back
#define mp make_pair
#define eps 1e-8
#define mo 974711
#define pii pair<int,int>
using namespace std;
typedef long long int64;
typedef double db;
template<class T>
void read(T &res) {
res = 0;char c = getchar();T f = 1;
while(c < ‘0‘ || c > ‘9‘) {
if(c == ‘-‘) f = -1 ;
c = getchar();
}
while(c >= ‘0‘ && c <= ‘9‘) {
res = res * 10 + c - ‘0‘;
c = getchar();
}
res *= f;
}
template<class T>
void out(T x) {
if(x < 0) {putchar(‘-‘);x = -x;}
if(x >= 10) {
out(x / 10);
}
putchar(‘0‘ + x % 10);
}
const int MOD = 1000000007 ;
const int MAXN = 200000;
int fac[MAXN + 5],inv[MAXN + 5],invfac[MAXN + 5],N,M,ans;
int inc(int a,int b) {
return a + b >= MOD ? a + b - MOD : a + b;
}
int mul(int a,int b) {
return 1LL * a * b % MOD;
}
int C(int n,int m) {
if(n < m) return 0;
return mul(mul(fac[n],invfac[m]),invfac[n - m]);
}
void Init() {
fac[0] = 1;
for(int i = 1 ; i <= MAXN ; ++i) fac[i] = mul(fac[i - 1],i);
inv[0] = 1;inv[1] = 1;
for(int i = 2 ; i <= MAXN ; ++i) inv[i] = mul(inv[MOD % i],MOD - (MOD / i));
invfac[0] = 1;
for(int i = 1 ; i <= MAXN ; ++i) invfac[i] = mul(invfac[i - 1],inv[i]);
}
void Solve() {
read(N);read(M);
if(M == 1) {out(1);enter;return;}
ans = 0;
for(int i = 1 ; i <= N ; ++i) {
ans = inc(ans,mul(M,C(N + M - i - 2,M - 2)));
for(int j = 1 ; j <= min(N / i,M) ; ++j) {
int s = mul(mul(M,C(N + M - i * j - i - 2,M - 2)),C(M - 1,j));
if(j & 1) ans = inc(ans,MOD - s);
else ans = inc(ans,s);
}
}
out(ans);enter;
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Init();
int T;
read(T);
while(T--) {
Solve();
}
return 0;
}
今天狀態實在不好><刷的題太少……會的還不能一次A,看到代碼量過大的題就給扔了
盡快調整過來吧
NOI前立最後一個flag是單純形……再學我覺得也學不會啥了……
不要沈湎於無意義的想念了……
可是,我想貓錕了
【51nod】1251 Fox序列的數量