Atcoder Beginner Contest 200 E. Minflip Summation(概率数论)
题面
一个字符串
T
T
T 是由一个包含 0
、1
、?
的字符串
S
S
S 循环连接
K
K
K 次获得的。
字符串
T
T
T 中的每个 ?
都可以换成 0
或 1
,假设
T
T
T 中一共有
q
q
q 个 ?
,那么把这
q
q
q 个 ?
全部换成 0
或 1
后能得到
2
q
2^q
2q 种不同的
T
′
T'
T′ ,现在求
∑
F
(
T
′
)
\sum F(T')
∑F(T′)
即把这
2
q
2^q
2q 个
F
(
T
′
)
F(T')
F(T′) 算出来并求和。
其中的
F
(
T
′
)
F(T')
F(T′) 定义为:用如下操作把
T
′
T'
T′ 的所有字符变得相同的最小操作数:
选择一个区间
[
l
,
r
]
[l,r]
[l,r] ,把
T
l
′
T'_l
Tl′ 到
T
r
′
T'_r
Tr′ 之间(包括端点)的所有字符都异或一.
- 1
≤
∣
S
∣
≤
1
0
5
1\leq |S|\leq 10^5
1≤∣S∣≤105
- 1
≤
K
≤
1
0
9
1\leq K\leq 10^9
1≤K≤109
题解
Subtask 1:“无惑”
首先,我们想想对于没有问号的一个序列
S
S
S 怎么做。
这应该是经典的差分思想了,把整个数组差分,每次操作相当于令两个不同位置分别异或 1 ,即每次操作能消掉两个 1。
我们再定量地分析,不妨设整个数组最终全变成序列的第一个数字(这样一定不会更优了),那么除了开头以外的每一段连续的 1 或 0 都会贡献一个 1(在差分数组上),然后如果最后一个数字与开头不同,那么一定还要在最后填一个 1,为了使数组全部相同。此时,容易发现 1 的总数一定是偶数,因为我们的值从
S
1
S_1
S1 又回到了
S
1
S_1
S1。所以答案就是 1 的个数除以二,即:
∑
i
=
1
∣
S
∣
(
S
i
=
=
S
i
%
∣
S
∣
+
1
?
0
:
1
)
2
\cfrac{\sum_{i=1}^{|S|}(S_i==S_{i\%|S|+1}?\;\;0:1)}{2}
2∑i=1∣S∣(Si==Si%∣S∣+1?0:1)
那么如果这个序列连续拼接
K
K
K 个呢?由于最后一个字符刚好和下一个串的第一个字符相邻,容易发现答案是
K
∗
∑
i
=
1
∣
S
∣
(
S
i
=
=
S
i
%
∣
S
∣
+
1
?
0
:
1
)
2
\cfrac{K*\sum_{i=1}^{|S|}(S_i==S_{i\%|S|+1}?\;\;0:1)}{2}
2K∗∑i=1∣S∣(Si==Si%∣S∣+1?0:1)
即之前的答案直接乘个
K
K
K 就是了。
Subtask 2:“释疑”
如果有问号呢?
一段连续的问号转变后可能千变万化,毫无内在联系,我们换个角度想。
由于每个问号互不相干,它们变成 1 或 0 的概率也互不相干,即:设
S
i
=
1
(
/
0
)
S_i=1(/0)
Si=1(/0) 为事件
A
A
A,
S
j
=
1
(
/
0
)
(
j
≠
i
)
S_j=1(/0)~(j\not=i)
Sj=1(/0) (j=i) 为事件
B
B
B ,则
P
(
A
∣
B
)
=
P
(
A
)
,
P
(
B
∣
A
)
=
P
(
B
)
P(A|B)=P(A),P(B|A)=P(B)
P(A∣B)=P(A),P(B∣A)=P(B)。
那么我们可以分开考虑每相邻两个数不同的概率,然后算出期望不同总数,(设
S
S
S 中问号个数为
q
q
q)答案即为
K
∗
E
(
∑
i
=
1
∣
S
∣
(
S
i
=
=
S
i
%
∣
S
∣
+
1
?
0
:
1
)
)
2
∗
2
q
K
\cfrac{K*E\Big(\sum_{i=1}^{|S|}(S_i==S_{i\%|S|+1}?\;\;0:1)\Big)}{2}*2^{qK}
2K∗E(∑i=1∣S∣(Si==Si%∣S∣+1?0:1))∗2qK
现在分情况考虑相邻的数,定义一个函数
f
(
i
,
j
)
=
f(i,j)=
f(i,j)=
- S
i
,
S
j
S_i,S_j
Si,Sj 其中一个是问号:
1
+
0
2
=
1
2
\frac{1+0}{2}=\frac{1}{2}
21+0=21 (一半的概率不同)
- S
i
,
S
j
S_i,S_j
Si,Sj 都是问号:
1
+
1
+
0
+
0
4
=
1
2
\frac{1+1+0+0}{4}=\frac{1}{2}
41+1+0+0=21 (四种情况,两种情况不同)
- S
i
,
S
j
S_i,S_j
Si,Sj 都不是问号:
S
i
=
=
S
j
?
0
:
1
S_i==S_j?\;\;0:1
Si==Sj?0:1
把分开的期望加起来,答案为:
K
∗
∑
i
=
1
∣
S
∣
f
(
i
,
i
%
∣
S
∣
+
1
)
2
∗
2
q
K
\cfrac{K*\sum_{i=1}^{|S|}f(i,i\%|S|+1)}{2}*2^{qK}
2K∗∑i=1∣S∣f(i,i%∣S∣+1)∗2qK
CODE
#include<map>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define MAXN 100005
#define ENDL putchar('\n')
#define LL long long
#define DB double
#define lowbit(x) ((-x) & (x))
#define eps 1e-9
LL read() {
LL f = 1,x = 0;char s = getchar();
while(s < '0' || s > '9') {if(s=='-')f = -f;s = getchar();}
while(s >= '0' && s <= '9') {x=x*10+(s-'0');s = getchar();}
return f * x;
}
const int MOD = 1000000007,inv2 = 1000000008/2;
int n,m,i,j,s,o,k;
char ss[MAXN];
int t[MAXN];
int qkpow(int a,LL b) {
int res = 1;
while(b > 0) {
if(b & 1) res = res *1ll* a % MOD;
a = a *1ll* a % MOD; b >>= 1;
}return res;
}
int main() {
scanf("%s",ss + 1);
n = strlen(ss + 1);
m = read();
if(n*1ll*m == 1) {printf("0\n");return 0;}
int cn = 0,cnt = 0;
LL cn2 = 0;
for(int i = 1;i <= n;i ++) if(ss[i] == '?') cnt ++;
for(int i = 1;i <= n;i ++) {
int ct = 0,j = i+1; if(j > n) j = 1;
if(ss[i] == '?') ct ++;
if(ss[j] == '?') ct ++;
if(ct) (cn += inv2) %= MOD;
else {
if(ss[i] != ss[j]) (cn += 1) %= MOD,cn2 ++;
}
}
cn = cn *1ll* m % MOD;
int ans = cn *1ll* qkpow(qkpow(qkpow(2,cnt),m),1) % MOD;
ans = ans *1ll* inv2 % MOD;
printf("%d\n",ans);
return 0;
}
Atcoder Beginner Contest 200 E. Minflip Summation(概率数论)的更多相关文章
- KYOCERA Programming Contest 2021(AtCoder Beginner Contest 200) 题解
KYOCERA Programming Contest 2021(AtCoder Beginner Contest 200) 题解 哦淦我已经菜到被ABC吊打了. A - Century 首先把当前年 ...
- AtCoder Beginner Contest 200
前言:果然自己连\(ABC\)都打不好了吗. 没看清题目,卡了巨久,排名一直跌,笔记本键盘坏了,心态崩了. 冷静. ------------------------------------------ ...
- AtCoder Beginner Contest 084 D - 2017-like Number【数论/素数/前缀和】
D - 2017-like Number Time limit : 2sec / Memory limit : 256MB Score : 400 points Problem Statement W ...
- AtCoder Beginner Contest 100 2018/06/16
A - Happy Birthday! Time limit : 2sec / Memory limit : 1000MB Score: 100 points Problem Statement E8 ...
- AtCoder Beginner Contest 053 ABCD题
A - ABC/ARC Time limit : 2sec / Memory limit : 256MB Score : 100 points Problem Statement Smeke has ...
- AtCoder Beginner Contest 076
A - Rating Goal Time limit : 2sec / Memory limit : 256MB Score : 100 points Problem Statement Takaha ...
- AtCoder Beginner Contest 068 ABCD题
A - ABCxxx Time limit : 2sec / Memory limit : 256MB Score : 100 points Problem Statement This contes ...
- AtCoder Beginner Contest 052
没看到Beginner,然后就做啊做,发现A,B太简单了...然后想想做完算了..没想到C卡了一下,然后还是做出来了.D的话瞎想了一下,然后感觉也没问题.假装all kill.2333 AtCoder ...
- AtCoder Beginner Contest 136
AtCoder Beginner Contest 136 题目链接 A - +-x 直接取\(max\)即可. Code #include <bits/stdc++.h> using na ...
随机推荐
- ngx_http_fastcgi_module 的那些事
是什么? 顾名思义,是Nginx用来处理FastCGI的模块.FastCGI是什么?这个以后再讲,可以说的是现在LNMP架构里面,PHP一般是以PHP-CGI的形式在运行,它就是一种FastCGI,我 ...
- WAVE音频格式及及转换代码
音频信号的读写.播放及录音 python已经支持WAV格式的书写,而实时的声音输入输出需要安装pyAudio(http://people.csail.mit.edu/hubert/pyaudio).最 ...
- Ubuntu-apt-get出错
ubuntu apt-get的时候出现: E: Could not open lock file /var/lib/dpkg/lock-frontend - open (13: Permission ...
- Vue3.0系列——「vue3.0性能是如何变快的?」
前言 先学习vue2.x,很多2.x内容依然保留: 先学习TypeScript,vue3.0是用TS重写的,想知其然知其所以然必须学习TS. 为什么学习vue3.0? 性能比vue2.x快1.2-2倍 ...
- cut命令、case与select语句
cut命令 常用参数: -c character 字符 -d delimiter 分隔符 -f field 域(列) --output-delimiter 输出分隔符 例: # echo 12345 ...
- 深入解析kubernetes controller-runtime
Overview controller-runtime 是 Kubernetes 社区提供可供快速搭建一套 实现了controller 功能的工具,无需自行实现Controller的功能了:在 Kub ...
- 「BUAA OO Unit 4 HW16」第四单元总结与课程回顾
「BUAA OO Unit 4 HW16」第四单元总结与课程回顾 目录 「BUAA OO Unit 4 HW16」第四单元总结与课程回顾 Part 0 第四单元作业架构设计 架构设计概要 AppRun ...
- Python简单实现自动评论、自动点赞、自动关注脚本
一些哔哔: 今天的这个脚本,是一个别人发的外包,交互界面的代码就不在这里说了,但是可以分享下自动评论.自动点赞.自动关注.采集评论和视频的数据是如何实现的 开发环境 python 3.8 运行代码py ...
- Python爬虫常用:谷歌浏览器驱动——Chromedriver 插件安装教程
我们在做爬虫的时候经常要使用谷歌浏览器驱动,今天分享下这个Chromedriver 插件的安装方法. 第一步:打开谷歌浏览器打开设置面板 嫌枯燥的小伙伴可以点击此处找管理员小姐姐领取免费资料 第二步: ...
- 关于cpu体系架构的一些有趣的故事分享
从排查一次匪夷所思的coredump,引出各种体系架构的差异. 本文中的所有内容来自学习DCC888的学习笔记或者自己理解的整理,如需转载请注明出处.周荣华@燧原科技 1 背景 从全世界有记载的第一台 ...