HDU - 7048 Link with EQ

题目链接:

https://acm.hdu.edu.cn/showproblem.php?pid=7048

思路分析:

  • 我们可以把状态分为两种。
  • 一种是有两端封闭的连续空区间(长度为 x x x),那么在这个连续空区间内能坐下的人数为: f [ x ] = f [ ⌊ ( x − 1 ) 2 ⌋ ] + f [ ( x − 1 ) − ⌊ ( x − 1 ) 2 ⌋ ] + 1 f[x] = f[\lfloor {\frac{(x - 1)}{2}} \rfloor] + f[(x - 1) - \lfloor {\frac{(x - 1)}{2}} \rfloor ]+ 1 f[x]=f[2(x1)]+f[(x1)2(x1)]+1 即在这个封闭区间选了一个人后,就会得到两个两端封闭的连续空区间,而两个两端封闭的连续空区间的长度分别为 ⌊ ( x − 1 ) 2 ⌋ \lfloor {\frac{(x - 1)}{2}} \rfloor 2(x1) ( x − 1 ) (x - 1) (x1) − - ⌊ ( x − 1 ) 2 ⌋ \lfloor {\frac{(x - 1)}{2}} \rfloor 2(x1)
  • 另一种是一端封闭的连续空区间(长度为 x x x),那么在这个连续空区间内能坐下的人数为: g [ x ] = f [ x − 1 ] + 1 g[x] = f[x - 1] + 1 g[x]=f[x1]+1 因为我们要使离其他人最短距离最大,那么我们只能去取这段连续空区间的另一端的端点,使得满足题目条件-。
  • 注意: f [ 1 ] = f [ 2 ] = 0 f[1] = f[2] = 0 f[1]=f[2]=0 g [ 1 ] = 0 g[1] = 0 g[1]=0
  • 最后我们只需要从桌子上一个人都没有开始即可,我们设能坐下的人为 h [ x ] h[x] h[x],则有: h [ x ] = Σ i = 1 x ( g [ i − 1 ] + g [ x − i ] + 1 ) h[x] = \Sigma_{i = 1}^{x}{(g[i - 1] + g[x - i] + 1)} h[x]=Σi=1x(g[i1]+g[xi]+1)(为空的时候我们随便坐哪都行,所以 i i i是指选的位置,选了之后,容易得到两端长度分别为 i − 1 i - 1 i1 x − i x - i xi
  • 那么期望则为: 1 x × h [ x ] \frac{1}{x}\times h[x] x1×h[x] 1 x \frac{1}{x} x1 是因为第一个选的位置有x种选择),并且h[x]可以化简为: h [ x ] = x + 2 Σ i = 1 x − 1 g [ i ] h[x] = x + 2\Sigma_{i = 1}^{x - 1}g[i] h[x]=x+2Σi=1x1g[i]
  • 所以我们要初始化f[x]、g[x]、sum[g[x]]即可。

代码如下:

#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn = 1e6 + 10;
const ll p = 1e9 + 7;
ll f[maxn], g[maxn], sum[maxn];
void init()
{
        f[1] = 0, f[2] = 0;
        for (int i = 3; i <= maxn; i++)
        {
                f[i] = (f[(i - 1) / 2] + f[(i - 1) - ((i - 1) / 2)] + 1) % p;
        }
        g[0] = 0, g[1] = 0;
        sum[1] = 0;
        for (int i = 2; i <= maxn; i++)
        {
                g[i] = (f[i - 1] + 1) % p;
                sum[i] = (sum[i - 1] + g[i]) % p;
                //前缀和
        }
}
ll qpow(ll a, ll b, ll p)
{
        ll ans = 1;
        while (b)
        {
                if (b & 1)
                {
                        ans = ans * a % p;
                }
                a = a * a % p;
                b >>= 1;
        }
        return ans % p;
}
//求逆元
int main()
{
        int x;
        int t;
        init();
        scanf("%d", &t);
        while (t--)
        {
                scanf("%d", &x);
                ll ans = (1 + 2 * qpow(x, p - 2, p) % p * sum[x - 1] % p) % p;
                printf("%lld\n", ans);
        }
        return 0;
}
  • 3
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值