hdu3666差分约束系统

题意：一个n*m的矩阵，两个数组a，b。矩阵中的每个元素L<=(Xij*Ai)/Bj<=U，问是否存在这样的两个数组。

这是一个差分约束的问题，将不等式变为Xij*Ai-U*Bj<=0   and   L*Bj-Xij*Ai<=0，但是由于是乘法形式，不符合差分约束Xi-Xj<=B的形式，故不等式两面同时取对数，即变为log(Ai)-log(Bj)<=log(U/Xij)   and   log(Bj)-log(Ai)<=log(Xij/L)， 然后按照差分约束的模板做就可以了。注意的是L，U是浮点型的，因为这个错误我找了好长时间。。。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
const int oo=1e8;
const int mn=161000;
const int mm=961000;
int edge,n,m;
double l,u;
int ver[mm],next[mm];
double dis[mn],cost[mm];
int head[mn],vis[mn],q[mn],outque[mn];
void addedge(int u,int v,double c)
{
    ver[edge]=v,cost[edge]=c,next[edge]=head[u],head[u]=edge++;
}
bool spfa()
{
    int i,u,v,l,r=0;
    double tmp;
    for(i=1;i<=n+m;i++) dis[i]=oo,vis[i]=1,outque[i]=0,q[r++]=i;
    dis[1]=0;
    for(l=0;l!=r;(++l>mn)?l=0:l)
        for(i=head[u=q[l]],vis[u]=0;i>=0;i=next[i])
            if(dis[v=ver[i]]>(tmp=dis[u]+cost[i]))
            {
                dis[v]=tmp;
               // cout<<"u="<<u<<" v="<<v<<endl;
                if(vis[v])  continue;
                if(++outque[v]>(int)sqrt((n+m)*1.0))    return false;
                vis[q[r++]=v]=1;
                //cout<<"l="<<l<<endl;
                if(r>=mn)   r=0;
            }
    return true;
}
int main()
{
    int x,i,j;
    while(~scanf("%d%d%lf%lf",&n,&m,&l,&u))
    {
        edge=0;
        memset(head,-1,sizeof(head));
        for(i=1;i<=n;i++)
            for(j=1;j<=m;j++)
            {
                scanf("%d",&x);
                addedge(j+n,i,log10(double(u/x)));
                addedge(i,j+n,log10(double(x/l)));
            }
        if(spfa()) puts("YES");
        else puts("NO");
    }
    return 0;
}