[USACO19JAN,Platinum]Train Tracking 2

虽然是简单的dp ,but真的太难想到了,而且我的码力。。。

Train Tracking 2

【题目描述】
每天特快列车都会经过农场。列车有N节车厢(1≤N≤105),每节车厢上有一个1到109之间的正整数编号;不同的车厢可能会有相同的编号。 平时,Bessie会观察驶过的列车,记录车厢的编号。但是今天雾实在太浓了,Bessie一个编号也看不见!幸运的是,她从城市里某个可靠的信息源获知了列车编号序列的所有滑动窗口中的最小值。具体地说,她得到了一个正整数K,以及N−K+1个正整数c1,…,cN+1−K,其中ci是车厢i,i+1,…,i+K−1之中编号的最小值。

帮助Bessie求出满足所有滑动窗口最小值的对每节车厢进行编号的方法数量。由于这个数字可能非常大,只要你求出这个数字对10^9+7取余的结果Bessie就满意了。

Bessie的消息是完全可靠的;也就是说,保证存在至少一种符合要求的编号方式。

【输入输出格式】
输入的第一行包含两个空格分隔的整数N和K。余下的行包含所有的滑动窗口最小值c1,…,cN+1−K,每行一个数。

输出一个整数:对每节车厢给予一个不超过109的正整数编号的方法数量对10^9+7取余的结果,满足车厢i,i+1,…,i+K−1之中编号的最小值等于ci,对于1≤i≤N−K+1均成立。

【输入输出样例】
输入
4 2
999999998
999999999
999999998
输出
3
【分析】
1:999999998 999999999 999999999 999999998
2:999999998 999999999 1000000000 999999998
3:999999998 1000000000 999999999 999999998

【解题思路】
首先你嘚想到相邻的i,ji,jij如果它们的最小值不一样,那么一定可以固定一个点

eg:iiii+k−1i+k-1i+k1的最小值4,j4,j4,jj+k−1j+k-1j+k1的最小值333iii小于jjj,即j=i+1j=i+1j=i+1,这个时候就固定了j+k−1j+k-1j+k1得是333

从中我们的启发是什么

把c值相同的连续区间单独拎出来

ci=ci+1=ci+2=...=cj=vci=ci+1=ci+2=...=cj=vci=ci+1=ci+2=...=cj=v

这一段代表的区间总长为j−i+kj−i+kji+k,证明一下:iiijjj区间长为j−i+1j-i+1ji+1没问题

那么jjj包含的区间长是kkk,加在一起再减去一个重复的jjj就是j−i+kj-i+kji+k,继续

所有的数都大于等于vvv,同时每kkk个中就有至少一个vvv

有一个直接的dp:设dp[i]dp[i]dp[i]表示最后一个vvviii位置时的合法方案数,ppp1e9−v1e9−v1e9v

那么有转移dp[i]=∑(j=i−k,i−1)pi−j−1∗f[j]dp[i]=∑( j=i−k, i−1)p^i−j−1*f[j]dp[i]=j=iki1pij1f[j]

这样是O(nk)O(nk)O(nk)的,会TLE

稍微改变一下错位相消就变成了:

dp[i]=(p+1)*f[i−1]−p^k∗f[i−k−1]

在这里要解释一下首先如果有最小值是xxx,那么方案数就是(1e9−x+1)j−i+1(1e9-x+1)^{j-i+1}(1e9x+1)ji+1减去不合法的方案数

[i,j][i,j][ij]区间不合法的方案数其实是所有取值全都大于了vvv,没有一个是vvv,其实也就是(1e9−x+1)j−i+1(1e9-x+1)^{j-i+1}(1e9x+1)ji+1

问题是为什么要乘以dp[i−k−1]dp[i-k-1]dp[ik1]呢?

你想想在iii不合法的方案数,在i−1i-1i1时是属于合法的方案数,所以影响iii的是i−k−1i-k-1ik1的所有方案
因为从i−k−1i-k-1ik1开始往前不管选什么都与iii无关,所以后面是不合法的,前面不管合不合法都算不合法,所以要相乘

接着就是考虑如何处理两个vvv不一样的交界处
说白了就是两个vvv谁更大,相交的部分就交给谁处理在这里插入图片描述
因为这就好比不等式解集的取法,如果这个香蕉区间既要满足最小值为5又要满足最小值为6那坑定是选满足最小值为6啊!

如果处理到iii时,iii有一部分区间被前面处理过,那么就要减去kkk

因为我们是循环iii是每次只加111所以当最小值不一样的时候,两个区间的香蕉部分就是k−1k-1k1,又因为这k−1k-1k1被前面处理过,也就意味着iii这个点为了满足最小值也是固定的,那么真正波动的有选择的长度就减去了kkk

同理如果处理到iii时,发现i+1i+1i+1的最小值更大,也就意味着有k−1k-1k1区间是后面处理,i也固定下来了。长度便又减去了kkk

最后就是注意头i=1i=1i=1和尾i=n−k+1i=n-k+1i=nk+1时的特判,如果i=1i=1i=1它就不能进行把k−1k-1k1区间交给前面的人处理,因为它前面没人;如果i=n−k+1i=n-k+1i=nk+1它就不能进行把k−1k-1k1区间交给后面的人处理,因为它后面也没人了。
在这里插入图片描述
在这里,要重点解释一下下面代码的dp初值设定,早上吃饭时突然灵机乍现,明白这样写的原因

为什么?dp[0]和dp[1]初值是1 为什么?window里i循环要从2开始
1)dp[0]的情况,就是这一段区间刚好就是1~k,
那么当我们算到dp[k]的时候,它是包含了从1~k所有都取大于min的值的不合法的方案数
而我们要减掉这些不合法的方案数,也就是每一个数的取值都是min+1~1e9,有k个
所以不合法方案数就是(1e9-(min-1)+1)^k也就是代码中的tp,
那么这个系数是由dp[0]来提供,因而要赋值成1

2)dp[1]的情况,就是这一段区间刚好是1~k+1
而我们能进入这个函数一定要满足这个区间所有数的最小值是一样的val
所以当我们算到dp[k+1]的时候
不合法的方案数=2~k全都不是最小值方案数*dp[1]的方案数
那dp[1]不应该是tp吗?
no~~因为我们的通向dp式dp[i-k-1]是乘以(1e9-min)的K次方,而2到k全都不是最小值方案数是(1e9-min)的K-1次方,dp[1]的所有方案数是(1e9-min)相乘刚好就是k次方
那么这时候的系数又是由dp[1]提供的,就只能赋值成1

同样这也是为什么我们的i是从2开始循环,为了考虑这两种情况就要给0,1腾位置,

【代码实现】

#include <cstdio>
const int mod = 1e9 + 7;
const int p = 1e9;
#define MAXN 100005
#define LL long long
int n, k;
LL c[MAXN];
LL result = 1;
LL dp[MAXN];LL qkpow ( LL x, LL y ) {LL res = 1;while ( y ) {if ( y % 2 ) res = ( res * x ) % mod;x = ( x * x ) % mod;y >>= 1;}return res % mod;
}LL window ( LL val, LL len ) {LL tmp = p - val;LL tp = qkpow ( tmp, k );dp[0] = dp[1] = 1;for ( int i = 2;i <= len + 1;i ++ ) {dp[i] = ( ( tmp + 1 ) * dp[i - 1] ) % mod;if ( i - k - 1 >= 0 ) dp[i] = ( ( dp[i] - ( ( tp * dp[i - k - 1] ) % mod ) + mod ) % mod ) % mod;}return dp[len + 1] % mod;
}int main() {scanf ( "%d %d", &n, &k );for ( int i = 1;i <= n - k + 1;i ++ )scanf ( "%lld", &c[i] );for ( int i = 1, j;i <= n - k + 1;i = j + 1 ) {LL len;j = i;while ( c[i] == c[j + 1] )j ++;len = j - i + k;if ( i != 1 && c[i - 1] > c[i] )len -= k;if ( j != n - k + 1 && c[j + 1] > c[i] )len -= k;if ( len > 0 )result = ( result * window ( c[i], len ) ) % mod;}printf ( "%lld", result % mod );return 0;
} 

在这里插入图片描述

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

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

相关文章

.netcore 堆栈调用方法小记

1|0背景上午临近午饭时&#xff0c;公司同事反馈验证码被攻击灌水。我们匆忙查询验证码明细&#xff0c;对已频繁出现的IP插入黑名单&#xff0c;但IP仍然隔断时间频繁变动&#xff0c;不得已之下只能先封禁对应公司id的验证码发送功能。年初时候&#xff0c;专门对SSO站点的发…

洛谷P2480:古代猪文(中国剩余定理)(欧拉定理)

传送门 文章目录题目描述解析总结代码题目描述 解析 简单来说&#xff0c;就是求&#xff1a; g∑C(d,n)(d是n的约数&#xff09;mod 999911659 可以先特判一下&#xff0c;999911659|g时&#xff0c;答案为0 否则&#xff0c;可以通过欧拉定理转化为&#xff1a; g∑C(d,n)(d…

杯子 + Kronican

杯子 Kronican 【题目描述】 重庆八中在80周年校庆的时候获捐n个杯子&#xff0c; 每个杯子有两个属性&#xff1a;一个是已装水量 ai&#xff0c;一个是可装水量 bi&#xff08;ai < bi&#xff09;。 从一个杯子向另一个杯子倒 x 单位体积的水需要花费的时间是 x 秒。 现…

微软开源的Trill是什么?

以下是一篇15年的文章的译文&#xff1a;https://dwainegilmer.wordpress.com/2015/01/28/microsoft-trill-for-streaming-analytics-from-microsoft-research/当今许多大数据应用程序套件的重点是数据存储。它们是围绕狭窄范围的数据集设想和设计的&#xff0c;通常是为了组织…

剪纸游戏(博弈论)(SG函数)

文章目录题目描述解析题目描述 解析 本题的关键就是SG函数的定义 尝试了一些自己直观上可能对但题解没有使用的约定方法&#xff08;当然最后证明都是错的 。。。&#xff09;&#xff0c;对SG的理解更深刻了一些 SG0的含义是无法再移动换句话说也就是再移动也对败局于事无补&…

[COCI2017-2018#5] Karte

[COCI2017-2018#5] Karte&#xff0c;简短的代码想到了就AC 这道题是SPJ放心搞 但是我的脑子里面的东西&#xff0c;不用我多说&#xff0c;你们就知道是水和面粉和成的 看招 题 【题目描述】 你有一副共有N张牌的牌&#xff0c;在第i张牌上会有一个数字ai表示在这张牌下面至…

YBTOJ:红与蓝(博弈论)

文章目录题目描述解析代码题目描述 解析 首先&#xff0c;这道题的情境对二人来说是不对称的&#xff0c;所以不太好使用SG函数来求解 但直观上也好考虑 利用树的递归性质可以求出每个节点的颜色是否确定 并确定根的颜色是否确定 如果确定是红就随便涂 确定是蓝就-1 关键在于不…

eShopOnContainers 看微服务④:Catalog Service

服务简介Catalog service&#xff08;目录服务&#xff09;维护着所有产品信息&#xff0c;包括库存、价格。所以该微服务的核心业务为&#xff1a;产品信息的维护库存的更新价格的维护架构模式先看代码结构&#xff08;下图&#xff09;。主要依赖&#xff1a;1、HealthCheck …

【CF 1195】Basketball Exercise/Submarine in the Rybinsk Sea (hard edition)/OpenStreetMap+二维单调队列滑动窗口模板

寡人认为C&#xff0c;E都是比较板的题 butD2也太ex了&#xff0c;大大是被那个mod精给弄疯了&#xff0c;我mod了那么多次还是炸了longlong orz 文章目录二维单调队列模板C&#xff1a;Basketball Exercise题目大意题解代码实现D2&#xff1a;Submarine in the Rybinsk Sea (…

YBTOJ:方程的解(组合数学)(插板法)

文章目录题目描述解析代码题目描述 解析 第一感觉&#xff1a;啥都没感觉出来。。。 直接拿动态规划高精做的 但是只能拿40 重新分析一下这道题&#xff1a; g&#xff08;x&#xff09;首先可以拿快速幂很容易的求出来 问题就转化为了**把g(x)个东西分成k份的方案数 其实答案…

开源库Magicodes.Storage正式发布

说明Magicodes.Storage&#xff0c;是心莱科技团队提供的统一存储库&#xff0c;相关库均使用.NET标准库&#xff08;netstandard2.0&#xff09;编写&#xff0c;支持.NET Framework以及.NET Core。我们希望&#xff0c;使用了Magicodes.Storage之后&#xff0c;开发者可以很快…

【CF 1188 A1,B,C】Add on a Tree // Count Pairs // Array Beauty

传送门 这些天风也温柔&#xff0c;题也温柔 开车啦&#xff01; 文章目录A1&#xff1a;Add on a Tree题意翻译题解证明代码实现B&#xff1a;Count Pairs题意翻译题解代码实现C&#xff1a;Array Beauty题目描述题解代码实现A1&#xff1a;Add on a Tree 题意翻译 给定一棵…

eShopOnContainers 知多少[5]:EventBus With RabbitMQ

1. 引言事件总线这个概念对你来说可能很陌生&#xff0c;但提到观察者&#xff08;发布-订阅&#xff09;模式&#xff0c;你也许就很熟悉。事件总线是对发布-订阅模式的一种实现。它是一种集中式事件处理机制&#xff0c;允许不同的组件之间进行彼此通信而又不需要相互依赖&am…

YBTOJ:数列方案(组合数学)

文章目录题目描述解析代码题目描述 解析 如果它不取等&#xff0c;那就和方程的解这道题一样了&#xff0c;但有了等号就很头疼 如何把等号去掉呢&#xff1f; 定义BiAiiB~i~A~i~iB i A i i那么我们就可以得到&#xff1a;0<B1<B2<...<Bm<mn0<B~1~<B~2~&…

【COCI 2018/2019 Round #2】Kocka

这道题也是一个ex的模拟题 不过他比Zamjena可爱 作为一个帅气的小哥哥&#xff0c;让我们一起&#xff0c; 开启你的模拟ex大门&#xff0c;C从入门到放弃&#xff01; 题目 题目描述 我又来了&#xff01;我又来了&#xff01; 在清晨来到儿童游乐园的时候&#xff0c;出题…

Matrix Equation

题意&#xff1a; 题目给出两个矩阵X,Y,现在有两种操作 Z X Y D X⊙Y 问是否存在一个矩阵C&#xff0c;使得ACB⊙C式子成立&#xff0c;问矩阵C能有多少个 题解&#xff1a; 这个式子在模2意义下的加法就等于异或 也就相当于 那现在有 将BC移到左边 然后将Ci,j的系数进…

eShopOnContainers 知多少[6]:持久化事件日志

1. 引言事件总线解决了微服务间如何基于集成事件进行异步通信的问题。然而只有事件总线正常运行&#xff0c;微服务之间基于事件的通信才得以运转。 而现实情况是&#xff0c;总有这样或那样的问题&#xff0c;导致事件总线不稳定或不可用&#xff0c;比如&#xff1a;网络中断…

多体问题

代码&#xff1a; function SunEarthMoon % M函数文件load planets; % 将planets.mat中的变量mass、position、velocity加载过来[sun, earth, moon] deal(18, 3, 25); % sun、earth、moon分别是18、3、25行 list [sun, earth, moon]; % 1行3列矩阵 G 6.67e-11; % gr…

【CF1179 A,B,C】Valeriy and Deque / Tolik and His Uncle / Serge and Dining Room

还好题很温柔&#xff0c;温柔得我差点没做完 文章目录A&#xff1a;Valeriy and Deque题意题解代码实现B&#xff1a;Tolik and His Uncle题目题解代码实现C&#xff1a;Serge and Dining Room题目题解代码实现A&#xff1a;Valeriy and Deque 题意 给定一个双端队列&#…

YBTOJ:比赛得分(期望)

文章目录题目描述解析代码题目描述 解析 不太难的题 显然本题在AB队员大小关系相反时其对答案的贡献互为相反数。 所以想到把B队队员sort一下后就可以二分找到大小关系相反的分界点 然后维护和与平方和两个前缀数组搞一搞即可O1求出贡献 总复杂度&#xff1a;nlognnlognnlogn …