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;
}