POJ 1384 Intervals (区间差分约束,根据不等式建图,然后跑spfa)
传送门:
http://acm.hdu.edu.cn/showproblem.php?pid=1384
Intervals
Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 4841 Accepted Submission(s): 1815
Write a program that:
> reads the number of intervals, their endpoints and integers c1, ..., cn from the standard input,
> computes the minimal size of a set Z of integers which has at least ci common elements with interval [ai, bi], for each i = 1, 2, ..., n,
> writes the answer to the standard output
Process to the end of file.
3 7 3
8 10 3
6 8 1
1 3 1
10 11 1
#include <iostream>
#include <cstdio>
#include<stdio.h>
#include<algorithm>
#include<cstring>
#include<math.h>
#include<memory>
#include<queue>
#include<vector>
using namespace std;
#define max_v 50010
#define INF 9999999
int tot;
int head[max_v];
int vis[max_v];
int dis[max_v];
int minn,maxx;
struct node
{
int u,v,val;
int next;
}Edge[max_v<<];
void addEdge(int u,int v,int val)
{
Edge[tot].u=u;
Edge[tot].v=v;
Edge[tot].val=val;
Edge[tot].next=head[u];
head[u]=tot++;
}
void spfa()
{
for(int i=minn-;i<=maxx;i++)
dis[i]=-INF;
queue<int> q;
dis[maxx]=;
vis[maxx]=;
q.push(maxx);
while(!q.empty())
{
int u=q.front();
q.pop();
vis[u]=;
for(int i=head[u];i!=-;i=Edge[i].next)
{
int v=Edge[i].v;
if(dis[v]<dis[u]+Edge[i].val)
{
dis[v]=dis[u]+Edge[i].val;
if(!vis[v])
{
vis[v]=;
q.push(v);
}
}
}
}
printf("%d\n",dis[minn-]);
return ;
}
int main()
{
int n,a,b,c;
while(~scanf("%d",&n))
{
tot=;
memset(head,-,sizeof(head));
memset(vis,,sizeof(vis));
maxx=;
minn=INF;
for(int i=;i<n;i++)
{
scanf("%d %d %d",&a,&b,&c);
a++;
b++;
minn=min(minn,a);
maxx=max(maxx,b);
addEdge(b,a-,c);
}
for(int i=minn;i<=maxx;i++)
{
addEdge(i,i-,);
addEdge(i-,i,-);
}
spfa();
}
return ;
}
今天学差分约束的时候又遇到了这个题,跑最长路初始化的时候搞错了,应该是dis[i]=-INF。。。。,卡了好一会
从一系列区间[a,b]中每个至少取ci个数构成集合s,问你集合s至少需要多少个数
区间约束问题,差分约束解决
f(x):[0,x]区间内取f(x)个数
隐藏关系:1>=f(i)-f(i-1)>=0,变形一下:
f(i)-f(i-1)>=0
f(i-1)-f(i)>=-1
maxx:所有区间的最大值
minx:所有区间的最小值
f(b)-f(a-1)>=c
f(i)-f(i-1)>=0 i属于[minx,maxx]
f(i-1)-f(i)>=-1 i属于[minx,maxx]
按照j->i建边,权值为c
所以:
(a-1)->b 建边 权值c
(i-1)->i 建边 权值0
i->(i-1) 建边 权值-1
>=是最小值,所以是最长路
所以是问你从minx-1到maxx的最长路
建图完毕之后spfa跑一遍最长路,不能用dj,因为存在负权边
#include<stdio.h>
#include<iostream>
#include<math.h>
#include<string.h>
#include<set>
#include<map>
#include<list>
#include<math.h>
#include<queue>
#include<algorithm>
using namespace std;
typedef long long LL;
#define INF 9999999999
#define me(a,x) memset(a,x,sizeof(a))
int mon1[]= {,,,,,,,,,,,,};
int mon2[]= {,,,,,,,,,,,,};
int dir[][]= {{,},{,-},{,},{-,}}; int getval()
{
int ret();
char c;
while((c=getchar())==' '||c=='\n'||c=='\r');
ret=c-'';
while((c=getchar())!=' '&&c!='\n'&&c!='\r')
ret=ret*+c-'';
return ret;
}
void out(int a)
{
if(a>)
out(a/);
putchar(a%+'');
} #define max_v 50005
int dis[max_v];
int vis[max_v]; struct node
{
int v,w;
node(int vv=,int ww=):v(vv),w(ww){}
};
vector<node> G[max_v];
int n,m; void init()
{
for(int i=;i<max_v;i++)
G[i].clear();
for(int i=;i<max_v;i++)
{
dis[i]=-INF;
vis[i]=;
}
} void spfa(int s)
{ queue<int> q;
q.push(s);
vis[s]=;
dis[s]=; while(!q.empty())
{
int u=q.front();
q.pop();
vis[u]=; for(int j=;j<G[u].size();j++)
{
int v=G[u][j].v;
int w=G[u][j].w;
if(dis[v]<dis[u]+w)//最长路
{
dis[v]=dis[u]+w;
if(vis[v]==)
{
vis[v]=;
q.push(v);
}
}
}
}
} int f(int u,int v)
{
for(int j=;j<G[u].size();j++)
{
if(G[u][j].v==v)
return ;
}
return ;
}
int main()
{
int a;
while(~scanf("%d",&a))
{
int x,y,w;
init();
int maxx=-INF,minx=INF;
for(int i=;i<=a;i++)
{
scanf("%d %d %d",&x,&y,&w);
x++,y++;//防止减1的时候出现负数
maxx=max(maxx,y);
minx=min(minx,x);
if(f(x-,y))//预防重边
G[x-].push_back(node(y,w));
}
for(int i=minx;i<=maxx;i++)
{
if(f(i-,i))//预防重边
G[i-].push_back(node(i,));
if(f(i,i-))//预防重边
G[i].push_back(node(i-,-));
}
spfa(minx-);
printf("%d\n",dis[maxx]);
}
return ;
}
/*
题目意思:
从一系列区间[a,b]中每个至少取ci个数构成集合s,问你集合s至少需要多少个数
区间约束问题,差分约束解决 分析:
f(x):[0,x]区间内取f(x)个数 则区间[a,b]内至少取c个数可以变形为:f(b)-f(a-1)>=c
隐藏关系:1>=f(i)-f(i-1)>=0,变形一下:
f(i)-f(i-1)>=0
f(i-1)-f(i)>=-1
maxx:所有区间的最大值
minx:所有区间的最小值 i属于[minx,maxx] 总结一下得到这样三个关系:
f(b)-f(a-1)>=c
f(i)-f(i-1)>=0 i属于[minx,maxx]
f(i-1)-f(i)>=-1 i属于[minx,maxx] 这三个表达式的形式都是一样的,xi-xj>=c
按照j->i建边,权值为c
所以:
(a-1)->b 建边 权值c
(i-1)->i 建边 权值0
i->(i-1) 建边 权值-1 我们要求的东西ans可以变形为:f(maxx)-f(minx-1)>=ans
>=是最小值,所以是最长路
所以是问你从minx-1到maxx的最长路
建图完毕之后spfa跑一遍最长路,不能用dj,因为存在负权边 */
POJ 1384 Intervals (区间差分约束,根据不等式建图,然后跑spfa)的更多相关文章
- POJ 1201 Intervals(差分约束 区间约束模版)
关于差分约束详情可阅读:http://www.cppblog.com/menjitianya/archive/2015/11/19/212292.html 题意: 给定n个区间[L,R], 每个区间至 ...
- poj 1201 Intervals(差分约束)
做的第一道差分约束的题目,思考了一天,终于把差分约束弄懂了O(∩_∩)O哈哈~ 题意(略坑):三元组{ai,bi,ci},表示区间[ai,bi]上至少要有ci个数字相同,其实就是说,在区间[0,500 ...
- POJ 1201 Intervals (经典) (差分约束)
<题目链接> 题目大意:给你$n$段区间,$a_i,b_i,c_i$ 表示在 $[a_i,b_i]$ 区间内至少要选择$c_i$个点.现在问你在满足这n个条件的情况下,最少要选多少个点? ...
- HDU 1384 Intervals【差分约束-SPFA】
类型:给出一些形如a−b<=k的不等式(或a−b>=k或a−b<k或a−b>k等),问是否有解[是否有负环]或求差的极值[最短/长路径].例子:b−a<=k1,c−b&l ...
- hdu 1384 Intervals (差分约束)
Problem - 1384 好歹用了一天,也算是看懂了差分约束的原理,做出第一条查分约束了. 题意是告诉你一些区间中最少有多少元素,最少需要多少个元素才能满足所有要求. 构图的方法是,(a)-> ...
- POJ 1201 Intervals【差分约束】
传送门:http://poj.org/problem?id=1201 题意: 有n个如下形式的条件:,表示在区间[, ]内至少要选择个整数点.问你满足以上所有条件,最少需要选多少个点? 思路:第一道差 ...
- 【题解】 POJ 1201 Intervals(差分约束)
懒得复制,戳我戳我 Solution: 这道题就是一个板子题 抽象成第\(a\)至第\(b\)间选择数的个数为\(c\),我们就可以用前缀和来表示,这样就可以得到不等式\(s[b]-s[a-1]> ...
- POJ 1201 Intervals (差分约束,最短路)
题意: 有一个集合Z,其元素都是整整数,但是数量未知.现有n个约束,形如 [a,b]=c 表示整数区间[a,b]中有c个元素在Z中出现.问集合Z最小可能含多少个元素? 思路: 对于所给的区间 cnt[ ...
- poj 1201 Intervals【差分约束+spfa】
设s为前缀和,首先显然的条件是\[ s_{bi}-s_{ai-1}>=c \],然后隐含的是\[ s_i-s_{i-1}>=0 s_i-s_{i-1}<=1 \] 然后根据差分约束, ...
- HDU 1384 Intervals(差分约束)
Intervals Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total ...
随机推荐
- LeetCode ImplementStrstr
class Solution { public: char *strStr(char *haystack, char *needle) { if (haystack == NULL || needle ...
- PhpStorm 自定义快捷键
PhpStorm 两个重要快捷键 1.CTRL+SHIFT+A 用于恢复隐藏项 2.文件之间的快速跳转 CTRL+SHIFT+N 3.自定义快捷键 第一步:打开左上角file-setting 第二步: ...
- ActiveMQ的用途
ActiveMQ 是Apache出品,最流行的,能力强劲的开源消息总线. ActiveMQ 是一个完全支持JMS1.1和J2EE 1.4规范的 JMS Provider实现. 消息队列的主要作用是为了 ...
- drupal7 转化 public:// 为实际url
file_create_url('public://xxx.png'); // 得到URL drupal_realpath('public://xxx.png'); // 得到系统路径(磁盘路径,如D ...
- JavaScript的进阶之路(四)理解对象2
对象的三个属性 原型属性 1.var v={}的原型是Object.prototype;继承了一个constructor属性指代Object()构造函数,实际的原型是constructor.proto ...
- C# Timer定时器用法
System.Timers.Timer timer1 = new System.Timers.Timer(); timer1.Elapsed += new System.Timers.ElapsedE ...
- Linux / mysql: is it safe to copy mysql db files with cp command from one db to another?
Copying is very simple for MyISAM and completely 100% risky (near suicidal) with InnoDB. From your q ...
- 131.005 Unsupervised Learning - Cluster | 非监督学习 - 聚类
@(131 - Machine Learning | 机器学习) 零. Goal How Unsupervised Learning fills in that model gap from the ...
- spring boot(7)-mybatis全注解化
关于配置数据库可以参考上一篇文章,这里只讲mybatis pom.xml <!-- 引入mybatis --> <dependency> <groupId>org. ...
- Windows server 2008系统各类版本的优缺点比较,Windows2008系统标准版 企业版 数据中心版 WEB版等
大家都知道Windows Server 2008 发行了多种版本,以支持各种规模的企业对服务器不断变化的需求.Windows Server 2008 有 5 种不同版本,另外还有三个不支持 Windo ...