hdu 3666 THE MATRIX PROBLEM

题意：给你一个n*m的矩阵，和一个区间上下界L，U，问你是否能找出n个数a1,a2,a3,......,an和m个数b1,b2,b3,......,bm，使得矩阵中的第i行的数乘以ai得到的值范围在[L,U]内，而且矩阵中第j列的数除以bj得到的值范围也在[L,U]内，能找到你就输出YES，否则NO。

根据题意我们可以简单得到下面的关系：

（1） L<=Aij * ai <=U （第一个要求）

（2） L<=Aij / bj <=U （第二个要求）

看到上面两个式子，似乎和差分约束的不等式关系不符，但是我们可以通过取对数使其变成两数相加或者两数相减，通过取对数可以得到下面变形式：

（1） log(L) <= log(Aij) +log(ai) <=log(U)

（2） log(L) <= log(Aij) - log(bj) <=log(U)

哈哈哈，两式相加可以得到： 2*log(L) <= 2*log(Aij) +log(ai) -log(bj) <=2*log(U)

                            移项可得： 2*log(L)-2*log(Aij) <=log(ai) - log(bj) <= 2*log(U)-2*log(Aij)

那么就可以在i到j+n之间建边了，边长就是上式中的值，具体看下代码就是了，但是此题卡队列，用一般的方法入队次数大于点数判环会TLE，在这里可以用以下方法优化：

 如果点数为n，那么其中一个点的入队次数大于sqrt(n) 就有环

#include<iostream>
#include<cstring>
#include<cmath>
#include<queue>
#include<vector>
#include<cstdio>
using namespace std;
#define inf 0x7f7f7f7f
int head[50000], n, num, m, cnt[50000];bool vis[50000];double dis[50000];
struct node
{
	int en, next;double len;
}E[500000];
void init()
{
	num = 0;
	memset(head, -1, sizeof(head));
}
void add(int st, int en ,double len)
{
	E[num].en = en;E[num].len = len;
	E[num].next = head[st];head[st] = num++;
}
bool spfa()
{
	int i;
	for (i = 0;i <= n+m;i++)
	{
		dis[i] = inf;vis[i] = false;
		cnt[i] = 0;
	}
	queue<int>q;
	dis[1] = 0.0;vis[1] = true;cnt[1] ++;
	q.push(1);
	while (!q.empty())
	{
		int x = q.front();
		q.pop();
		vis[x] =false;
		for (i = head[x];i != -1;i = E[i].next)
		{
			int v= E[i].en;
			if (dis[v] >dis[x] + E[i].len)
			{
				dis[v] = dis[x] + E[i].len;
				if (!vis[v])
				{
					cnt[v]++;q.push(v);
					if (cnt[v] >(int)sqrt(n+m)) return false;
					vis[v] = true;
				}
			}
		}
	}
	return true;
}
int main()
{
	double l, u, i, j,x;
	while (scanf("%d%d%lf%lf", &n, &m, &l, &u) != EOF)
	{
		init();
		u = 2.0*log(u), l = 2.0*log(l);
		for (i = 0;i < n;i++)
		{
			for (j = 0;j < m;j++)
			{
				scanf("%lf",&x);x = 2.0 * log(x);
				add(j , i+m, u - x);
				add(i + m, j, x - l);
			}
		}
		bool flag = spfa();
		if (flag == 0) printf("NO\n");
		else printf("YES\n");
	}
	return 0;
}