BZOJ1098: [POI2007]办公楼biu

1098: [POI2007]办公楼biu

Time Limit: 20 Sec  Memory Limit: 162 MB
Submit: 777  Solved: 326
[Submit][Status]

Description

FGD开办了一家电话公司。他雇用了N个职员，给了每个职员一部手机。每个职员的手机里都存储有一些同事的电话号码。由于FGD的公司规模不断扩大，旧的办公楼已经显得十分狭窄，FGD决定将公司迁至一些新的办公楼。 FGD希望职员被安置在尽量多的办公楼当中，这样对于每个职员来说都会有一个相对更好的工作环境。但是，为了联系方便起见，如果两个职员被安置在两个不同的办公楼之内，他们必须拥有彼此的电话号码。

Input

第一行包含两个整数N（2<=N<=100000）和M（1<=M<=2000000）。职员被依次编号为1,2,……,N. 以下M行，每行包含两个正数A和B(1<=A

Output

包含两行。第一行包含一个数S，表示FGD最多可以将职员安置进的办公楼数。第二行包含S个从小到大排列的数，每个数后面接一个空格，表示每个办公楼里安排的职员数。

Sample Input

7 16
1 3
1 4
1 5
2 3
3 4
4 5
4 7
4 6
5 6
6 7
2 4
2 7
2 5
3 5
3 7
1 7

Sample Output

3
1 2 4

HINT

 

FGD可以将职员4安排进一号办公楼，职员5和职员7安排进2号办公楼，其他人进3号办公楼。

题解：
我想到了答案就是反图的联通分量，可只会暴力的暴力。还以为复杂度很小，结果毫无悬念TLE+MLE+RE了
。
后来看了题解，觉得写个链表好麻烦，今天写感觉好简单。
巧妙的bfs+链表动态删点，所以每个点只入队一次，而且当前已经在队列里的就从链表里删去，复杂度？
不好估计，总之感觉很快。
还有一种想法是wyx528神牛在他的题解中提到的：
最直接的想法就是构造这个图的反图如去求连通块的个数及每个连通块的大小。但是构建反图需要 O(n2)
的时间，会超时。我们不妨这样想，n 个点 m 条边的无向图，记第 i 个连与 d[i]个点连通，则最小的
d [i]一定不会超过2m/n，这个很显然。我们记 d[i]最小的点为点 p，则与点 p 不相邻的 n－d[i]个点
一定都要跟点 p 放在一 起，这样我们就把这 n－d[i]个点合并了，对于剩下的 d[i]个点只要暴力求就
可以了。

但是实现起来比较不容易，我就没有写这个代码。

代码：

 1 #include<cstdio>
 2 #include<cstdlib>
 3 #include<cmath>
 4 #include<cstring>
 5 #include<algorithm>
 6 #include<iostream>
 7 #include<vector>
 8 #include<map>
 9 #include<set>
10 #include<queue>
11 #include<string>
12 #define inf 1000000000
13 #define maxn 150000
14 #define maxm 2100000
15 #define eps 1e-10
16 #define ll long long
17 #define pa pair<int,int>
18 #define for0(i,n) for(int i=0;i<=(n);i++)
19 #define for1(i,n) for(int i=1;i<=(n);i++)
20 #define for2(i,x,y) for(int i=(x);i<=(y);i++)
21 using namespace std;
22 inline int read()
23 {
24     int x=0,f=1;char ch=getchar();
25     while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
26     while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}
27     return x*f;
28 }
29 int n,m,tot=0,ans=0,a[maxn],head[maxn],v[maxn],fa[maxn],s[maxn],next[maxn],pre[maxn],q[maxn];
30 bool vis[maxn];
31 struct edge{int go,next;}e[2*maxm];
32 inline void insert(int x,int y)
33 {
34     e[++tot].go=y;e[tot].next=head[x];head[x]=tot;
35     e[++tot].go=x;e[tot].next=head[y];head[y]=tot;
36 }
37 inline int find(int x){return fa[x]==x?x:fa[x]=find(fa[x]);}
38 inline void del(int x)
39 {
40     next[pre[x]]=next[x];
41     pre[next[x]]=pre[x];
42 }
43 int main()
44 {   
45     freopen("input.txt","r",stdin);
46     freopen("output.txt","w",stdout);
47     n=read();m=read();int x,y;
48     for1(i,m)x=read(),y=read(),insert(x,y);
49     for1(i,n)fa[i]=i;
50     for0(i,n)pre[i]=i-1,next[i]=i+1;
51     while(1)
52     {
53     if(next[0]==n+1)break;    
54     int l=0,r=1,x,y;q[1]=next[0];del(q[1]);
55     while(l<r)
56     {
57         x=q[++l];
58         for(int i=head[x];i;i=e[i].next)v[y=e[i].go]=1;
59         for(int i=next[0];i!=n+1;i=next[i])if(!v[i])q[++r]=i,fa[i]=x,del(i);
60         for(int i=head[x];i;i=e[i].next)v[y=e[i].go]=0;
61     }
62     }
63     for1(i,n)s[find(i)]++;
64     for1(i,n)if(s[i])a[++ans]=s[i];
65     printf("%d\n",ans);
66     sort(a+1,a+ans+1);
67     for1(i,ans)printf("%d ",a[i]);
68     return 0;
69 }

View Code