【重庆市NOIP模拟赛】数据

最新推荐文章于 2018-11-01 18:05:59 发布
最新推荐文章于 2018-11-01 18:05:59 发布
阅读量1k 收藏 1
点赞数 1
分类专栏： DP 文章标签： DP优化
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/C20181220_xiang_m_y/article/details/76033724
版权

数据
时间限制: 1 Sec 内存限制: 128 MB
提交: 58 解决: 31
[提交][状态][我的提交]
题目描述
Mr_H 出了一道信息学竞赛题，就是给 n 个数排序。输入格式是这样的：
试题有若干组数据。每组数据的第一个是一个整数 n，表示总共有 n 个数待排序；接下来 n 个整数，分别表示这 n 个待排序的数。
例如：3 4 2 –1 4 1 2 3 4，就表示有两组数据。第一组有 3 个数(4,2,-1)，第二组有 4个数(1,2,3,4)。可是现在 Mr_H 做的输入数据出了一些问题。例如：2 1 9 3 2 按理说第一组数据有 2 个数(1,9)，第二组数据有 3 个数，可是“3”后面并没有出现三个数，只出现了一个数“2”而已！
现在 Mr_H 需要对数据进行修改，改动中“一步”的含义是对文件中的某一个数+1 或-1，写个程序，计算最少需要多少步才能将数据改得合法。

输入
第一行一个整数 m，表示 Mr_H 做的输入数据包含的整数个数。第二行包含 m 个整数 a[i]，每个整数的绝对值不超过 10000。

输出
一个整数，表示把数据修改为合法的情况下，最少需要多少步。

样例输入
Copy (如果复制到控制台无换行，可以先粘贴到文本编辑器，再复制)

4
1 9 3 2
样例输出
2
提示
对于 20%的数据，m<=10, |a[i]|<=5;

对于 60%的数据，m<=5000, |a[i]|<=10000

对于 100%的数据，m<=100000, |a[i]|<=10000

法一：图论。

  • 每个点i可以零损耗地到达第i+a[i]+1个点，所以连一条权值为零的单向边。
  • 同时，每个点可以通过加减来左移或右移，而每移动一次就是修改一次，所以从每个点向两边相邻的点连一条权值为一的边。
  • 最后答案就是从1到n+1的最短路径。

注意三点：

  1. 点1如果为负数要先把它变为0再连边
  2. 点1不向两边连，点2不向左连
  3. 连向的点可能超过n+1,也是合法的
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
int n,a,nn,dis[110005],u,v,ans;
queue<int>q;
bool inq[110005];
int fir[110005],nxt[660010],to[660010],w[660010],tot;
void line(int x,int y,int z)
{
	nxt[++tot]=fir[x];
	fir[x]=tot;
	to[tot]=y;
	w[tot]=z;
}
void spfa()
{
	memset(dis,0x3f,sizeof dis);
	dis[1]=0;
	q.push(1);
	while(!q.empty())
	{
		u=q.front();
		for(int i=fir[u];i;i=nxt[i])
		{
			v=to[i];
			if(dis[v]>dis[u]+w[i])
			{
				dis[v]=dis[u]+w[i];
				if(!inq[v]) inq[v]=1,q.push(v);
			}
		}
		inq[u]=0;
		q.pop();
	}
}
int main()
{
	scanf("%d",&n);
	scanf("%d",&a);
	if(a<0) ans=-a,line(1,2,0);
	else line(1,a+2,0);
	nn=max(nn,a+2);
	for(int i=2;i<=n;i++)
	{
		scanf("%d",&a);
		if(a>=0) line(i,i+a+1,0);
		nn=max(nn,i+a+1);
		line(i,i+1,1);
		line(i+1,i,1);
	}
	for(int i=n+1;i<nn;i++) line(i,i+1,1),line(i+1,i,1);
	spfa();
	printf("%d",dis[n+1]+ans);
}

法二 DP的（优先队列）堆优化

设f[i]表示把前i个数合法化需要的最小操作次数,则：
f[i]=f[j]+abs(a[j+1]-(i-j-1));{0<=j< i, 1<=i<=n;}

这里的abs(a[j+1]+j+1-i)便是整个题最磨人的地方了。

尽管我们很容易想到分类讨论，但却发现两个严肃的问题：

  • a[j+1]可能为负数，所以a[j+1]+j+1并不具有严格的增减性，不能使用单调队列优化

  • 当i< a[j+1]+j+1时， 我们不知道什么时候该pop,因为i在递增，此时不满足的a[j+1]+j+1可能以后就满足了。

so，我们采用了一种巧妙的优化。见代码：

#include<cstdio>
#include<queue>
#include<cctype>
using namespace std;
struct node
{
	int x,y;
	node(){}
	node(int a,int b){x=a;y=b;}
	bool operator < (const node& p)const{return x>p.x;}
};
priority_queue<node>q;
int n,a[100005],minx=10000000,f[100005];
inline void get(int &a)
{
	char c;a=0;int f=1;
	while(!isdigit(c=getchar())) if(c=='-') f=-1;
	while(isdigit(c)) a=a*10+c-'0',c=getchar();
	a*=f;
}
int main()
{
	get(n);
	for(int i=1;i<=n;i++) get(a[i]);
	q.push(node(a[1]+1,a[1]+1));
	for(int i=1;i<=n;i++)
	{
		while(!q.empty()&&q.top().y<=i)
		{
			minx=min(q.top().x-2*q.top().y,minx);
			q.pop();
		}
		f[i]=minx+i;
		if(!q.empty()) f[i]=min(f[i],q.top().x-i);
		q.push(node(f[i]+a[i+1]+i+1,a[i+1]+i+1));
	}
	printf("%d",f[n]);
}

pay attention to the rule!!!
这个优先队列并不是按照a[i+1]+i+1来排序的，而是f[i]+a[i+1]+i+1!!!
首先，对于a[j+1]+j+1<=i的情况，因为前面的i满足，后面的i肯定也满足，所以直接用一个minx来储存，然后直接pop。
而对于a[j+1]+j+1>i的情况，直接取最优的f[j]+a[j+1]+j+1则可。

但是，同志们可能会发现一个问题，既然是按照f[i]+a[i+1]+i+1排序，
那么如果在a[j+1]+j+1>i时直接退出，那么如果后面还有满足a[j+1]+j+1<=i的值怎么办呢？

这就是巧妙之处了。

令f[j]=a,a[j+1]+j+1=b；当b>i时记作a1,b1;当b<=i时记作a2,b2
那么我们的f[i]=min(a2-b2+i,a1+b1-i);
这里我们取了最优的a1+b1-i,那么我们如何说明后面的a2-b2+i不会比现在更优呢？
注意优先队列的规则，我们得知：a1+b1<=a2+b2
而我们已知：b2<=i
所以 2 * b2 <= 2 * i
所以 b2 <= 2*i-b2
所以 a2+b2<=a2+2*i-b2
所以 a1+b1<=a2+b2<=a2-b2+2*i
所以 a1+b1-i<=a2-b2+i !!!
证明至此，我们可知后面的解不会比当前更优，直接退出。

法三：DP的线段树优化

众所周知，对于这种含绝对值的DP,线段树是最在行（难写）的了，（基本上就没写对过）。。。。鉴于本人目前对于线段树优化一知半解。。
这里。。。就不加解析了。。。等孤学成归来在好好补上吧。。这里就直接上代码了

于是过了两年我学成归来了，做到原题根本没有想到上面两种做法直接开打线段树优化。。
f [ i ] = f [ j ] + a b s ( a [ j + 1 ] − ( i − j − 1 ) ) ; 0 < = j < i , 1 < = i < = n ; f[i]=f[j]+abs(a[j+1]-(i-j-1));{0<=j< i, 1<=i<=n;} f[i]=f[j]+abs(a[j+1]−(i−j−1));0<=j<i,1<=i<=n;
根据 i i i可以得到abs取正负时 j j j的范围，可以看出时两个区间，线段树维护 f [ j ] + a [ j + 1 ] + j + 1 f[j]+a[j+1]+j+1 f[j]+a[j+1]+j+1和 f [ j ] − ( a [ j + 1 ] + j + 1 ) f[j]-(a[j+1]+j+1) f[j]−(a[j+1]+j+1)的最小值即可。

这个代码是从后往前做的，大致思路相同。

#include<bits/stdc++.h>
#define maxn 100005
using namespace std;
const int inf = 1<<30;
int n,a[maxn],f[maxn];
struct Segment_Tree{
	int mn[maxn<<2];
	void insert(int i,int l,int r,int x,int v){
		if(l==r) {mn[i]=v;return;}//special
		int mid=(l+r)>>1;
		if(x<=mid) insert(i<<1,l,mid,x,v);
		else insert(i<<1|1,mid+1,r,x,v);
		mn[i]=min(mn[i<<1],mn[i<<1|1]);
	}
	int query(int i,int l,int r,int x,int y){
		if(x<=l&&r<=y) return mn[i];
		int mid=(l+r)>>1,ret=inf;
		if(x<=mid) ret=min(ret,query(i<<1,l,mid,x,y));
		if(y>mid) ret=min(ret,query(i<<1|1,mid+1,r,x,y));
		return ret;
	}
	void ins(int x,int v){insert(1,1,n+1,x,v);}
	int ask(int l,int r){return query(1,1,n+1,l,r);}
}t1,t2;
int main()
{
	freopen("data.in","r",stdin);
	freopen("data.out","w",stdout);
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	t1.ins(n+1,-(n+1)),t2.ins(n+1,n+1);
	for(int i=n;i>=1;i--){
		f[i]=inf;
		if(a[i]>=0) f[i]=min(f[i],t1.ask(i+1,min(i+a[i]+1,n+1))+i+a[i]+1);
		if(i+a[i]+1<=n) f[i]=min(f[i],t2.ask(max(i+a[i]+2,i+1),n+1)-(i+a[i]+1));
		t1.ins(i,f[i]-i),t2.ins(i,f[i]+i);
	}
	printf("%d\n",f[1]);
}

Just do it！

Master.Yi
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
2015重庆市NOIP模拟赛题目+数据
10-04
2015年10月4日重庆市NOIP模拟赛题目+数据
【DP+线段树优化】[CQBZOJ2933]数据
weixin_30408165的博客
10-11 106
题目描述 Mr_H 出了一道信息学竞赛题，就是给 n 个数排序。输入格式是这样的： 试题有若干组数据。每组数据的第一个是一个整数 n，表示总共有 n 个数待排序；接下来 n 个 整数，分别表示这 n 个待排序的数。 例如：3 4 2 –1 4 1 2 3 4，就表示有两组数据。第一组有 3 个数(4,2,-1)，第二组有 4 个数(1,2,3,4)。可是现在 Mr_H...
【重庆市NOIP模拟赛】业务
SunIsMe的博客
08-28 1071
题目描述 Mr_H 谋得一份兼职——货车司机，从此以后他将会开着货车穿行在 C 国的各大城市之间。 C 国中有 n 座城市（编号为 1~n），并且有 m 条双向公路，每条公路连接两座不同的城市。货车 从任意一座城市出发都可以抵达任意另一座城市。在每条公路上，都有一个收费站，通过的车辆需要 交纳一定过路费。可能有多条公路连接相同的两座城市。 为了增加财政收入，C 国还在每座城市也设
数据--dp，线段树优化
weixin_30276935的博客
10-11 113
题目描述 Mr_H 出了一道信息学竞赛题，就是给 n 个数排序。输入格式是这样的：试题有若干组数据。每组数据的第一个是一个整数 n，表示总共有 n 个数待排序；接下来 n 个整数，分别表示这 n 个待排序的数。 例如：3 4 2 –1 4 1 2 3 4，就表示有两组数据。第一组有 3 个数(4,2,-1)，第二组有 4个数(1,2,3,4)。可是现在 Mr_H 做的输入数...
[BZOJ2933]数据
UISG103的博客
07-24 370
题目描述 Mr_H 出了一道信息学竞赛题，就是给 n 个数排序。输入格式是这样的： 试题有若干组数据。每组数据的第一个是一个整数 n，表示总共有 n 个数待排序；接下来 n 个 整数，分别表示这 n 个待排序的数。 例如：3 4 2 –1 4 1 2 3 4，就表示有两组数据。第一组有 3 个数(4,2,-1)，第二组有 4 个数(1,2,3,4)。可是现在 Mr_H 做的输入数据出了一
NOIP模拟赛 试题
07-15
此为NOIP模拟赛的试题。 共三题，时间2小时； 难度符合NOIP标准 （还是自己出的呢） 勘误： 1.A题 输入样例中两字符串无空格 2.C题 n是位置的个数，不是猫的个数
noip模拟赛day1
10-28
【NOIP模拟赛Day1详解】——为竞赛之路铺砖添瓦 全国青少年信息学奥林匹克联赛（NOIP）是一项旨在发掘和培养青少年编程才能的重要赛事。对于参赛者来说，赛前的充分准备至关重要，尤其是模拟赛的参与，可以帮助考生...
noip模拟赛第二试
10-27
【标题】"noip模拟赛第二试"涉及的是一个针对全国青少年信息学奥林匹克竞赛（NOIP）的模拟测试环节。NOIP是中国计算机学会举办的一项面向中学生的编程竞赛，旨在培养和选拔优秀的计算机科学后备人才。模拟赛是为正式...
NOIP2017提高组数据
10-25
NOIP2017提高组数据，非常好的程序设计竞赛资源，值得下载
NOIP近7年真题（含测试数据和测试软件）
02-26
全国青少年信息学奥林匹克联赛（National Olympiad in Informatics in Provinces，简称NOIP）自1995年至2017年已举办23次。每年由中国计算机学会统一组织。 NOIP在同一时间、不同地点以各省市为单位由特派员组织。全国统一大纲、统一试卷。初、高中或其他中等专业学校的学生可报名参加联赛。联赛分初赛和复赛两个阶段。初赛考察通用和实用的计算机科学知识，以笔试形式进行。复赛为程序设计，须在计算机上调试完成。参加初赛者须达到一定分数线后才有资格参加复赛。联赛分普及组和提高组两个组别，难度不同，分别面向初中和高中阶段的学生。
noip往年题目和数据
03-21
NOIP 95---12年数据还有试题和题解，十分完善
【总结】20151031重庆市NOIP模拟赛
蒟蒻的奋斗史
10-31 1385
世上最最悲催的事莫过于手残打错字母… 考完之后还是一句话我是SB… 今日的题其实都想到了正解，然而就木有AC一道题… 第一题当时想了很久，不知道为什么会对这个模型如此生疏（可能是好久都木有做组合数学的题了），还是在写第三题的时候才想出来。一阵狂喜就手残打错字母了……
#bzoj2934#【重庆市NOIP模拟赛】业务（SPFA / Dijk）
ElaineZXY的博客
07-24 304
2934: 【重庆市NOIP模拟赛】业务 时间限制: 1 Sec  内存限制: 128 MB 题目描述 Mr_H 谋得一份兼职——货车司机，从此以后他将会开着货车穿行在 C 国的各大城市之间。 C 国中有 n 座城市（编号为 1~n），并且有 m 条双向公路，每条公路连接两座不同的城市。货车 从任意一座城市出发都可以抵达任意另一座城市。在每条公路上，都有一个收费站，通过的车辆需要
【jzoj】2018.2.7NOIP普及组——某【B&C】组模拟赛
weixin_30871905的博客
02-08 125
前言 ……终于改完了，像之前小L一样崩溃。今天C组和B组一起做题，所以…… 正题 题目1：教主的花园（jzoj1792） 一平面直角坐标系，在x轴的位置建立一堵墙，墙上有n道门，给出门的位置，询问两个坐标的距离（曼哈顿）。 输入输出（建议无视） Input 输入的第1行为一个正整数N，为屏障上入口的个数。 第2行有N个整数，...
noip 2018 模拟赛13
Caristra的博客
11-01 268
T1T_1T1​——nim(3833) Description: 有nnn堆石子，第iii堆的石子个数为iii，两人可以每次可以在其中一堆中拿任意数量的石子，无法再继续拿为败。求必胜的方案数。 n≤101000n\le 10^{1000}n≤101000 Solution: 由于出题人的疯狂，强行练一练高精度， 此题可以说是博弈题，也可以说是找规律题（你甚至可以在大样例中找规律） 具体的规律是：...
9.22
weixin_30756499的博客
09-22 61
1.字典实例：建立学生学号成绩字典，做增删改查遍历操作。 #字典实例：建立学生学号成绩字典，做增删改查遍历操作。 cj={'39':89,'33':90,'38':80,'40':79} cj['36']=78#增加 print(cj) cj.pop('39')#删除 print(cj) cj['40']=69#修改值 print(cj) print(cj.get(...
【NOIP2017模拟赛】思维+转化+图论 徒然Children（好题）
Starlight的博客
11-09 751
题目描述 全世界都在谈恋爱，只有我深爱着学习。——前言。 2333身为单身狗却意外喜欢的吃狗粮的yuki坚信着：“自古红蓝出cp，黑白天生是夫妻，最是销魂红绿配，天然金紫成双对，生死相随红与黄，白紫一逢春满堂，千里缘牵白与绿，红黑生来为相聚。”的cp配对原则，却被徒然children中绿黄的御宅的春天组（山根隆夫×栗原千代）莫名戳中萌点。可见死宅真是一个有前途的职业。在本季末尾，一群cp们正在举
【DP】 花店橱窗布置
Master.Yi的博客
11-05 1647
花店橱窗布置(flower.cpp) 时间限制: 1 Sec  内存限制: 64 MB 题目描述 假设你想以最美观的方式布置花店的橱窗。你有F束花，每束花的品种都不一样，同时，你至少有同样数量的花瓶，被按顺序摆成一行。花瓶的位置是固定的，并从左至右，从1至V顺序编号，V是花瓶的数目，编号为1的花瓶在最左边，编号为V的花瓶在最右边。花束则可以移动，并且每束花用1至F的整数唯一标识。标识
noip模拟赛 测试数据
最新发布
09-13
NOI（全国青少年信息学奥林匹克竞赛）模拟赛的测试数据是指用于评测参赛选手的程序的输入和对应的输出。测试数据是非常重要的，因为它决定了参赛选手的程序能否正确地解决问题。 在NOI模拟赛中，测试数据具有以下...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值