[CF632F]Magic Matrix

题目

传送门 to luogu

思路

注意 $a_{ij}=a_{ji}$ 这个条件很舒服，让我们有思路将 $a_{ij}$ 试做图上 $i$ 和 $j$ 的无向边权值。抛开矩阵背景吧！

显然 $x\le y$ 时 $\max (n,x)\le \max (n,y)$ 。一个必要条件浮出水面：

$$a_{il}\le \max (a_{ik},a_{kl})\le \max [\max (a_{ij},a_{jk}),a_{kl}]=\max (a_{ij},a_{jk},a_{kl})$$

可以看出，这是有传递性的。于是，若 $i,j$ 之间存在一条路径，只由边权小于 $a_{ij}$ 的边构成，则 $\mathtt{N}\mathtt{O}\mathtt{T}\mathtt{M}\mathtt{A}\mathtt{G}\mathtt{I}\mathtt{C}$ 。

继续，我们可以证明这是充分的。题目中的条件本就是 $a_{ij}\le \max (a_{ik},a_{kj})$ ，难道这不是一种路径？

所以，我们获得了充要条件：权值小于 $a_{ij}$ 的边不能使 $i,j$ 联通。这不是最小生成树的过程吗？

于是我们获得了 $\mathcal{O}(n^2\log n)$ 的算法，解决了此题。瓶颈在于排序。

---

$2021/3/1\mathtt{u}\mathtt{p}\mathtt{d}\mathtt{a}\mathtt{t}\mathtt{e}$：为什么不用 $\mathrm{p}\mathrm{r}\mathrm{i}\mathrm{m}$ 呢？我当时是脑子被门夹了还是什么？

如果一个点的 $dis$ 被更新了两次，那就是 $\mathtt{N}\mathtt{O}\mathtt{T}\mathtt{M}\mathtt{A}\mathtt{G}\mathtt{I}\mathtt{C}$ 的，多简单！

代码

#include <cstdio>
#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;
inline int readint(){
	int a = 0; char c = getchar(), f = 1;
	for(; c<'0'||c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c&&c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}
template < typename T >
void getMax(T&a,const T&b){if(a<b)a=b;}
template < typename T >
void getMin(T&a,const T&b){if(b<a)a=b;}

const int MaxN = 2501;
namespace ufs{
	int fa[MaxN];
	void init(int n){
		for(int i=1; i<=n; ++i)
			fa[i] = i;
	}
	int find(int a){
		if(a != fa[a])
			fa[a] = find(fa[a]);
		return fa[a];
	}
	bool linked(int a,int b){
		return find(a) == find(b);
	}
	void link(int a,int b){
		fa[find(a)] = find(b);
	}
}

struct Node{
	int x, y, val;
	bool operator < (const Node &that) const {
		return val < that.val;
	}
} node[MaxN*MaxN];

int n;
void kruskal(){
	sort(node+1,node+n*n+1);
	node[n*n+1].val = -1;
	ufs::init(n);
	for(int l=1,r=0; l<=n*n; l=r+1){
		while(node[r+1].val == node[l].val)
			++ r;
		if(node[l].val != 0)
		for(int i=l; i<=r; ++i)
			if(ufs::linked(node[i].x,node[i].y)){
				puts("NOT MAGIC"); return ;
			}
		for(int i=l; i<=r; ++i)
			ufs::link(node[i].x,node[i].y);
	}
	puts("MAGIC");
}

int main(){
	n = readint();
	for(int i=1; i<=n; ++i)
		for(int j=1; j<=n; ++j){
			Node &now = node[(i-1)*n+j];
			now.x = i, now.y = j;
			now.val = readint();
			if(i == j && now.val){
				puts("NOT MAGIC");
				return 0;
			}
		}
	kruskal();
	return 0;
}

后记

更有趣的性质：考虑 $a_{ij},a_{jk},a_{ki}$ 三条边。不妨设 $a_{ij}$ 最大，因为本质是三元环上的最大边。那么我们有 $a_{ij}\le \max (a_{jk},a_{ki})$ ，可是 $a_{ij}\ge \max (a_{jk},a_{ki})$ ，故而 a i j ∈ { a j k , a k i } a_{ij}\in\{a_{jk},a_{ki}\} aij​∈{ajk​,aki​} 。

可以看到，我们的每个三元环都是一个“等腰三角形”，且腰比底更长（至少不更短）。这是充要条件。

不过我没有以此想出算法。或许你可以试一试！