HDU-3666-THE MATRIX PROBLEM（差分约束系统）

最新推荐文章于 2020-11-23 18:21:09 发布

永远鲜红の幼月

最新推荐文章于 2020-11-23 18:21:09 发布

阅读量245

点赞数 1

分类专栏：程序设计竞赛记录

本文链接：https://blog.csdn.net/qq_40482358/article/details/87530501

版权

程序设计竞赛记录专栏收录该内容

339 篇文章 6 订阅

订阅专栏

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3666

题目大意：一个n*m的矩阵和两个数L，U，对于矩阵中的每个元素，存在两个数组an，bn，使得L<=arr[i][j]*a[i]/b[j]<=U；如果存在这样两个数组，输出YES，否则输出NO。

思路：这道题一开始我想的是暴力（想错了），就是找出：L<=arr[i][j]*a[i]<=U和L<=arr[i][j]/b[j]<=U，这个肯定是错的。。跟题目都不一样，亏我还找了半天。。

这道题的正解是差分约束系统，我们上面得到的式子，对所有数取对数就能得到差分约束系统的板子了，即：log(L/arr[i][j])<=log(a[i])-log(b[j])<=log(U/arr[i][j]).

然后就是直接跑板子了。

AC：

// luogu-judger-enable-o2
//#pragma comment(linker, "/STACK:1024000000,1024000000")
  
#include<stdio.h>
#include<string.h> 
#include<math.h> 
   
//#include<map>  
#include<set>
#include<deque> 
#include<queue> 
#include<stack> 
#include<bitset>
#include<string> 
#include<fstream>
#include<iostream> 
#include<algorithm> 
using namespace std; 
  
#define ll long long 
#define Pair pair<int,int>
//#define max(a,b) (a)>(b)?(a):(b)
//#define min(a,b) (a)<(b)?(a):(b)
#define clean(a,b) memset(a,b,sizeof(a))// 水印
//std::ios::sync_with_stdio(false);
//  register
const int MAXN=4e2+10;
const int INF32=0x3f3f3f3f;
const ll INF64=0x3f3f3f3f3f3f3f3f;
const int mod=998244353;
const double PI=acos(-1.0);

struct node{
	int v,nxt;double val;
	node(int _v=0,double _val=0,int _nxt=0){
		v=_v;val=_val;nxt=_nxt;
	}
}edge[MAXN*MAXN<<2];
int head[MAXN<<2],ecnt;
int stk[MAXN*MAXN<<2],tot;
double dis[MAXN<<2];
int vis[MAXN<<2],sum[MAXN<<2];
int n,m,L,U,flag;

void intt(){
	clean(head,-1);
	ecnt=0;tot=0;
}
void add(int u,int v,double val){
	edge[ecnt]=node(v,val,head[u]);
	head[u]=ecnt++;
}
void dfs_spfa(int str){
	int back=n+m;
	for(int i=1;i<=back;++i){
		dis[i]=999999999;vis[i]=0;sum[i]=0;
	}dis[str]=0;vis[str]=1;sum[str]++;
	stk[tot++]=str;
	while(tot){
		int u=stk[--tot];vis[u]=0;
		for(int i=head[u];i+1;i=edge[i].nxt){
			int v=edge[i].v;
			if(dis[v]>dis[u]+edge[i].val){
				dis[v]=dis[u]+edge[i].val;
				if(vis[v]==0){
					vis[v]=1;
					sum[v]++;
					stk[tot++]=v;
					if(sum[v]>back){
						flag=1;return ;
					}
				}
			}
		}
	}
}
int main(){
    while(~scanf("%d%d%d%d",&n,&m,&L,&U)){
    	intt();int val;
    	for(int i=1;i<=n;++i){
    		for(int j=1;j<=m;++j){
    			scanf("%d",&val);
    			add(n+j,i,log((U*1.0)/val));
    			add(i,n+j,-log((L*1.0)/val));
			}
		}flag=0;
		dfs_spfa(1);
		if(flag)	printf("NO\n");
		else	printf("YES\n");
	}
}


/*

*/