北京区域赛I题,Uva7676,A Boring Problem,前缀和差分

A Boring Problem

题解

其实这题不难,只要想到了前缀和差分就基本OK了.

我们要求的是第iii项的式子:

F(i)=(a1+a2+...+ai)k+(a2+...+ai)k+...+(ai)kF(i)=(a_1+a_2+...+a_i)^k+(a2+...+a_i)^k+...+(a_i)^kF(i)=(a1+a2+...+ai)k+(a2+...+ai)k+...+(ai)k

Si=a1+a2+...+ai,S0=0S_i = a_1+a_2+...+a_i,S_0=0Si=a1+a2+...+ai,S0=0

F(i)=(Si−S0)k+(Si−S1)k+...+(Si−Si−1)kF(i)=(S_i-S_0)^k+(S_i-S_1)^k+...+(S_i-S_{i-1})^kF(i)=(SiS0)k+(SiS1)k+...+(SiSi1)k

二项式定理展开:

F(i)=∑t=0kCktSit(−S0)k−t+∑t=0kCktSit(−S1)k−t+...+∑t=0kCktSit(−Si−1)k−tF(i)=\sum_{t=0}^kC_k^tS_i^t(-S_0)^{k-t}+\sum_{t=0}^kC_k^tS_i^t(-S_1)^{k-t}+...+\sum_{t=0}^kC_k^tS_i^t(-S_{i-1})^{k-t}F(i)=t=0kCktSit(S0)kt+t=0kCktSit(S1)kt+...+t=0kCktSit(Si1)kt

整理得:

F(i)=∑t=0kCktSit(−1)k−t(S0k−t+S1k−t+...+Si−1k−t)F(i) = \sum_{t=0}^kC_k^tS_i^t(-1)^{k-t}(S_0^{k-t}+S_1^{k-t}+...+S_{i-1}^{k-t})F(i)=t=0kCktSit(1)kt(S0kt+S1kt+...+Si1kt)

再记

SS[i][j]=S0i+S1i+...+Sj−1iSS[i][j] =S_0^i+S_1^i+...+S_{j-1}^iSS[i][j]=S0i+S1i+...+Sj1i

那么

F(i)=∑t=0kCktSit(−1)k−t(SS[k−t][i−1])F(i) = \sum_{t=0}^kC_k^tS_i^t(-1)^{k-t}(SS[k-t][i-1])F(i)=t=0kCktSit(1)kt(SS[kt][i1])

注意到SSSSSS可以O(nk)O(nk)O(nk)预处理出来,SSS可以O(n)O(n)O(n)预处理出来,而F(i)F(i)F(i)就可以O(k)O(k)O(k)得到.

注意!S00=1S_0^0 = 1S00=1

代码

#include <iostream>
#include <algorithm>
#include <cstring>
#define pr(x) std::cout << #x << ':' << x << std::endl
#define rep(i,a,b) for(int i = a;i <= b;++i)typedef long long LL;
const int N = 100010;
const LL P = 1e9+7;
int T,n,k;
char s[N];
long long S[101][N],SS[101][N];
long long C[110][110];
void init() {C[0][0] = 1;for(int i = 1;i <= 100;++i) {C[i][0] = 1;for(int j = 1;j <= i;++j) {C[i][j] = (C[i-1][j-1] + C[i-1][j]) % P;}}
}
int main() {std::ios::sync_with_stdio(false);init();std::cin >> T;while(T--) {std::cin >> n >> k;std::cin >> s;for(int i = 0;i <= n;++i) S[0][i] = 1;for(int i = 1;i <= n;++i) S[1][i] = (s[i-1]-'0') + S[1][i-1] ;for(int i = 2;i <= k;++i) for(int j = 1;j <= n;++j)S[i][j] = S[1][j] * S[i-1][j] % P;SS[0][0] = 1;for(int i = 0;i <= k;++i) {for(int j = 1;j <= n;++j) SS[i][j] = (SS[i][j-1] + S[i][j])% P;}for(int i = 1;i <= n;++i) {long long ans = 0;for(int j = 0;j <= k;++j) {long long res = C[k][j]*S[j][i]%P*SS[k-j][i-1]%P;if((k-j)%2==0) ans = (ans + res) % P;else ans = (ans - res + P) % P;}if(i != 1) std::cout << " ";std::cout << ans;}std::cout << std::endl;}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/320970.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

P4331-[BalticOI2004]Sequence数字序列【左偏树】

正题 题目链接:https://www.luogu.com.cn/problem/P4331 题目大意 给出一个序列aaa&#xff0c;求一个单调上升的序列bbb使得∑i1n∣ai−bi∣\sum_{i1}^n|a_i-b_i|∑i1n​∣ai​−bi​∣最小。 解题思路 巧妙的解法 首先我们让所有的ai−ia_i-iai​−i这样我们求的bbb序列就…

通俗易懂,什么是.NET?什么是.NET Framework?什么是.NET Core?

什么是.NET&#xff1f;什么是.NET Framework?本文将从上往下&#xff0c;循序渐进的介绍一系列相关.NET的概念&#xff0c;先从类型系统开始讲起&#xff0c;我将通过跨语言操作这个例子来逐渐引入一系列.NET的相关概念&#xff0c;这主要包括&#xff1a;CLS、CTS(CLI)、FCL…

【离散化】【DP】命运石之门的选择

命运石之门的选择 题目大意&#xff1a; 有n个盒子&#xff0c;高度为ai&#xff0c;可以数值刷盒子&#xff0c;也可以横着刷&#xff0c;但如果前面没盒子了&#xff0c;就要停下&#xff0c;问刷完这些盒子最少要刷多少次 原题&#xff1a; 题目描述 在某一条不知名世界…

P2153 晨跑,费用流裸题

晨跑 题目连接 https://www.luogu.org/problemnew/show/P2153 题解 求最大不相交路径数,并在路径数最大前提下,求总路程最短. 太裸了. 求不相交路径数:将除1,n1,n1,n两点外的所有点拆分,中间连一条容量为111,费用为000的边.然后所有的原边u→vu \rightarrow vu→v视作从u…

P3812-[模板]线性基

正题 题目链接:https://www.luogu.com.cn/problem/P3812 题目大意 给出nnn个数&#xff0c;求在其中选出若干个数使得它们的异或和最大。 解题思路 序列aaa的线性基bbb满足以下性质 aaa中的任何一个数都可以由bbb中的若干个数异或得到bbb中的任何一个数都不可由bbb中的若干个…

.Net Core小技巧 - Hosted Services + Quartz实现定时任务调度

背景之前一直有朋友问&#xff0c;.Net Core Linux环境有没有类似Windows服务的东西。其实是有的&#xff0c;我了解的方法有两种&#xff1a;#1 创建一个ASP.Net Core的Web项目&#xff08;如Web API&#xff09;&#xff0c;然后通过添加中间件&#xff08;Middleware&#…

初一模拟赛总结(4.7)

成绩&#xff1a; rank是有算其他大佬的 rankrankranknamenamenamescorescorescoreT1T1T1T2T2T2T3T3T3T4T4T4333fyfyfy230230230100100100100100100202020101010444tjhtjhtjh210210210100100100000100100100101010666wjjwjjwjj160160160100100100303030303030000777hkyhkyhky1…

P2053 SCOI2007 修车,费用流好题

修车 题目链接 https://www.luogu.org/problemnew/show/P2053 题解 每个人每次只能修一辆车,且这个人修的最后一辆车所花时间为111倍的修这辆车的时间,修倒数第iii辆车所花的时间是iii倍修这辆车所花的时间. 000号点代表源点,编号为1−M1-M1−M的点代表维修工人,0→[1,M]0 …

P4570-[BJWC2011]元素【线性基,贪心】

正题 题目链接:https://www.luogu.com.cn/problem/P4570 题目大意 给出nnn个物品有aia_iai​和bib_ibi​。要求选出一个bib_ibi​和最大的子集满足其中的aia_iai​不能由其中的其他aia_iai​异或得到 解题思路 我们发现对于一个集合能否加入一个物品其实就是判断aia_iai​是…

在.NET Core微服务中使用HostBuilder和Generic Host

基于控制台的服务中探索一个简单模式&#xff0c;以解决跨领域问题。“通用” Host和HostBuilder是随.NET Core 2.1发布而发布的新功能组件的组件。它们的一个用例是通过提供用于添加横切关注点&#xff08;例如依赖注入&#xff0c;配置和日志记录&#xff09;的模式来简化基于…

【最小环】【Floyed】观光旅游(ssl 1763)

观光旅游 ssl 1763 题目大意&#xff1a; 在一个图中找到最小环 原题&#xff1a; 题目背景 在桑给巴尔岛的Adelton城镇上有一个旅游机构。它们决定在提供许多的其它吸引之外&#xff0c;再向客人们提供旅游本镇的服务。 为了从提供的吸引服务中尽可能地获利&#xff0c;…

P2604 ZJOI2010 网络扩容,费用流裸题

网络扩容 题目链接 https://www.luogu.org/problemnew/show/P2604 题解 对于每条边u→vu \rightarrow vu→v,我们将按照容量CCC,费用000建边,这些算是免费边. 然后我们再对同一对点u,vu,vu,v建立一条容量为INFINFINF,费用为扩容费用WWW的边,这些边算是扩容流量要走的边,有费…

P3292-[SCOI2016]幸运数字【线性基,LCA,倍增】

正题 题目链接:https://www.luogu.com.cn/problem/P3292 题目大意 nnn个点的一棵树&#xff0c;每个点都点权。每次询问一条路径&#xff0c;选择若干个点的异或和最大。 解题思路 路径上的如何进行计算&#xff0c;我们知道我们可以用倍增来计算权值和。我们可以把每个线性基…

EntityFramework Core进行读写分离最佳实践方式,了解一下?

本来打算写ASP.NET Core MVC基础系列内容&#xff0c;博客评论有园友提出讲讲读写分离&#xff0c;这个问题提的好。大多数情况下&#xff0c;对于园友在评论中提出的问题&#xff0c;如果是值得深究或者大多数同行比较关注的问题我都会私下去看看&#xff0c;然后进行对应解答…

P3980 NOI2008志愿者招募

志愿者招募 题目链接 https://www.luogu.org/problemnew/show/P3980 题解 这道题很神奇,这种建图方法很有启发性. 我们平时做的题都点都是是一对一的,而这道题的点的对应关系是一对多(即一个志愿者对应一段连续的区间,也就是多个时间点)的,直接按照传统的网络流建图方法来…

小X的加法难题

小X的加法难题 题目大意&#xff1a; 一个AB的式子&#xff0c;但之间有一些没用的空格&#xff0c;要你求出结果 原题&#xff1a; 解题思路&#xff1a; 字符输入&#xff0c;有数字就存起来&#xff0c;用加号和换行符来分开 代码&#xff1a; #include<cstdio>…

P3857-[TJOI2008]彩灯【线性基】

正题 题目链接:https://www.luogu.com.cn/problem/P3857 题目大意 nnn个彩灯&#xff0c;mmm个开关能使得某些彩灯取反&#xff0c;求有多少种彩灯样式。 解题思路 其实就是mmm个数种若干个数异或起来有多少不同的数。 又是一道考线性基性质的题目&#xff0c;因为线性基中任…

【桶排】小 X 的密码破译

小 X 的密码破译 题目大意&#xff1a; 有一堆式子&#xff0c;去重再排序后&#xff0c;按规则求出结果 原题&#xff1a; 解题思路&#xff1a; 把每一个求出来的数丢到桶里&#xff0c;然后跑一遍就可以了 代码&#xff1a; #include<cstdio> using namespace s…

Uva12325 Zombie's Treasure Chest [二分区间+模拟退火]

Zombie’s Treasure Chest 题目链接 https://cn.vjudge.net/problem/UVA-12325 题意 两种物品无穷多个,第一种物品重量s1s_1s1​,价值v1v_1v1​,第二种物品重量s2s_2s2​,价值v2v_2v2​,背包重nnn,求能装的最大价值之和. 数据全都是2e92e92e9.也就是两种物品的完全背包. 题…

Identity Server 4 预备知识 -- OpenID Connect 简介

我之前的文章简单的介绍了OAuth 2.0 (在这里: 要用Identity Server 4 -- OAuth 2.0 超级简介, 还不是很全.这篇文章我要介绍一下 OpenID Connect.OAuth 2.0 不是身份认证协议OAuth 2.0 不是身份认证(Authentication)协议. 为什么有人会认为OAuth 2.0具有身份认证的功能? 这是因…