码迷,mamicode.com
首页 > 其他好文 > 详细

[CF] 402 E. Strictly Positive Matrix

时间:2018-07-08 12:43:28      阅读:121      评论:0      收藏:0      [点我收藏+]

标签:head   含义   elm   比较   turn   接受   class   cout   pre   

一个矩阵,自乘无限次后能否全为正数?

如果n比较小,可以二分一下,但是这里n很大,乘一次都无法接受

可以考虑实际含义:矩阵看成邻接矩阵,那么0就是没有边,其余就是有边。

我们知道邻接矩阵自乘k次就相当于在原图走k步,无限次后是否有0?也就是图能否强连通。

 

判断就好,整个是环的情况题目限制不存在。

 

#include<iostream>
#include<cstdio>

using namespace std;

const int MAXN=2048;

struct Edge{
  int next,to;
}e[MAXN*MAXN*2];
int ecnt,head[MAXN];
inline void add(int x,int y){
  e[++ecnt].to = y;
  e[ecnt].next = head[x];
  head[x] = ecnt;
}

inline int rd(){
  int ret=0,f=1;char c;
  while(c=getchar(),!isdigit(c))f=c==-?-1:1;
  while(isdigit(c))ret=ret*10+c-0,c=getchar();
  return ret*f;
}

int n;

int bl[MAXN],cnt;
int ins[MAXN],sta[MAXN],top;
int dfn[MAXN],low[MAXN],tot;
void tarjan(int x){
  dfn[x]=low[x]=++tot;ins[x]=1;sta[++top]=x;
  for(int i=head[x];i;i=e[i].next){
    int v=e[i].to;
    // cout<<x<<" "<<v<<endl<<endl;
    if(!dfn[v]){
      tarjan(v);
      low[x]=min(low[x],low[v]);
    }else if(ins[v]){
      low[x]=min(low[x],dfn[v]);
    }
  }
  if(dfn[x]==low[x]){
    int elm;cnt++;
    do{
        elm = sta[top--];
        ins[elm] = false;
        bl[elm] = cnt;
    }while(elm != x);
  }
}

int main(){
  n=rd();int x;
  for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
      if(rd()) add(i,j);
  for(int i=1;i<=n;i++) if(!dfn[i]) tarjan(i);
  if(cnt!=1) return puts("NO"),0;
  puts("YES");
  return 0;
}

 

[CF] 402 E. Strictly Positive Matrix

标签:head   含义   elm   比较   turn   接受   class   cout   pre   

原文地址:https://www.cnblogs.com/ghostcai/p/9279609.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!