POJ-1201-Intervals，POJ-1716-Integer Intervals（区间约束）

题目链接：http://poj.org/problem?id=1201

题目大意：给出n个区间，每个区间[a,b]能最少能够选择c个数，问，最后能选多少个数，使这些数对所有的区间都满足。

思路：区间约束的一道例题了，但是我感觉我写的好玄学？？麻瓜代码，看别人是最长路过的，我怎么还是最短路？？反向最短路谜之AC？还只用了200+ms？？搞得我好迷啊，建图就是差分约束条件建图。

对于所有的点，首先满足：0<=d[i]-d[i-1]<=1，（选不选该点）

然后对于每个区间，满足d[r]-d[l-1]>=c，然后*-1得：d[l-1]-d[r]<=-c,（区间选择一些点）

得出约束条件：d[i-1]-d[i]<=0；d[i]-d[i-1]<=1；d[l-1]-d[r]<=-c，然后就直接跑板子了。反向最短路强行AC！！？？

补充，后来发现最长路就是反向建边的最短路，所以这个就是歪打正着？

AC：

//#pragma comment(linker, "/STACK:1024000000,1024000000")
  
#include<stdio.h>
#include<string.h> 
#include<math.h> 
   
#include<map>  
//#include<set>
#include<deque> 
#include<queue> 
#include<stack> 
#include<bitset>
#include<string> 
#include<fstream>
#include<iostream> 
#include<algorithm> 
using namespace std; 
  
#define ll long long 
#define Pair pair<int,int>
//#define max(a,b) (a)>(b)?(a):(b)
//#define min(a,b) (a)<(b)?(a):(b)
#define clean(a,b) memset(a,b,sizeof(a))// 水印
//std::ios::sync_with_stdio(false);
//  register
const int MAXN=5e4+10;
const int INF32=0x3f3f3f3f;
const ll INF64=0x3f3f3f3f3f3f3f3f;
const ll mod=1e9+7;
const double PI=acos(-1.0);

struct node{
	int v,val,nxt;
	node(int _v=0,int _val=0,int _nxt=0){
		v=_v;val=_val;nxt=_nxt;
	}
}edge[MAXN<<2];
int head[MAXN],ecnt;
int st[MAXN],tot;
int dis[MAXN],vis[MAXN],sum[MAXN];
int n,flag;

void intt(){
	clean(head,-1);
	ecnt=0;tot=0;
}
void add(int u,int v,int val){
	edge[ecnt]=node(v,val,head[u]);
	head[u]=ecnt++;
}
void dfs_spfa(int str){
	clean(dis,INF32);dis[str]=0;
	clean(vis,0);vis[str]=1;
	clean(sum,0);sum[str]++;
	st[tot++]=str;
	while(tot){
		int u=st[--tot];vis[u]=0;
		//printf("\n u: %d - ",u);
		for(int i=head[u];i+1;i=edge[i].nxt){
			int v=edge[i].v;
			if(dis[v]>dis[u]+edge[i].val){
				//printf("v: %d ",v);
				dis[v]=dis[u]+edge[i].val;
				if(vis[v]==0){
					vis[v]=1;
					st[tot++]=v;
					sum[v]++;
					if(sum[v]>n){
						flag=1;return ;
					}
				}
			}
		}
	}
}
int main(){
	intt();
	scanf("%d",&n);
	int a,b,c,ed=-1,str=INF32;
	for(int i=1;i<=n;++i){
		scanf("%d%d%d",&a,&b,&c);
		ed=max(ed,b);str=min(str,a-1);
		add(b,a-1,-c);
	}
	for(int i=str;i<ed;++i){
		add(i,i+1,1);
		add(i+1,i,0);
	}flag=0;
	dfs_spfa(ed);
//	printf("\n");
//	for(int i=1;i<=ed;++i){
//		printf("%d ",dis[i]);
//	}printf("\n");
	printf("%d\n",abs(dis[str]));
}
/*
8
4 9 5
8 9 1
5 6 1
3 9 5
7 9 1
2 5 2
3 4 1
2 5 3

输出6。
*/

我又去讨论区看了看，里面有用线段树写的，贪心的，网络流的（神TMD网络流），强连通分量变形的，IDA*+舞蹈链过的。

还有一个直接暴力双向BFS+随机搜索（险过）？？

还有一个说如何位运算优化搜索会比较快？？什么东西嘛这都是？？！！

神仙神仙，这道题等等再看看吧。。。

最后我发现下面的评论有些LZSB，让楼主吃药什么的。。，看了一下他们的做题数，Solved，1K+。看来这个神仙也写不出来。

算了算了，神仙打架，我就不参与了，我就膜拜膜拜，也不想实现了。

update：

POJ-1716-Integer Intervals

题目链接：http://poj.org/problem?id=1716

题目大意：n个区间内，每个区间选2个数，问最少选多少个数。

思路，原题代码，-c变成-2即可。

AC：

//#pragma comment(linker, "/STACK:1024000000,1024000000")
  
#include<stdio.h>
#include<string.h> 
#include<math.h> 
   
#include<map>  
//#include<set>
#include<deque> 
#include<queue> 
#include<stack> 
#include<bitset>
#include<string> 
#include<fstream>
#include<iostream> 
#include<algorithm> 
using namespace std; 
  
#define ll long long 
#define Pair pair<int,int>
//#define max(a,b) (a)>(b)?(a):(b)
//#define min(a,b) (a)<(b)?(a):(b)
#define clean(a,b) memset(a,b,sizeof(a))// 水印
//std::ios::sync_with_stdio(false);
//  register
const int MAXN=1e4+10;
const int INF32=0x3f3f3f3f;
const ll INF64=0x3f3f3f3f3f3f3f3f;
const ll mod=1e9+7;
const double PI=acos(-1.0);

struct node{
	int v,val,nxt;
	node(int _v=0,int _val=0,int _nxt=0){
		v=_v;val=_val;nxt=_nxt;
	}
}edge[MAXN<<4];
int head[MAXN],ecnt;
int st[MAXN],tot;
int dis[MAXN],vis[MAXN],sum[MAXN];
int n,m,flag;

void intt(){
	clean(head,-1);
	ecnt=0;tot=0;
}
void add(int u,int v,int val){
	edge[ecnt]=node(v,val,head[u]);
	head[u]=ecnt++;
}
void dfs_spfa(int str){
	clean(dis,INF32);dis[str]=0;
	clean(vis,0);vis[str]=1;
	clean(sum,0);sum[str]++;
	st[tot++]=str;
	while(tot){
		int u=st[--tot];vis[u]=0;
		//printf("\n u: %d - ",u);
		if(sum[u]>n+1)	continue;
		for(int i=head[u];i+1;i=edge[i].nxt){
			int v=edge[i].v;
			if(dis[v]>dis[u]+edge[i].val){
				dis[v]=dis[u]+edge[i].val;
				//printf("v: %d %d",v,dis[v]);
				if(vis[v]==0){
					vis[v]=1;
					st[tot++]=v;
					sum[v]++;
					if(sum[v]>n){
						flag=1;return ;
					}
				}
			}
		}
	}
}
int main(){
	while(~scanf("%d",&n)){
		if(n==-1)	break;intt();
		int a,b,str=INF32,ed=-1;
		for(int i=1;i<=n;++i){
			scanf("%d%d",&a,&b);
			add(b+1,a,-2);
			str=min(str,a);ed=max(ed,b+1);
		}
		for(int i=str;i<ed;++i){
			add(i,i+1,1);
			add(i+1,i,0);
		}flag=0;
		dfs_spfa(ed);
//		printf("\n");
//		for(int i=str;i<=n+1;++i){ 
//			printf("%d ",dis[i]);
//		}printf("\n");
		printf("%d\n",abs(dis[str]));
	}
}
/*

*/