HDU 3666 THE MATRIX PROBLEM 差分约束

http://acm.hdu.edu.cn/showproblem.php?pid=3666

题目大意：

给你个N*M的矩阵，问是否存在一个序列a[1……N]和b[1……m]，使得矩阵中的每个元素L<=C[i][j] * a[i] /b[j]<=U

思路：

对于这种不等式有没有解的，可以想到差分约束。。。。。

L<=C[i][j] * a[i] /b[j]<=U 变形为

L/C[i][j] <=a[i] /b[j]<=U /C[i][j]  

两边取对数得：

log(L/C[I][J] )<=log(a[i])-log(b[j])<=log(UC[I][J] )

然后就是建图啦~

敲完代码后TLE。。。百度之，看是否负环判断人家是开了个根号。。。。

改了。。

无数次WA。。。。

检查老半天发现是dis数组我写int.....给我个豆腐好吗。。

改成double就AC了(queue(sqrt(n+m) 800MS)

后想起以前卡队列直接用栈过。。然后改了下（仍然是超过n+m）。1900+MS险过。。。。

而stack（sqrt(n+m))500MS。。。。。

#include <cstdio>#include<cstring>#include<stack>#include<cmath>#include <algorithm>using namespace std;const int MAXN=(400<<1)+10;const int MAXM=MAXN*MAXN+10;int n,m;int head[MAXN],len;struct edge{    int to,next;    double val;    }e[MAXM];void add(int from,int to,double val){    e[len].to=to;    e[len].val=val;    e[len].next=head[from];    head[from]=len++;}bool spfa(){    bool vis[MAXN];    int cnt[MAXN];    double dis[MAXN];    stack<int> q;    int tot=n+m;    int num=sqrt(double(n+m+1));    for(int i=0;i<tot;i++)    {        cnt[i]=1;        vis[i]=true;        dis[i]=0;        q.push(i);    }    while(!q.empty())    {        int cur=q.top();        q.pop();        vis[cur]=false;        for(int i=head[cur];i!=-1;i=e[i].next)        {            int id=e[i].to;            if(e[i].val+dis[cur] < dis[id])            {                dis[id]=e[i].val+dis[cur];                if(!vis[id])                {                    if(++cnt[id]>num)                        return false;                    q.push(id);                    vis[id]=true;                }            }        }    }    return true;}int main (){    double L,U;    while(~scanf("%d%d%lf%lf",&n,&m,&L,&U))    {        memset(head,-1,sizeof(head));        len=0;        for(int i=0;i<n;i++)        {            for(int j=0;j<m;j++)            {                int ci;                scanf("%d",&ci);                add(i,j+n,log(U/ci));                add(j+n,i,-log(L/ci));            }        }        if(spfa())            printf("YES\n");        else            printf("NO\n");    }    return 0;}