hdu 3666 THE MATRIX PROBLEM

题目意思就是是否存在ai,bj，使得l<=cij*(ai/bj)<=u (1<=i<=n,1<=j<=m)成立

首先，把cij除到两边:l'<=ai/bj<=u'，如果差分约束的话，应该是ai-bj的形式，于是可以取对数

log(l')<=log(ai)-log(bj)<=log(u')

 

把log(ai)和log(bj)看成两个点ai和bj，化成求最短路的形式：dis[ai]-dis[bj]<=log(u'),dis[bj]-dis[ai]<=-log(l')

然后判负环就行，深搜和广搜都可以

 

注意的是，如果spfa队列判负环：

(1)不必判断某个点入队次数大于N,只要判断是否大于sqrt(1.0*N)

(2)或者所有点的入队次数大于T*N，即存在负环，一般T取2

N为所有点的个数

 

代码：

spfa队列版本

#include<iostream> #include<memory.h> #include<string> #include<cstdio> #include<algorithm> #include<math.h> #include<stack> #include<queue> using namespace std; const int MAX=805; struct node { int v,next; double c; }g[MAX*MAX]; int adj[MAX]; int n,m,e; double dis[MAX],l,u; bool vis[MAX]; int cnt[MAX]; void add(int u,int v,double c) { g[e].v=v; g[e].c=c; g[e].next=adj[u]; adj[u]=e++; } bool spfa() { int i,u,v; int limit=(int)sqrt(1.0*(n+m)); //int limit=2*(n+m); queue<int>que; memset(vis,0,sizeof(vis)); memset(cnt,0,sizeof(cnt)); for(i=0;i<=n+m;i++) { dis[i]=0; que.push(i); } while(!que.empty()) { u=que.front(); que.pop(); vis[u]=false; for(i=adj[u];i!=-1;i=g[i].next) { v=g[i].v; if(dis[v]>dis[u]+g[i].c) { dis[v]=dis[u]+g[i].c; if(!vis[v]) { vis[v]=true; if(++cnt[v]>limit)//或者cnt[v]>sqrt(1.0*(n+m)) return false; que.push(v); } } } } return true; } int main() { int i,j; double t; while(scanf("%d%d %lf %lf",&n,&m,&l,&u)!=EOF) { e=0; memset(adj,-1,sizeof(adj)); for(i=1;i<=n;i++) { for(j=1;j<=m;j++) { scanf("%lf",&t); add(j+n,i,log(u/t)); add(i,j+n,-log(l/t)); } } if(spfa()) puts("YES"); else puts("NO"); } return 0; } 

 

spfa深搜版本

#include<iostream> #include<memory.h> #include<string> #include<cstdio> #include<algorithm> #include<math.h> #include<stack> #include<queue> using namespace std; const int MAX=805; struct node { int v,next; double c; }g[MAX*MAX]; int adj[MAX]; int n,m,e; double dis[MAX],l,u; bool vis[MAX],inStack[MAX]; inline void add(int u,int v,double c) { g[e].v=v; g[e].c=c; g[e].next=adj[u]; adj[u]=e++; } bool spfa(int u) { int i,v; if(inStack[u]) return false; inStack[u]=true; vis[u]=true; for(i=adj[u];i!=-1;i=g[i].next) { v=g[i].v; if(dis[v]>dis[u]+g[i].c) { dis[v]=dis[u]+g[i].c; if(!spfa(v)) { return false; } } } inStack[u]=false; return true; } bool ok() { int i,u,v,cnt=0; memset(vis,0,sizeof(vis)); memset(inStack,0,sizeof(inStack)); //memset(cnt,0,sizeof(cnt)); for(i=0;i<=n+m;i++) { dis[i]=0; //que.push(i); } for(i=1;i<=n+m;i++) { if(!vis[i]) { if(!spfa(i)) { return false; } } } return true; } int main() { int i,j; double t; while(scanf("%d%d %lf %lf",&n,&m,&l,&u)!=EOF) { e=0; memset(adj,-1,sizeof(adj)); for(i=1;i<=n;i++) { for(j=1;j<=m;j++) { scanf("%lf",&t); add(j+n,i,log(u/t)); add(i,j+n,-log(l/t)); } } if(ok()) puts("YES"); else puts("NO"); } return 0; }