2020百度之星初赛一 Rotate(期望)

Problem Description
我们有一个圈,从内到外一共被分成了 nn 个环,中间是空的。

我们把从外到内第 ii 层环平分成 a[i]a[i] 份,其中 a[i]a[i] 是偶数,我们把这 a[i]a[i] 份黑白染色,第奇数个染成黑色,第偶数个染成白色。

现在我们旋转每一层,每一层都会等概率随机到一个中止位置。

问黑色的联通块数目的期望。两块黑色的区域有交点即算联通。层之间的旋转是相互独立的。

Input
第一行一个正整数 test(1 \le test \le 10)test(1≤test≤10) 表示数据组数。

对于每组数据,第一行一个正整数 n(1 \le n \le 10)n(1≤n≤10)。

接下来一行 nn 个正整数 a[i](2 \le a[i] \le 1000)ai,a[i]a[i] 为偶数,另外保证 aa 序列不降。

Output
对于每组数据,一行一个数表示答案,由于答案 A/BA/B 中的 ABAB 可能很大,请输出 A/B \mod 10^9+7A/Bmod10
9
+7,假设 A/BA/B 为最简分数,A/B \mod 10^9+7 = A * B^{-1} \mod 10^9+7A/Bmod10
9
+7=A∗B
−1
mod10
9
+7,B^{-1}B
−1
为满足 B^{-1}*B \mod 10^9+7 = 1B
−1
∗Bmod10
9
+7=1 的整数。

Sample Input
3
2
2 2
2
2 6
4
10 230 666 1000
Sample Output
1
2
500000256

样例解释
对于第一组样例,第一个环和第二个环各有一半是黑的,他们的黑色部分必然有交点,所以期望的联通块数目为 1。
对于第二组样例,旋转第一个环,只在几个特定角度上第一个环的黑色部分和第二个环的三个黑色部分都有交点,大部分情况下都只和两个黑色部分有交点,因为算期望,所以几个特定角度对答案没有影响,答案为 2。

思路:
感觉很有意思的一道期望题,唯一需要考虑只是每个环的划分数。

只需要 a [ i ] a[i] a[i]是单调的,递增递减问题不大。
不妨把本题看做单调递增,那么不存在内环向相连外环,一个黑块连接大于一个黑块的情况。

把黑块看作点,那么最后构成了一个森林(a单调的结论得来),结果就是总的点数减去边数。

点数为 ∑ a [ i ] 2 ∑\frac{a[i]}{2} 2a[i]

假设外环有 a [ i ] a[i] a[i]个,内环有 a [ i − 1 ] a[i-1] a[i1]个,
则里面外环的每个黑块,连接了内环 ( a [ i − 1 ] a [ i ] + 1 ) 2 \frac{(\frac{a[i-1]}{a[i]}+1)}{2} 2(a[i]a[i1]+1)个黑块,比如2和4,那么外环的黑块连接了内环3个块(可能出现连接为2个的情况,但因为只有一个这种情况,所以概率为0),每个块有一半几率是黑块,相当于3/2个黑块。

而外环有 a [ i ] 2 \frac{a[i]}{2} 2a[i]个黑块,所以相邻两个环的边数为 a [ i ] + a [ i − 1 ] 4 \frac{a[i]+a[i-1]}{4} 4a[i]+a[i1]

用点数减去边数,得到结果 a [ 1 ] + a [ n ] 4 \frac{a[1]+a[n]}{4} 4a[1]+a[n]

upd:
在这里插入图片描述
群里大佬的解释,是一个意思,但是应该更加严谨一些。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <string>
#include <cmath>
#include <iostream>
#include <stack>
#include <vector>
#include <queue>
using namespace std;

typedef long long ll;

const int maxn = 1005;
const int mod = 1e9 + 7;
ll a[maxn];

ll qpow(ll x,ll n) {
    ll res = 1;
    while(n) {
        if(n & 1) res = res * x % mod;
        x = x * x % mod;
        n >>= 1;
    }
    return res;
}

int main() {
    int T;scanf("%d",&T);
    while(T--) {
        int n;scanf("%d",&n);
        for(int i = 1;i <= n;i++) {
            scanf("%lld",&a[i]);
        }
        printf("%lld\n",(a[1] + a[n]) * qpow(4,mod - 2) % mod);
    }
    return 0;
}
  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 4
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值