HDU3666_差分约束

HDU3666——差分约束

[题目链接]http://acm.hdu.edu.cn/showproblem.php?pid=3666

  总结：本题是差分约束的典型题目之一。差分约束主要是将a-b<=w一类的不等式约束问题转换为图论中求解最短（最长）路径的方法，实际上属于图论类型的题目。对于这类题，关键点是从题目中找出合适的约束不等式来建图。
　　本题中，对于矩阵中的第i行j列元素Cij，我们都可以列出不等式L<=Cij*Ai/Bj<=U，移项后我们可以得到关于Ai和Bj的两个不等式，由于我们需要的是a-b<=w一类的不等式，所以我们可以考虑到用对数将除法转换为减法，这样就找到了需要的约束不等式。
　　差分约束中，若我们要求a-b的最大值，则要求解的是从b到a的最短路。而求解最短路时我们可能会得到三种结果：
　　1.求出了最短路，长度为w。表示a-b的最大值为w。
　　2.b到a不连通。表示从约束条件中找不到约束关系，即a-b的最大值有无穷多解。
　　3.图中b到a之间存在负值圈。此时最短路w可以无穷小，即a-b最大值无解。
　　本题我们要判断是否存在序列{An}和{Bm}满足约束条件，也就是判断图中是否有负值圈，用spfa算法能很容易完成。
　　代码如下：

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <string>
#include <vector>
#include <algorithm>
#include <queue>
#include <math.h>
#include <map>
#include <stdlib.h>
#include <ctype.h>
using namespace std;
#define ll long long
const int maxn=1000,INF=1e9+7;
int N,M,L,U;
double lgL,lgU;
struct Edge{
	int from,to;double cost;
	Edge(int a,int b,double c){
		from=a,to=b,cost=c;
	}
};
vector<Edge> edges;
vector<int> G[maxn];
void AddEdge(int from,int to,double cost){
	edges.push_back(Edge(from,to,cost));
	G[from].push_back(edges.size()-1);
}
void build(void){
	edges.clear();
	for(int i=0;i<=N+M;++i) G[i].clear();
	for(int i=1;i<=N+M;++i) AddEdge(0,i,0);
	for(int i=1;i<=N;++i)
		for(int j=1;j<=M;++j){
			int x;scanf("%d",&x);
			double lgc=log(1.0*x);
			AddEdge(N+j,i,lgU-lgc);
			AddEdge(i,N+j,lgc-lgL);
		}
}
int inq[maxn],cnt[maxn];
double d[maxn];
bool spfa(void){						
	memset(inq,0,sizeof(inq));
	memset(cnt,0,sizeof(cnt));
	for(int i=0;i<=N+M;++i) d[i]=INF;
	deque<int> q;					// SLF优化
	q.push_back(0);d[0]=0;
	inq[0]=1;
	while(!q.empty()){
		int u=q.front();q.pop_front();
		inq[u]=0;
		for(int i=0;i<G[u].size();++i){
			Edge& e=edges[G[u][i]];
			if(d[e.to]>e.cost+d[u]){
				d[e.to]=e.cost+d[u];
				if(!inq[e.to]){
					inq[e.to]=1;
					if(!q.empty()&&d[e.to]<d[q.front()])
						q.push_front(e.to);
					else q.push_back(e.to);
					if(++cnt[e.to]>N+M) return false;
				}
			}
		}
	}
	return true;
}
int main(void){
	while(scanf("%d%d%d%d",&N,&M,&L,&U)!=EOF){
		lgL=log(L),lgU=log(U);
		build();
		printf("%s\n",spfa()?"YES":"NO");
	}
	return 0;
}