Codeforces Round #345 (Div. 1) C. Table Compression

掉回div2了悲伤啊= =

题意是给出一个n*m的矩阵，要求出一个新的矩阵，同一行（列）中的任意两数大小关系不变（大，小，等于三种关系），同时使得新矩阵的最大元素最小，输出任意方案。

大体思路就是把矩阵转化成有向图，然后在图中找最长路。

把矩阵的每个元素都对应一个节点，对每一行（列）排序之后对相邻的两个元素的对应节点连一条边。这里有个问题是相邻两个元素可能相等，这个时候可以用并查集来缩点。然后就把矩阵转化成有向图了。

至于这个最长路嘛。。。刚开始脑残地按照最短路bfs写得，然后就写煞笔了= =赛后想想水dp啊，真是脑残无药医。主要还是刚开始惯性思维有向边是从小的元素指导大一看就像是最短路，后来把边的方向反过来就是明显的dp了。

$$dp[x]=max(dp[x],dp[child]+1)$$
注意这里的dp实现可以和上面的缩点结合到一起，搜到这个节点时直接把同一层（元素相等）的结点的儿子都搜完了，然后把最大值记录起来作为这一层的答案。这样就不用在建图时候专门费劲的缩点了.

/* ***********************************************
Author        :axp
Created Time  :2016/3/7 18:25:23
TASK          :C.cpp
LANG          :C++
************************************************ */

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <cmath>
#include <cstdlib>
#include <ctime>
using namespace std;

const int N = 1e6+10;
int n,m;
int ar[N];
int pos[N];
int cnt;
vector<int> e[N];
vector<int> se[N];
int ans[N];
queue<int> q;
int dp[N];
int p[N];

int find(int x)
{
    return p[x]==x?x:p[x]=find(p[x]);
}

bool cmp(int x,int y)
{
    return ar[x]<ar[y];
}

void work(int x,int y)
{
    //cout<<x<<' '<<y<<endl;
    if(ar[x]<ar[y])swap(x,y);
    if(ar[x]==ar[y])
    {
        p[find(x)]=find(y);
        return;
    }
    //cout<<x<<' '<<y<<endl;
    e[x].push_back(y);
}

int dfs(int x)
{
    x=find(x);
    if(dp[x])return dp[x];
    dp[x]=1;
    for(int i=0;i<se[x].size();i++)
        for(int j=0;j<e[se[x][i]].size();j++)
        {
            int to=e[se[x][i]][j];
            dp[x]=max(dp[x],dfs(to)+1);
        }
    return dp[x];
}

int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            cnt++;
            scanf("%d",&ar[cnt]);
        }
    for(int i=1;i<=cnt;i++)p[i]=i;

    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
            pos[j]=(i-1)*m+j;
        sort(pos+1,pos+1+m,cmp);
        for(int j=1;j<m;j++)
            work(pos[j],pos[j+1]);
    }
    for(int j=1;j<=m;j++)
    {
        for(int i=1;i<=n;i++)
            pos[i]=(i-1)*m+j;
        sort(pos+1,pos+1+n,cmp);
        for(int i=1;i<n;i++)
            work(pos[i],pos[i+1]);
    }

    for(int i=1;i<=cnt;i++)
    {
        int fa=find(i);
        se[fa].push_back(i);
    }

    for(int i=1;i<=cnt;i++)
        ans[i]=dfs(find(i));
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            printf("%d%c",ans[(i-1)*m+j],j==m?'\n':' ');
    return 0;
}