题目
多组询问，给定一个(n*m)的矩阵(C)和一个区间([L,R])，
问是否存在一个长度为(n)的序列(A)和一个长度为(m)的序列(B)，
使得所有
[frac{A_i}{B_j}C_{i,j}in [L,R]
]
不需要输出具体方案
分析
这其实是一个不等式，变量是(A_i)和(B_j)，这样可以化简成
[frac{L}{C_{i,j}}leq frac{A_i}{B_j}leq frac{R}{C_{i,j}}
]
这样不够好做，考虑把它化成对数的形式，即是
[log(L)-log(C_{i,j})leq log(A_i)-log(B_j)leq log(R)-log(C_{i,j})
]
按照这个建图判断有没有负环即可，有负环就是无解
但是bfs版的spfa会TLE，考虑SLF优化，只跑了700ms（但是这样更容易被卡掉）
代码
#include <cstdio>
#include <cctype>
#include <cmath>
#include <deque>
#define rr register
using namespace std;
const int N=811; deque<int>q;
int v[N],cnt[N],as[N],S,n,m,et; double dis[N],l,r;
struct node{int y; double w; int next;}e[N*N>>1];
inline signed iut(){
rr int ans=0,f=1; rr char c=getchar();
while (!isdigit(c)) f=(c=='-')?-f:f,c=getchar();
while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();
return ans*f;
}
inline void add(int x,int y,double w){
e[++et]=(node){y,w,as[x]},as[x]=et;
}
inline bool spfa(){
while (!q.empty()) q.pop_front();
for (rr int i=1;i<=S;++i)
dis[i]=1e18,v[i]=cnt[i]=0;
q.push_back(S); dis[S]=0,v[S]=1;
while (!q.empty()){
rr int x=q.front(); q.pop_front();
for (rr int i=as[x];i;i=e[i].next)
if (dis[e[i].y]>dis[x]+e[i].w){
dis[e[i].y]=dis[x]+e[i].w;
if (++cnt[e[i].y]>=S) return 0;
if (!v[e[i].y]){
v[e[i].y]=1;
if (!q.empty()&&dis[e[i].y]<dis[q.front()]) q.push_front(e[i].y);
else q.push_back(e[i].y);
}
}
v[x]=0;
}
return 1;
}
signed main(){
while (scanf("%d",&n)==1){
m=iut(),S=n+m+1,et=0,
l=log(iut()),r=log(iut());
for (rr int i=1;i<=S;++i) as[i]=0;
for (rr int i=1;i<S;++i) add(S,i,0);
for (rr int i=1;i<=n;++i)
for (rr int j=1;j<=m;++j){
rr double x=log(iut());
add(i,j+n,r-x),add(j+n,i,x-l);
}
puts(spfa()?"YES":"NO");
}
return 0;
}