博客
关于我
E. Strictly Positive Matrix [强联通+矩阵幂转图论+缩点] 好题
阅读量:537 次
发布时间:2019-03-08

本文共 2294 字,大约阅读时间需要 7 分钟。

E. Strictly Positive Matrix
time limit per test
1 second
memory limit per test
256 megabytes
input
standard input
output
standard output

You have matrix a of size n × n. Let's number the rows of the matrix from 1 to n from top to bottom, let's number the columns from 1 to nfrom left to right. Let's use aij to represent the element on the intersection of the i-th row and the j-th column.

Matrix a meets the following two conditions:

  • for any numbers i, j (1 ≤ i, j ≤ n) the following inequality holds: aij ≥ 0;
  • .

Matrix b is strictly positive, if for any numbers i, j (1 ≤ i, j ≤ n) the inequality bij > 0 holds. You task is to determine if there is such integer k ≥ 1, that matrix ak is strictly positive.

Input

The first line contains integer n (2 ≤ n ≤ 2000) — the number of rows and columns in matrix a.

The next n lines contain the description of the rows of matrix a. The i-th line contains n non-negative integers ai1, ai2, ..., ain(0 ≤ aij ≤ 50). It is guaranteed that .

Output

If there is a positive integer k ≥ 1, such that matrix ak is strictly positive, print "YES" (without the quotes). Otherwise, print "NO" (without the quotes).

Examples
input
Copy
21 00 1
output
NO
input
Copy
54 5 6 1 21 2 3 4 56 4 1 2 41 1 1 1 14 4 4 4 4
output
YES

题意:给定一个矩阵A,其中A[i][j]>=0 && sigma[i][i]>0 问,是否存在K,使得B=A^K,其中B[i][j]>0

思路:矩阵乘法转图论. A^k大于0,代表所有元素直接是联通.

B=A^K 中 b[i][j]代表i到j长度为k的路条数

#include
#define bug cout <<"bug"<
edge[MAX_N];int n,cur,cnt,top,mp[MAX_N],DFN[MAX_N],LOW[MAX_N],s[MAX_N];void tarjan(int u){ DFN[u]=LOW[u]=cur++; s[++top]=u; for(int i=0;i<(int)edge[u].size();i++){ int v=edge[u][i]; if(!DFN[v]){ tarjan(v); LOW[u]=min(LOW[u],LOW[v]); } else if(DFN[v] && !mp[v]) LOW[u]=min(LOW[u],DFN[v]); } if(DFN[u]==LOW[u]){ int v=-1; cnt++; while(u!=v){ v=s[top--]; mp[v]=cnt;// cout << v<<" "; }// puts("************"); }}int main(void){ read(n); for(int i=1;i<=n;i++){ for(int j=1;j<=n;j++){ int x;read(x); if(i==j) continue; if(x>0) edge[i].push_back(j); } } for(int i=1;i<=n;i++) if(!DFN[i]) tarjan(i);// cout << cnt << endl; if(cnt==1) cout <<"YES"<

转载地址:http://jkkiz.baihongyu.com/

你可能感兴趣的文章
MySQL_西安11月销售昨日未上架的产品_20161212
查看>>
Mysql——深入浅出InnoDB底层原理
查看>>
MySQL“被动”性能优化汇总
查看>>
MySQL、HBase 和 Elasticsearch:特点与区别详解
查看>>
MySQL、Redis高频面试题汇总
查看>>
MYSQL、SQL Server、Oracle数据库排序空值null问题及其解决办法
查看>>
mysql一个字段为空时使用另一个字段排序
查看>>
MySQL一个表A中多个字段关联了表B的ID,如何关联查询?
查看>>
MYSQL一直显示正在启动
查看>>
MySQL一站到底!华为首发MySQL进阶宝典,基础+优化+源码+架构+实战五飞
查看>>
MySQL万字总结!超详细!
查看>>
Mysql下载以及安装(新手入门,超详细)
查看>>
MySQL不会性能调优?看看这份清华架构师编写的MySQL性能优化手册吧
查看>>
MySQL不同字符集及排序规则详解:业务场景下的最佳选
查看>>
Mysql不同官方版本对比
查看>>
MySQL与Informix数据库中的同义表创建:深入解析与比较
查看>>
mysql与mem_细说 MySQL 之 MEM_ROOT
查看>>
MySQL与Oracle的数据迁移注意事项,另附转换工具链接
查看>>
mysql丢失更新问题
查看>>
MySQL两千万数据优化&迁移
查看>>