2021杭电多校(7) 补题

本文介绍了两道编程竞赛题目,第一题涉及组合数学与隔板法的运用,解析了如何通过组合公式计算特定条件下拿球的方案数。第二题探讨了序列填充问题的解决策略,利用动态规划计算满足特定条件的座位填充分布的期望值。两题均展示了数学思维在编程中的重要性。
摘要由CSDN通过智能技术生成

传送门

1. Link with Balls

题意:

2 n 2n 2n个桶,对于第 2 x − 1 2x-1 2x1个桶,可以从里面得到 k x kx kx个球,对于第 2 x 2x 2x个桶,最多可以拿 x x x个球,求最后恰好拿到 m m m个球的方案数。

题解:

对于第 2 x − 1 2x-1 2x1个桶,可以拿 k x kx kx,对于第 2 x − 2 2x-2 2x2个桶,可以拿 [ 0 , x − 1 ] [0,x-1] [0,x1] ,那么这两个结合在一起,就可以拿任意数量的球,

那么我们可以让相邻的桶两两结合: ( 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)...(2n2,2n1),(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=0nCmi+n1n1=Cm+nnCm1n

代码:

#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[i1][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[ii/21][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[i1][1]+dp[ni][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);
    }
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值