【差分约束系统】 note

>>>>题目

【题目描述】

最近有一款很火的游戏,叫做八分音符酱,它和马里奥很相似,不过它的跳跃距离是由你的声音大小来控制
的。不过我们现在对玩法就行一些修改:
现有一共有n 个柱子,两个相邻的柱子之间的初始水平距离为1,蠢蠢的jyb 现在在最矮的柱子上,他每次
只能向恰好比这个柱子高的另一个柱子跳跃,最后要跳到最高的柱子上。

jyb 需要从第二个柱子跳到第一个柱子,再跳到第三个柱子。
jyb 的最大声音为d,代表他能在满足 |posi-posj| + |heighti-heightj | ≤d  的两个柱子之间跳跃。

假设我们可以在不改变它们位置的相对顺序的前提下水平移动柱子,调整他们的水平位置(但相邻间隔至少为1,且为整数)。

jyb 想问你:在能从最矮柱子跳到最高柱子的前提下,它们(最矮柱子跳到最高柱子) 的最大水平距离是多少。

【输入描述】

第1 行,1 个整数T, 表示数据组数,对于每组数据:
第1 行,1 个整数n;m,表示一共有多少个柱子。
接下来1 行,有n 个数,hi 表示柱子的高度。保证柱子高度互不相同

【输出描述】

对于每组数据,输出在能从最矮柱子跳到最高柱子的前提下,它们的最大水平距离是多少。如果不能,输出-1

【输入样例1】

2
3 4
3 2 4
3 4
3 2 6

【输出样例1】

2

-1

【输入样例2】

2
5 10
4 2 1 8 10
5 2
10 8 2 1 4

【输出样例】

12

-1

【数据范围及约定】

对于30% 的数据,1 n 100,1 d 1000,1 hi 1000;
•对于100% 的数据,1 T 100,1 n 103,1 d 106,1 hi 106。

>>>>分析

我们不能随意地将水平位置拉开,因为必须要能够从最矮的柱子跳到最高的柱子。

其实就相当于有了很多限制条件,显然这是一个差分约束系统。

先根据柱子的高度排序(初始化时用结构体记录高度和位置)

我们约定,在相邻两柱子,从下标小的向大的建边(以下用u(下标小),v(下标大)表示)

我们有两个约束条件:

1..根据题意:pos[v]-pos[u]+height[v]-height[u]≤d

2..两个相邻的柱子之间距离最少为1::pos[v]-pos[u]≥1

整理后得到

pos[v]-pos[u]≤d-(height[v]-height[u])

于是建边分为两种:

1.从i 到i+1建边,长度为 d-(height[v]-height[u])

2.从i+1到i建边,长度为-1

>>>>代码

#include<bits/stdc++.h>
#define maxn 1005
using namespace std;
struct node
{
int to,nxt,w;
}e[maxn<<];
struct node1
{
int pos,h;
}a[maxn];
int T,n,d,tot,flag,st,ed;
int dis[maxn],head[maxn],in[maxn];
bool vis[maxn];
queue<int> q;
bool cmp(const node1 &x,const node1 &y)
{
return x.h<y.h;
}
void add(int u,int v,int w)
{
e[++tot].to=v;
e[tot].nxt=head[u];
e[tot].w=w;
head[u]=tot;
}
int spfa(int st,int ed)//找最短路
{
memset(dis,,sizeof(dis));
memset(vis,,sizeof(vis));
memset(in,,sizeof(in));
q.push(st),vis[st]=;
in[st]++,dis[st]=;
while(!q.empty())
{
int now=q.front();
q.pop(),vis[now]=;
for(int i=head[now];i;i=e[i].nxt)
{
int v=e[i].to;
if(dis[v]>dis[now]+e[i].w)
{
dis[v]=dis[now]+e[i].w;
if(!vis[v])
{
vis[v]=;
q.push(v);
in[v]++;
if(in[v]>=n+) return -;//判负环
}
}
}
}
return dis[ed];
}
int main()
{
// freopen("note.in","r",stdin);
// freopen("note.out","w",stdout);
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&n,&d);
for(int i=;i<=n;++i)
{
scanf("%d",&a[i].h);
a[i].pos=i;
}
tot=,flag=;
memset(head,,sizeof(head));
sort(a+,a+n+,cmp);
//找到起点和终点
st=min(a[].pos,a[n].pos);
ed=max(a[].pos,a[n].pos);
for(int i=;i<=n-;++i)
{
add(i+,i,-);
//从标号小的边向大的边连
int u=min(a[i].pos,a[i+].pos);
int v=max(a[i].pos,a[i+].pos);
if(d<v-u+a[i+].h-a[i].h)
{
flag=; break;
}
add(u,v,d-a[i+].h+a[i].h);
}
if(flag) printf("%d\n",spfa(st,ed));
else printf("-1\n");
}
}
/*
2
3 4
3 2 4
3 4
3 2 6
*/

完结撒花!

【差分约束系统】 note的更多相关文章

  1. 【差分约束系统】【spfa】UVALive - 4885 - Task

    差分约束系统讲解看这里:http://blog.csdn.net/xuezhongfenfei/article/details/8685313 模板题,不多说.要注意的一点是!!!对于带有within ...

  2. 【POJ 1275】 Cashier Employment(差分约束系统的建立和求解)

    [POJ 1275] Cashier Employment(差分约束系统的建立和求解) Cashier Employment Time Limit: 1000MS   Memory Limit: 10 ...

  3. UVA11478 Halum [差分约束系统]

    https://vjudge.net/problem/UVA-11478 给定一个有向图,每条边都有一个权值.每次你可以选择一个结点v和一个整数d,把所有以v为终点的边的权值减小d,把所有以v为起点的 ...

  4. BZOJ 2330: [SCOI2011]糖果 [差分约束系统] 【学习笔记】

    2330: [SCOI2011]糖果 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 5395  Solved: 1750[Submit][Status ...

  5. ACM/ICPC 之 差分约束系统两道(ZOJ2770-POJ1201)

    当对问题建立数学模型后,发现其是一个差分方程组,那么问题可以转换为最短路问题,一下分别选用Bellmanford-SPFA解题 ZOJ2770-Burn the Linked Camp //差分约束方 ...

  6. POJ1201 Intervals(差分约束系统)

    与ZOJ2770一个建模方式,前缀和当作点. 对于每个区间[a,b]有这么个条件,Sa-Sb-1>=c,然后我就那样连边WA了好几次. 后来偷看数据才想到这题还有两个隐藏的约束条件. 这题前缀和 ...

  7. UVA 11374 Halum (差分约束系统,最短路)

    题意:给定一个带权有向图,每次你可以选择一个结点v 和整数d ,把所有以v为终点的边权值减少d,把所有以v为起点的边权值增加d,最后要让所有的边权值为正,且尽量大.若无解,输出结果.若可无限大,输出结 ...

  8. Burn the Linked Camp(bellman 差分约束系统)

    Burn the Linked Camp Time Limit: 2 Seconds      Memory Limit: 65536 KB It is well known that, in the ...

  9. zoj 2770 Burn the Linked Camp (差分约束系统)

    // 差分约束系统// 火烧连营 // n个点 m条边 每天边约束i到j这些军营的人数 n个兵营都有容量// Si表示前i个军营的总数 那么 1.Si-S(i-1)<=C[i] 这里 建边(i- ...

随机推荐

  1. JavaScript中创建对象的几种模式

    代码如下: <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF ...

  2. 网易blog服务器关闭了,我曾经的大部分博文迁移至此,留作纪念。欢迎大家去我的简书blog~

    ARM时代过去了,另一个AI时代正在上演~ 这一次我看上了计算机视觉,一个依旧在发展,论文数量不断上升的技术领域,还有很多待研究的技术突破点,我可以玩的时间长一点了. 活到老,学到了.学以致用乐趣无穷 ...

  3. ABP入门系列之2——ABP模板项目

    进入官网下载模板项目 依次按下图选择: 输入验证码开始下载 下载提示: 二.启动项目 使用VS2017打开项目,还原Nuget包: 设置以Web结尾的项目,设置为启动项目: 打开Web.config, ...

  4. vue 脚手架关于路由的一点理解

    https://router.vuejs.org/zh/ 可以先翻翻文档看看介绍啊,一般我不怎么喜欢看文档,都是直接看人家案例,在回头看文档的,所以学习速度慢很多,希望我以后成为一个爱学习的妹子,比较 ...

  5. Max Sum Plus Plus HDU - 1024

    Max Sum Plus Plus     HDU - 1024 Now I think you have got an AC in Ignatius.L's "Max Sum" ...

  6. Vagrant Box下载缓慢解决方法

    box 搜索页面:https://atlas.hashicorp.com/boxes/search example:homestead 1,选中box和版本,先根据提示安装 2,获取box下载地址,采 ...

  7. zabbix_agentd重装后启动时IPC和共享内存段问题

    zabbix_agentd不知为啥被干掉后重装了zabbix,zabbix用户组id也变了. 重装zabbix后导致zabbix_agentd无法启动,两个问题 问题1: zabbix_agentd ...

  8. Laravel5.5学习笔记

    安装composer 下载安装脚本 php -r "copy('https://install.phpcomposer.com/installer', 'composer-setup.php ...

  9. 【Django视图与网址进阶004】

    一.在网页上做加减法 1. 采用 /add/?a=4&b=5 这样GET方法进行 django-admin.py startproject zqxt_views cd zqxt_views p ...

  10. linux准备工具VMware以及ubuntu

    链接:https://pan.baidu.com/s/1j7DZCmzOcC2rPAciHtZgkA 密码:et5s