需要注意的是先缩点再连边,然后是连缩的点,tuopu的时候要注意写法,当入度减为0后才能入队!连边的时候为了降低复杂度,所以只要相邻的一对关系连起来就可以了,这样可以串到整个图!
#include<bits/stdc++.h>
using namespace std;
const int maxn=1e6+10;
int n,m;
int a[maxn],fa[maxn],b[maxn],f[maxn],dp[maxn],degree[maxn];
struct node{
int x,y;
}d[maxn];
vector<int>e[maxn];
bool cmp(node x,node y){
return x.x<y.x;
}
void init(){
for(int i=0;i<n*m;i++){
fa[i]=i;degree[i]=0;e[i].clear();
}
}
int find(int x){
if(fa[x]==x) return x;
return fa[x]=find(fa[x]);
}
void add(int x,int y){
e[x].push_back(y);degree[y]++;
}
int main(){
cin>>n>>m;
init();
//input
for(int i=0;i<n;i++){
for(int j=0;j<m;j++){
scanf("%d",&a[i*m+j]);
b[i*m+j]=a[i*m+j];
}
}
//suo dian
for(int i=0;i<n;i++){
for(int j=0;j<m;j++){
d[j].x=b[i*m+j];d[j].y=i*m+j;
}
sort(d,d+m,cmp);
for(int j=1;j<m;j++){
if(d[j].x==d[j-1].x){
fa[find(d[j].y)]=find(d[j-1].y);
}
}
}
for(int j=0;j<m;j++){
for(int i=0;i<n;i++){
d[i].x=b[i*m+j];d[i].y=i*m+j;
}
sort(d,d+n,cmp);
for(int j=1;j<n;j++){
if(d[j].x==d[j-1].x){
fa[find(d[j].y)]=find(d[j-1].y);
}
}
}
for(int i=0;i<n;i++){
for(int j=0;j<m;j++){
d[j].x=b[i*m+j];d[j].y=i*m+j;
}
sort(d,d+m,cmp);
for(int j=1;j<m;j++){
if(d[j].x!=d[j-1].x){
add(find(d[j-1].y),find(d[j].y));
}
}
}
for(int j=0;j<m;j++){
for(int i=0;i<n;i++){
d[i].x=b[i*m+j];d[i].y=i*m+j;
}
sort(d,d+n,cmp);
for(int j=1;j<n;j++){
if(d[j].x!=d[j-1].x){
add(find(d[j-1].y),find(d[j].y));
}
}
}
/*puts("********");
for(int i=0;i<n;i++){
for(int j=0;j<m;j++){
printf("%d ",(int)e[i*m+j].size());
}
puts("");
}
puts("*******");*/
//tuopu
queue<int>q;
for(int i=0;i<n*m;i++){
if((fa[i]==i)&&(degree[i]==0)) {
q.push(i);dp[i]=1;
// cout<<"i = "<<i<<endl;
}
}
while(!q.empty()){
int now=q.front();q.pop();
for(auto v : e[now]){
degree[v]--;
if(degree[v]==0) {
q.push(v);dp[v]=dp[now]+1;
}
}
}
for(int i=0;i<n;i++){
for(int j=0;j<m;j++){
printf("%d ",dp[find(i*m+j)]);
}
puts("");
}
}