hdu 3666 THE MATRIX PROBLEM(判负环)

解题思路：

(U,L,mat[i][j])>0; U>=L;

因为题目要判定 L<=mat[i][j]*ai/bj<=U; 那么也就是判定L/mat[i][j]<=ai/bj<=U/mat[i][j] 即判定log(L/mat[i][j])<=log(ai)-log(bj)<=log(U/mat[i][j])的存在;

因为减号具有传递性（比如a-b=5,b-c=5由传递性得a-c=10），图也具有传递性（可以求出一颗最短路径树的图），还具有矛盾性（负环，下面会详细介绍），我们就可以把它抽象成图来用判定负环来判定矛盾，有矛盾即不存在；

log(ai)-log(bj)<=log(U/mat[i][j]) 根据最短路的三角不等式，可以抽象成bj到ai有一边且权值为log(U/mat[i][j])的边；

log(L/mat[i][j])<=log(ai)-log(bj) 也就是 log(bj)-log(ai)<=-log(L/mat[i][j])，即也可以抽象成ai到bj有一边且权值为-log(L/mat[i][j])的边，明显-log(L/mat[i][j])>=-log(U/mat[i][j])；

如果整个图有了负环，也就是说在负环上会出现，如果a到b的最短路为 t(t>0)，则b到a最短路至少会小于-t 这与我们题设是矛盾的(-log(L/mat[i][j])>=-log(U/mat[i][j]));

所以建图后只要判断是否有负环就可以了；

我这里用了两个版本的spfa算法，一个是dfs版的，一个是bfs版的，发现当最短路初始化为全0的时候dfs的效率要略高于bfs；

代码如下：

#include<stdio.h>
#include<math.h>
#include<string.h>
#include<queue>
using namespace std;

#define M 400000

int n,m,L,U;
int head[810],ant,to[M],next[M];
double cost[M];
void init()
{
    ant=0;
    memset(head,-1,sizeof(head));
}
void add(int a,int b,double c){
    to[ant]=b;
    cost[ant]=c;
    next[ant]=head[a];
    head[a]=ant++;
}

int vis[810];
double dp[810];
bool spfa()
{
    int i,u,times=0;
    queue<int> que;
    memset(vis,0,sizeof(vis));
    for(i=1;i<=n+m;i++){
         dp[i]=0;
         vis[i]=1;
         que.push(i);
    }
    while(!que.empty()){
        u=que.front();que.pop();
        for(i=head[u];i!=-1;i=next[i]){
            if(dp[u]+cost[i]<dp[to[i]]){
                times++;
                if(times>2*(n+m)) return true;
                dp[to[i]]=dp[u]+cost[i];
                if(!vis[to[i]]){
                    vis[to[i]]=1;
                    que.push(to[i]);
                }
            }
        }
        vis[u]=0;
    }
    return false;
}
int red[810];
bool dfs(int u)
{
    int i;
    if(vis[u]) return true;
    vis[u]=1;
    red[u]=1;
    for(i=head[u];i!=-1;i=next[i]){
        if(dp[u]+cost[i]<dp[to[i]]){
            dp[to[i]]=dp[u]+cost[i];
            if(dfs(to[i])) return true;
        }
    }
    vis[u]=0;
    return false;
}
bool dspfa()
{
    int i;
    memset(red,0,sizeof(red));
    memset(vis,0,sizeof(vis));
    for(i=1;i<=n+m;i++) dp[i]=0;
    for(i=1;i<=n+m;i++){
        if(!red[i]){
            if(dfs(i)) return true;
        }
    }
    return false;
}
int main()
{
    int i,j,k;
    while(scanf("%d%d%d%d",&n,&m,&L,&U)!=EOF)
    {
        init();
        for(i=1;i<=n;i++){
            for(j=1;j<=m;j++){
                scanf("%d",&k);
                add(j+n,i,log(U*1.0/k));
                add(i,j+n,-log(L*1.0/k));
            }
        }
        puts(dspfa()?"NO":"YES");
    }
    return 0;
}