bzoj 1006: [HNOI2008]神奇的国度 弦图最小染色 最大势算法

题意

给你一个弦图，问最少使用多少种颜色给图染色使得每条边连接的两个点均不同色。
n<=10000,m<=1000000

分析

详细的算法介绍可以去看cdq的论文。
具体来讲就是先用最大势算法把弦图的完美消除序列求出来，然后按完美消除序列倒序来染色，每个点染上能染的最小颜色，最后颜色数就是答案。
一个简单的证明：设颜色数为k，那么证明我们找到了一个大小为k的团，所以颜色数必然不小于k。但由于这里就有一个k的解，所以最小染色就是k。
最大势算法的具体实现cdq没有讲太清楚，这里有jcvb大大的讲解。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;

const int N=10005;

int n,m,last[N],col[N],cnt,lab[N];
bool cho[N],vis[N];
struct edge{int to,next;}e[N*200];
vector<int> vec[N];

int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

void addedge(int u,int v)
{
    e[++cnt].to=v;e[cnt].next=last[u];last[u]=cnt;
    e[++cnt].to=u;e[cnt].next=last[v];last[v]=cnt;
}

int main()
{
    n=read();m=read();
    for (int i=1;i<=m;i++)
    {
        int x=read(),y=read();
        addedge(x,y);
    }
    for (int i=1;i<=n;i++) vec[0].push_back(i);
    int mx=0;
    for (int i=n;i>=1;i--)
    {
        int x=0;
        while (1)
        {
            if (!vec[mx].size()) mx--;
            while (vec[mx].size()&&vis[vec[mx][vec[mx].size()-1]]) vec[mx].pop_back();
            if (vec[mx].size())
            {
                x=vec[mx][vec[mx].size()-1];
                break;
            }
        }
        vis[x]=1;
        for (int j=last[x];j;j=e[j].next)
            if (!vis[e[j].to]) lab[e[j].to]++,vec[lab[e[j].to]].push_back(e[j].to),mx=max(mx,lab[e[j].to]);
        for (int j=last[x];j;j=e[j].next)
            if (vis[e[j].to]) cho[col[e[j].to]]=1;
        for (int j=1;;j++)
            if (!cho[j])
            {
                col[x]=j;
                break;
            }
        for (int j=last[x];j;j=e[j].next)
            if (vis[e[j].to]) cho[col[e[j].to]]=0;
    }
    int ans=0;
    for (int i=1;i<=n;i++)
        if (!cho[col[i]]) cho[col[i]]=1,ans++;
    printf("%d",ans);
    return 0;
}