【NOIP2016提高A组集训第16场11.15】兔子

题目大意

在一片草原上有N个兔子窝，每个窝里住着一只兔子，有M条路径连接这些窝。更特殊地是，至多只有一个兔子窝有3条或更多的路径与它相连，其它的兔子窝只有1条或2条路径与其相连。换句话讲，这些兔子窝之前的路径构成一张N个点、M条边的无向连通图，而度数大于2的点至多有1个。
兔子们决定把其中K个兔子窝扩建成临时避难所。当危险来临时，每只兔子均会同时前往距离它最近的避难所躲避，路程中花费的时间在数值上等于经过的路径条数。为了在最短的时间内让所有兔子脱离危险，请你安排一种建造避难所的方式，使最后一只到达避难所的兔子所花费的时间尽量少。

n,k≤1000 m≤1500

分析

根据题目，只有一个点度数大于2。画一下图发现这张图是一堆环和一堆链，比如：

其中度数大的点是那个中心点。

这道题显然满足二分性质，那么二分答案，先确定一个len，然后贪心地构造最小的答案。

首先搞链的情况。
对于一个洞，它可以覆盖掉链上2len+1个点，那么就一直这样覆盖，直到剩余部分不足2len+1。然后分类讨论：
假设剩余部分长度为t，如果t>l，那么如果在外面放洞，是不能覆盖完多余部分，所以应该继续放一个进来，那么会多出来一点长度。记所有这种情况的多出长度的最大值为Mx。
如果t≤l，那么如果在外面放洞，是存在使剩余部分被覆盖的方案的。所以先不在里面放，记这堆东西的最大值为Mi

那么现在比较Mx和Mi的大小。如果Mx≥Mi，表示我第一种情况多出的部分可以把第二种情况缺少的部分填满，并且多出的部分还可以给每个环覆盖掉2Mx+1的长度。
如果Mx < Mi，表示我多出的部分不够填，所以应该在中心点上再放一个洞，把缺的部分填满。（为什么一定要在中心放洞？如果在环上其它填满，对这个环的贡献依然是2len+1，而对其它环的贡献不足2len+1，如果在中心点填满，对每个环的贡献都是2len+1）
这样就处理掉了链的答案，并且得出对每个环的贡献，那么把环割掉贡献的长度后，剩余部分又是一条链，用相同的方法贪心去填满即可。

时间复杂度$O(nlogn)$

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>

using namespace std;

const int N=1005,M=3005;

typedef long long LL;

int n,m,K,l,r,mid,h[N],e[M],nxt[M],top,tot,Dep[N],D[N],len[N],ltot,R,o[M];

char c;

bool bz[N];

int read()
{
    int x=0,sig=1;
    for (c=getchar();c<'0' || c>'9';c=getchar()) if (c=='-') sig=-1;
    for (;c>='0' && c<='9';c=getchar()) x=x*10+c-48;
    return x*sig;
}

void add(int x,int y,int c)
{
    e[++tot]=y; nxt[tot]=h[x]; h[x]=tot; D[x]++; o[tot]=tot+c;
}

void init(int x,int F)
{
    top++;
    bz[x]=1;
    bool loop=0,b=0;
    for (int i=h[x];i;i=nxt[i])
    {
        if (i==F) continue;
        if (!bz[e[i]])
        {
            init(e[i],o[i]); b=1;
        }else if (e[i]==R)
        {
            len[ltot++]=top; loop=1; break;
        }
    }
    if (!loop && !b) Dep[tot++]=top;
    top--;
}

bool Check(int l)
{
    int cnt=0,Mx=0,Mi=0,R;
    for (int i=0;i<tot;i++)
    {
        int t=Dep[i]%(2*l+1);
        cnt+=Dep[i]/(2*l+1);
        if (t>l)
        {
                        Mx=max(Mx,2*l-t); cnt++;
                }
        else Mi=max(Mi,t);
    }
    if (Mx>=Mi) R=Mx*2+1;else
    {
        R=l*2+1; cnt++;
    }
    for (int i=0;i<ltot;i++)
    if (R<len[i])
    {
        int t=len[i]-R;
        cnt+=(t+2*l)/(2*l+1);
    }
    return cnt<=K;
}

int main()
{
    freopen("rabbit.in","r",stdin); freopen("rabbit.out","w",stdout);
    n=read(); m=read(); K=read();
    while (m--)
    {
        int x=read(),y=read();
        add(x,y,1); add(y,x,-1);
    }
    for (R=1;D[R]<=2;R++);
    tot=0;
    init(R,-1);
    for (int i=0;i<tot;i++) Dep[i]--;
    for (l=0,r=n-1,mid=r>>1;l<r;mid=l+r>>1)
        if (Check(mid)) r=mid;else l=mid+1;
    printf("%d\n",l);
        fclose(stdin); fclose(stdout);
        return 0;
}