hdu3666差分约束系统

来源：互联网发布：淘宝评论设置方法编辑：程序博客网时间：2024/05/16 14:50

题意：一个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;}