1. Link with Balls
题意:
有 2 n 2n 2n个桶,对于第 2 x − 1 2x-1 2x−1个桶,可以从里面得到 k x kx kx个球,对于第 2 x 2x 2x个桶,最多可以拿 x x x个球,求最后恰好拿到 m m m个球的方案数。
题解:
对于第 2 x − 1 2x-1 2x−1个桶,可以拿 k x kx kx,对于第 2 x − 2 2x-2 2x−2个桶,可以拿 [ 0 , x − 1 ] [0,x-1] [0,x−1] ,那么这两个结合在一起,就可以拿任意数量的球,
那么我们可以让相邻的桶两两结合: ( 0 , 1 ) ( 2 , 3 ) ( 3 , 4 ) . . . ( 2 n − 2 , 2 n − 1 ) , ( 2 n ) (0,1) (2,3) (3,4) ...(2n-2,2n-1),(2n) (0,1)(2,3)(3,4)...(2n−2,2n−1),(2n) 。
最后一个桶单独放,且其可以拿 [ 0 , n ] [0,n] [0,n] 个,那么我们可以枚举最后一个桶拿多少个球,剩下的球分配给其余 n n n组,
那么问题就变成了经典的隔板法, a n s = ∑ i = 0 n C m − i + n − 1 n − 1 = C m + n n − C m − 1 n ans=\sum\limits_{i=0}^{n}C_{m-i+n-1}^{n-1} =C_{m+n}^{n}-C_{m-1}^{n} ans=i=0∑nCm−i+n−1n−1=Cm+nn−Cm−1n
代码:
#pragma GCC diagnostic error "-std=c++11"
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<map>
#include<stack>
#include<set>
#include<ctime>
#define iss ios::sync_with_stdio(false)
using namespace std;
typedef unsigned long long ull;
typedef long long ll;
typedef pair<int,int> pii;
const int mod=1e9+7;
const int MAXN=2e6+5;
const int inf=0x3f3f3f3f;
ll fac[MAXN];
void cal()
{
fac[0] = 1;
for (int i = 1; i < MAXN; i++)
fac[i] = fac[i - 1] * i % mod;
}
ll qpow(ll a, ll b)
{
ll ans = 1;
for (; b; b >>= 1, a = a * a % mod)
if (b & 1)
ans = ans * a % mod;
return ans;
}
ll C(ll n, ll m)
{
if (n < m)
return 0;
return fac[n] * qpow(fac[n - m] * fac[m] % mod, mod - 2) % mod;
}
int main()
{
cal();
int t;
scanf("%d", &t);
while(t--)
{
int n, m;
scanf("%d%d", &n, &m);
ll ans = (C(m + n, n) - C(m - 1, n)+mod)%mod;
printf("%lld\n", ans);
}
}
2. Link with EQ
题意:
有 n n n 个空座位,第一个先随便选一个座位,随后每个人坐的位置都必须满足与其他人的距离的最小值最大,且不能相邻坐。问最后坐的人数的期望。
题解:
设 d p [ i ] [ 1 ] dp[i][1] dp[i][1]表示连续的 i i i 个空位,有一端已经有人坐了,还能坐多少人, d p [ i ] [ 2 ] dp[i][2] dp[i][2]表示连续的 i i i 个空位,两端都已经有人坐了,还能坐多少人。
对于 d p [ i ] [ 1 ] dp[i][1] dp[i][1],那么只能坐另一端,那么转移就是 d p [ i ] [ 1 ] = d p [ i − 1 ] [ 2 ] + 1 dp[i][1]=dp[i-1][2]+1 dp[i][1]=dp[i−1][2]+1。
对于 d p [ i ] [ 2 ] dp[i][2] dp[i][2] ,只能坐中间 , d p [ i ] [ 2 ] = d p [ i / 2 ] [ 2 ] + d p [ i − i / 2 − 1 ] [ 2 ] + 1 dp[i][2]=dp[i/2][2]+dp[i-i/2-1][2]+1 dp[i][2]=dp[i/2][2]+dp[i−i/2−1][2]+1
预处理出所有 d p dp dp后,对于每次询问的 n n n ,可以枚举第一个人坐 i i i ,那么人数就是 d p [ i − 1 ] [ 1 ] + d p [ n − i ] [ 1 ] + 1 dp[i-1][1]+dp[n-i][1]+1 dp[i−1][1]+dp[n−i][1]+1
考虑到 t t t 组,所以可以前缀和预处理,求出总和后除以 n n n即为期望
代码:
#pragma GCC diagnostic error "-std=c++11"
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<map>
#include<stack>
#include<set>
#include<ctime>
#define iss ios::sync_with_stdio(false)
using namespace std;
typedef unsigned long long ull;
typedef long long ll;
typedef pair<int,int> pii;
const int mod=1e9+7;
const int MAXN=1e6+5;
const int inf=0x3f3f3f3f;
ll dp[MAXN][3];
ll pre[MAXN];
ll qpow(ll a,ll b)
{
ll res = 1;
while(b)
{
if(b&1)
res = res * a % mod;
a = a * a % mod;
b >>=1;
}
return res;
}
int main()
{
for (int i = 2; i <= 1e6;i++)
{
dp[i][1] = (dp[i - 1][2] + 1)%mod;
if(i>=3)
{
dp[i][2] = (dp[i / 2][2] + dp[i - i / 2 - 1][2] + 1)%mod;
}
}
for (int i = 1; i <= 1e6;i++)
{
pre[i] = (pre[i - 1] + dp[i][1])%mod;
}
int t;
scanf("%d", &t);
while(t--)
{
int n;
scanf("%d", &n);
ll ans=(pre[n-1]+pre[n-1]+n)%mod*qpow(n,mod-2)%mod;
printf("%lld\n", ans%mod);
}
}