题目链接

大意

给定一个\(N\times N\)的矩阵,问是否满足以下三个条件:

  1. \(A_{i,i}=0\)
  2. \(A_{i,j}=A_{j,i}\)
  3. 对于任意的\(i,j,k\),满足\(A_{i,j}\ge max(A_{i,k},A_{k,j})\)

如果满足条件输出MAGIC,否则输出NOT MAGIC

思路

首先,对于第一二两个条件,直接判就行了。

考虑第三个条件:

对于这个矩形,有个直观的想法就是建成一个完全图。

点\(i\)和点\(j\)间的边权就是\(A_{i,j}\),那么第三个条件就转化为了:

对于任意一个三元环的三条边\(A,B,C\)



就有以下式子:

\[\begin{cases}
A\ge max(B,C)\\
B\ge max(A,C)\\
C\ge max(A,B)
\end{cases}
\]

化简一下即为,\(A,B,C\)中的最大值等于\(A,B,C\)中的次大值。

考虑如何快速地判断每一个三元环:

我们从小到大加边,权值相同的边一起加入,那么在任意时刻,图像中的每个连通块一定都是一个完全图,否则就不满足条件。

  • 证明: 若图中的某个连通块不是一个完全图,那么一定会缺少一些边,形成如下图的例子:



    而又由于1,3之间的边一定会严格大于\(A,B\),则该三元环的最大值不等于次大值,即不满足上述条件。

对于判断加边后是否每个连通块都是完全图,我们可以用最小生成树的思想来做:

我们知道加完所有边后的图像是一个完全图,那么我们在加边时只用判断一下该边连接的两个点是否在同一连通块内就行了。如果在同一连通块内,就说明之前形成了一个缺边的非完全图,就直接不合法了。

注:权值相同的边一定要一起加入呀。

代码

其实一个就是Kruskal算法。

#include<map>
#include<cstdio>
#include<vector>
#include<algorithm>
using namespace std;
const int MAXN=2505;
int N,A[MAXN][MAXN],Len;
int Fa[MAXN],B[MAXN*MAXN],Cnt,NCnt;
struct Edge{int x,y,z;}s[MAXN*MAXN];
vector<int>P[MAXN*MAXN];map<int,int>Mp;
bool cmp(Edge A,Edge B){return A.z<B.z;}
int Find(int x){return Fa[x]==x?x:Fa[x]=Find(Fa[x]);}
int main(){
scanf("%d",&N);
for(int i=1;i<=N;i++)
for(int j=1;j<=N;j++)
scanf("%d",&A[i][j]);
for(int i=1;i<=N;i++){
if(A[i][i]){
printf("NOT MAGIC\n");
return 0;
}
for(int j=i+1;j<=N;j++){
if(A[i][j]!=A[j][i]){
printf("NOT MAGIC\n");
return 0;
}Len++;
s[Len].x=i;
s[Len].y=j;
s[Len].z=A[i][j];
B[++Cnt]=A[i][j];
}
}
sort(B+1,B+Cnt+1);
sort(s+1,s+Len+1,cmp);B[0]=-1;
for(int i=1;i<=Cnt;i++)
if(B[i]!=B[i-1])Mp[B[i]]=++NCnt;
for(int i=1;i<=Len;i++){
s[i].z=Mp[s[i].z];
P[s[i].z].push_back(i);
}
for(int i=1;i<=N;i++)Fa[i]=i;
for(int i=1;i<=NCnt;i++){
int size=P[i].size();
for(int j=0;j<size;j++){
int u=P[i][j];
int x=Find(s[u].x);
int y=Find(s[u].y);
if(x==y){
printf("NOT MAGIC\n");
return 0;
}
}
for(int j=0;j<size;j++){
int x=Find(s[P[i][j]].x);
int y=Find(s[P[i][j]].y);
if(x!=y)Fa[x]=y;
}
}
printf("MAGIC\n");
}

【CF632F】Magic Matrix(生成树 脑洞)的更多相关文章

  1. 【思维题 经典模型】cf632F. Magic Matrix

    非常妙的经典模型转化啊…… You're given a matrix A of size n × n. Let's call the matrix with nonnegative elements ...

  2. CF 1042 E. Vasya and Magic Matrix

    E. Vasya and Magic Matrix http://codeforces.com/contest/1042/problem/E 题意: 一个n*m的矩阵,每个位置有一个元素,给定一个起点 ...

  3. Educational Codeforces Round 9 F. Magic Matrix 最小生成树

    F. Magic Matrix 题目连接: http://www.codeforces.com/contest/632/problem/F Description You're given a mat ...

  4. Codeforces 632F Magic Matrix(bitset)

    题目链接  Magic Matrix 考虑第三个条件,如果不符合的话说明$a[i][k] < a[i][j]$ 或 $a[j][k] < a[i][j]$ 于是我们把所有的$(a[i][j ...

  5. codeforces 632F. Magic Matrix (最小生成树)

    You're given a matrix A of size n × n. Let's call the matrix with nonnegative elements magic if it i ...

  6. codeforces 632F. Magic Matrix

    题目链接 给一个n*n的矩阵, 问是否对角线上的元素全都为0, a[i][j]是否等于a[j][i], a[i][j]是否小于等于max(a[i][k], a[j][k]), k为任意值. 前两个都好 ...

  7. Codeforces 632F - Magic Matrix(暴力 bitset or Prim 求最小生成树+最小瓶颈路)

    题面传送门 开始挖老祖宗(ycx)留下来的东西.jpg 本来想水一道紫题作为 AC 的第 500 道紫题的,结果发现点开了道神题. 首先先讲一个我想出来的暴力做法.条件一和条件二直接扫一遍判断掉.先将 ...

  8. CF1042E Vasya and Magic Matrix

    感觉不会期望. 首先把所有格子按照权值从小到大排一下序,这样一共有$n * m$个元素,每个元素有三个属性$x, y, val$. 下文中的下标均为排序后的下标. 这样子我们就可以推出公式: $f_i ...

  9. Vasya and Magic Matrix CodeForces - 1042E (概率dp)

    大意:给定n*m矩阵, 初始位置(r,c), 每一步随机移动到权值小于当前点的位置, 得分为移动距离的平方, 求得分期望. 直接暴力dp的话复杂度是O(n^4), 把距离平方拆开化简一下, 可以O(n ...

随机推荐

  1. Win10 开启 Hyper-V 及简单使用

    简介 Windows 10 上内置了 Hyper-V.Hyper-V 提供硬件虚拟化,每个虚拟机都在虚拟硬件上运行. 系统要求 Windows 10 企业版.专业版或教育版.家庭版.移动版.移动企业版 ...

  2. Pytest_Hook钩子函数总结(14)

    前言 pytest 的钩子函数有很多,通过钩子函数的学习可以了解到pytest在执行用例的每个阶段做什么事情,也方便后续对pytest二次开发学习.详细文档可以查看pytest官方文档https:// ...

  3. concat模糊查询

    <if test="name!=null"> name like concat('%',concat(#{name},'%')) </if> choose ...

  4. c# - 按引用内存地址传参 和 按输出传参 的具体使用

    1.前言 传递参数,不需要返回值,对懒人很舒服哟,缺点是不好定位数据 2.操作 using System; namespace ConsoleApp1.letVlaueGo { public clas ...

  5. PowerShell 管道符之Select的使用方法【二】

    这次讲解Select中的第二个方法:String 在我们的ISE编辑器中输入如下命令 Select-String - 可以了解到,原来这是正则表达式,它提供了一些额外的正则方法.具体如何使用,可以自行 ...

  6. CentOS7防火墙firewalld 和 CentOS6防火墙iptables的一些配置命令

    CentOS7 防火墙 一.防火墙的开启.关闭.禁用.查看状态命令 (1)启动防火墙:systemctl start firewalld (2)关闭防火墙:systemctl stop firewal ...

  7. Spring Boot Admin,贼好使!

    Spring Boot Admin(SBA)是一个开源的社区项目,用于管理和监控 Spring Boot 应用程序.应用程序可以通过 http 的方式,或 Spring Cloud 服务发现机制注册到 ...

  8. winfrom 双缓冲

    在窗体load函数中 this.DoubleBuffered = true; //控件,需要反射的方式设置 Type dgvType = this.dgv.GetType(); PropertyInf ...

  9. [USB波形分析] 全速USB波形数据分析(一)

    在之前的文章一次CAN波形分析之旅里,根据示波器采集的波形数据,详细地分析了CAN通信.今天来分析USB数据,还是同样的流程,但是这次使用matplotlib来协助分析. USB基本波形 USB通过一 ...

  10. VS2017:win32项目与win32控制台应用程序的转换方法

    原文:https://www.cnblogs.com/asuser/articles/12297251.html 刚开始使用VS2017新建项目工程时,有时把应用类型的工程建成控制台类型的工程,在编译 ...