hdu 3666 THE MATRIX PROBLEM

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

差分约束，关键在于找不等式，然后建图

这个题 spfa 判断负环时，如果以更新次数大于等于N判断有负环，会超时

只要判断更新次数大于sqrt(N)就可以，原理不知道

代码：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<string>
#include<set>
#include<map>
#include<cmath>
#include<algorithm>
#include<vector>
#include<cmath>
#include<queue>
#include<stack>

//#define ull unsigned long long

using namespace std;


const int INF=0x3f3f3f3f;
const int MOD=1000000007;
const double eps=1e-6;
const int N=10005;
const int M=1000005;
int head[N],I;
struct node
{
    int j,next;
    double d;
}edge[M];
double dist[N];
queue<int>qt;
bool had[N];
int in[N];
void add(int i,int j,double d)
{
    edge[I].j=j;
    edge[I].d=d;
    edge[I].next=head[i];
    head[i]=I++;
}
bool spfa(int n)
{
    while(!qt.empty())
    qt.pop();

    for(int i=0;i<n;++i)
    {
        dist[i]=0.0;
        had[i]=true;
        in[i]=1;
        qt.push(i);
    }
    while(!qt.empty())
    {
        int x=qt.front();
        qt.pop();
        had[x]=false;
        for(int t=head[x];t!=-1;t=edge[t].next)
        {
            int j=edge[t].j;
            if(dist[j]>dist[x]+edge[t].d)
            {
                dist[j]=dist[x]+edge[t].d;
                if(!had[j])
                {
                    had[j]=true;
                    ++in[j];
                    qt.push(j);
                    if(in[j]>sqrt(1.0*n))
                    return false;
                }
            }
        }
    }
    return true;
}
void init(int n,int m,double L,double U)
{
    memset(head,-1,sizeof(head));I=0;
    for(int i=0;i<n;++i)
    for(int j=0;j<m;++j)
    {
        double ftmp;
        scanf("%lf",&ftmp);
        add(i,j+n,log(ftmp/L));
        add(j+n,i,log(U/ftmp));
    }
}
int main()
{
    //freopen("data.in","r",stdin);
    int n,m;
    double L,U;
    while(scanf("%d %d %lf %lf",&n,&m,&L,&U)!=EOF)
    {
        init(n,m,L,U);

        if(spfa(n+m))
        printf("YES\n");
        else
        printf("NO\n");
    }
    return 0;
}