其实直接暴力O(n3)DP+O2O(n^3)DP+O_2O(n3)DP+O2​优化能过…


CODE O(n3)O(n^3)O(n3)

先来个O(n3)O(n^3)O(n3)暴力DP(开了O2O_2O2​)100分代码(极限数据0.5s0.5s0.5s)

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int MAXN = 2005, INF = 0x3f3f3f3f;
int n, m, a[MAXN], f[2][MAXN]; char s[MAXN];
inline void chkmin(int &x, const int &y) { if(y < x) x = y; }
int main () {
freopen("ni.in", "r", stdin);
freopen("ni.out", "w", stdout);
scanf("%d%d", &n, &m);
for(int i = n; i >= 1; --i) {
scanf("%s", s + 1);
for(int j = 1; j <= m; ++j)
if(s[j] == '*' && !a[j]) a[j] = i;
}
int now = 0;
memset(f[0], 0x3f, sizeof f[0]); f[0][a[1]] = 0; //f[i][j]表示i列消得剩下j个的最小步数,暴力刷表
for(int i = 1; i <= m; ++i) {
now ^= 1; memset(f[now], 0x3f, sizeof f[now]);
for(int j = a[i]; ~j; --j) if(f[now^1][j] != INF)
for(int k = j>>1; ~k; --k)
chkmin(f[now][max(a[i+1]+3*k-2*j, 0)], f[now^1][j]+j-k);
}
printf("%d\n", f[now][0]);
}

正解

如上面的暴力DP状态定义一样,只是换成了填表,便于考虑决策单调性.

化简后的方程式为

f[i][j]=f[i−1][2A[i]−2j−3k]+A[i]−j−kf[i][j]=f[i-1][2A[i]-2j-3k]+A[i]-j-kf[i][j]=f[i−1][2A[i]−2j−3k]+A[i]−j−k

单调性证明传送门

这句话:

再观察一下可以发现,在jjj确定的时候,F[i−1][2A[i]−2j−3k]+A[i]−j−kF[i−1][2A[i]−2j−3k]+A[i]−j−kF[i−1][2A[i]−2j−3k]+A[i]−j−k的取值是单峰的。因为F[i−1][2A[i]−2j−3k]F[i−1][2A[i]−2j−3k]F[i−1][2A[i]−2j−3k]随k的减小是不增的,A[i]−j−kA[i]−j−kA[i]−j−k是单增的,因此它是单峰的。

感受一下,大概可以这么想:F[i−1][2A[i]−2j−3k]F[i−1][2A[i]−2j−3k]F[i−1][2A[i]−2j−3k]是不降的,我们感受一下,它的导函数应该不会是什么奇怪的形状,而是单调的…

并且A[i]−j−kA[i]−j−kA[i]−j−k是关于kkk的一次函数,导数为常函数.

所以:

"导函数单调 +++ 一次函数(导函数为常函数) →\to→ 单峰函数(导函数单调)."画出来如下



这样的话f′(x)+g′(x)f'(x)+g'(x)f′(x)+g′(x)如果没有穿过xxx轴,就仍是单调函数.

若穿过了xxx轴,比如先为负,然后为000,然后为正,那么体现在图像上就是单峰了…(好像说了一堆废话)

然后就单调队列优化了.具体代码见下(粘来的代码)

CODE O(n2)O(n^2)O(n2)

#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 2005
using namespace std;
const int inf = 0x3f3f3f3f;
int n,m,a[maxn],f[2][maxn];
char s[maxn];
inline void chkmin(int &a,int b){if(a>b) a=b;}
int main()
{
freopen("ni.in","r",stdin);
freopen("ni.out","w",stdout);
scanf("%d%d",&n,&m);
for(int i=n;i>=1;i--){
scanf("%s",s+1);
for(int j=1;j<=m;j++) if(s[j]=='*'&&!a[j]) a[j]=i;
}
int now=0;memset(f[0],0x3f,sizeof f[0]);
f[0][0]=0;
for(int i=1;i<=m;i++,now=!now){
memset(f[!now],0x3f,sizeof f[!now]);
for(int j=0;j<=a[i];j++){
int id=(a[i]-j)>>1;
for(int k=id;k>=0;k--){
int x=k+(a[i]-j-2*k)*2;
if(x>a[i-1]) x=a[i-1];//emmm...
if(f[!now][j]>=f[now][x]+k+a[i]-j-2*k) id=k,f[!now][j]=f[now][x]+k+a[i]-j-2*k;
else break;
}
}
}
printf("%d",f[now][0]);
}

Newnode's NOI(P?)模拟赛 第二题 dp决策单调优化的更多相关文章

  1. 「模拟赛20190327」 第二题 DP+决策单调性优化

    题目描述 小火车虽然很穷,但是他还是得送礼物给妹子,所以他前往了二次元寻找不需要钱的礼物. 小火车准备玩玩二次元的游戏,游戏当然是在一个二维网格中展开的,网格大小是\(n\times m\)的,某些格 ...

  2. 4.13 省选模拟赛 树 树形dp 卷积 NTT优化dp.

    考试的时候 看到概率 看到期望我就怂 推了一波矩阵树推自闭了 发现 边权点权的什么也不是. 想到了树形dp 维护所有边的断开情况 然后发现数联通块的和再k次方过于困难. 这个时候 应该仔细观察一下 和 ...

  3. NOI.AC NOIP模拟赛 第二场 补记

    NOI.AC NOIP模拟赛 第二场 补记 palindrome 题目大意: 同[CEOI2017]Palindromic Partitions string 同[TC11326]Impossible ...

  4. 6.28 NOI模拟赛 好题 状压dp 随机化

    算是一道比较新颖的题目 尽管好像是两年前的省选模拟赛题目.. 对于20%的分数 可以进行爆搜,对于另外20%的数据 因为k很小所以考虑上状压dp. 观察最后答案是一个连通块 从而可以发现这个连通块必然 ...

  5. [NOI P模拟赛] 传统艺能(子序列自动机、矩阵乘法,线段树)

    (2:00)OID:"完了,蓝屏了!"(代码全消失) 众人欢呼 OID:开机,"原题测试--" (30min later)OID 开始传统艺能: " ...

  6. [NOI.AC]NOI2019省选模拟赛 第二场

    传送门 Solution A. 一共有\(T\)组数据 每次询问你\([l,r]\)中有多少个数能被他的所有数位整除(如果数位中含有\(0\)忽略掉) 数位dp,咕咕咕 B. 题面略 考虑一个个只有两 ...

  7. Newnode's NOI(P?)模拟赛 第三题 (主席树优化建图 + tarjan)

    题目/题解戳这里 这道题题目保证a,b,ca,b,ca,b,c各是一个排列-mdzz考场上想到正解但是没看到是排列,相等的情况想了半天-然后写了暴力60分走人- 由于两两间关系一定,那么就是一个竞赛图 ...

  8. 【XJOI】【NOI考前模拟赛7】

    DP+卡常数+高精度/  计算几何+二分+判区间交/  凸包 首先感谢徐老师的慷慨,让蒟蒻有幸膜拜了学军的神题.祝NOI2015圆满成功 同时膜拜碾压了蒟蒻的众神QAQ 填填填 我的DP比较逗比……( ...

  9. noip模拟赛 水题

    题目描述 LYK出了道水题. 这个水题是这样的:有两副牌,每副牌都有n张. 对于第一副牌的每张牌长和宽分别是xi和yi.对于第二副牌的每张牌长和宽分别是aj和bj.第一副牌的第i张牌能覆盖第二副牌的第 ...

随机推荐

  1. 20191204-使用nginx解决ajax测试调用接口跨域问题

    问题描述 之前要测试一个http的接口,在postman中测试成功,但使用ajax调用却跨域.于是通过使用ngin反向代理的方式来解决ajax调用跨域问题 测试页面的内容 <html> & ...

  2. oracle——学习之路(SQL基础)

    使用create语句创建表 create table 表名 ( 列名 类型 [null     |     not null], 列名 类型 [null     |     not null] ) 在 ...

  3. 接口踩坑:Status (blocked:other)

    1.请求接口时出现 Status (blocked:other) 2.原因分析:安装了Adblock 3.解决办法 1)关掉Adblock2)修改接口名称,不能用 ad 或者 XX ad XX 等名称 ...

  4. Web开发Flask框架学习笔记

    Python 是一种跨平台的[计算机程序设计语言],是一种面向对象的动态类型语言,Python是纯粹的自由软件,源代码和解释器CPython遵循 GPL(GNU General Public Lice ...

  5. 北电之死:谁谋杀了华为的对手?——银湖资本(Silver Lake)董事总经理爱德华·詹德出任CEO,既不了解华为,也不重视中国,直截了当地否决了收购华为

    作者:戴老板:微信公众号:饭统戴老板(ID: worldofboss) 2003年5月,北京SARS疫情紧张,摩托罗拉集团总裁迈克·扎菲罗夫斯基(Mike Zafirovski)却准备不走寻常路,决定 ...

  6. Codeforces Round #406 (Div. 2) A MONSTER

    A. The Monster time limit per test 1 second memory limit per test 256 megabytes input standard input ...

  7. 使用Enablebuffering多次读取Asp Net Core 请求体

    使用Enablebuffering多次读取Asp Net Core 请求体 1 .Net Core 2.X时代 使用EnableRewind倒带 public IActionResult Index( ...

  8. 【原创】大数据基础之Gobblin(2)持久化kafka到hdfs

    gobblin 0.10 想要持久化kafka到hdfs有很多种方式,比如flume.logstash.gobblin,其中flume和logstash是流式的,gobblin是批处理式的,gobbl ...

  9. 04 Go语言之包

    1.为什么有包这个概念? 1)开发中,往往要在不同的文件中调用其他文件的函数 2)Go代码最小粒度单位是”包” 3)go的每一个文件都属于一个包,通过package管理 4)go以包的形式管理文件和项 ...

  10. 转:GitHub团队项目合作流程

    转自:https://www.cnblogs.com/schaepher/p/4933873.html GitHub团队项目合作流程   已在另一篇博客中写出关于以下问题的解决,点此进入: 同步团队项 ...