bzoj2067: [Poi2004]SZN

我的未来我决定 提交于 2019-11-27 00:41:11

bzoj2067: [Poi2004]SZN

一开始没看出来是贪心,还以为是树规,多亏ooo提醒一句,然后刚了一个半小时搞出来了。

首先‘最长线最短’二分没错了,想了想他确实是单调的,最长线越长,用的线就越短(注意这里的最长线只是不超过,并不是一定要达到)。

二分最长线长度,对于已知的最长线长度len,考虑如何求解最少线数。

树里有什么特殊的点吗?叶子节点,叶子节点一定是一条线的一个端点。所以从叶子节点开始,一条线最优肯定是直接连接两个叶子节点,证明自己YY(其实是我不会)。

设f[i]表示节点i向上的线的长度。显然叶子节点的f为1。对于非叶子节点考虑如何合并:

f相加<=len的两个儿子可以合并,但并不是随便合并就是最优的。可以将儿子存到一个队列里,按f值从小到大排序,用两个指针实现合并。如果f[head]+f[tail]<=len直接合并head++,tail--,ans++。如果大于len直接把tail独立成一条线吗?这样并不是最优的,可以把tail的值记录加到f[x]到上一层合并。但是只能有一个点和x接上,显然选f最小的,其他的点就只能独立成一条线了。注意最后如果head==tail要特判是和x接上还是独立成一条线

这样直接合并到根就求出了最长线长度为len时的最少线数。

复杂度我不会算,有点玄但是跑的还挺快。

  1 //19.00~20.00 :93
  2 //20.00~20.25 :100
  3 #include<algorithm>
  4 #include<iostream>
  5 #include<cstring>
  6 #include<cstdio>
  7 #define int LL
  8 #define MAXN 20010
  9 #define LL long long
 10 #define ma(x) memset(x,0,sizeof(x))
 11 using namespace std;
 12 struct edge
 13 {    
 14     int u,v,nxt;
 15     #define u(x) ed[x].u
 16     #define v(x) ed[x].v
 17     #define n(x) ed[x].nxt
 18 }ed[MAXN*4];
 19 int first[MAXN],num_e;
 20 #define f(x) first[x]
 21 const int INF=0x7ffffffffff;
 22 int n,tans;
 23 int f[MAXN];
 24 
 25 pair<int,int> q[MAXN];int head,tail;
 26 void dfs(int x,int fa,int len)
 27 {
 28     int tem=0;
 29     for(int i=f(x);i;i=n(i))
 30     if(v(i)!=fa){dfs(v(i),x,len);}
 31     
 32     f[x]=0;if(x!=1)f[x]++;head=1,tail=0;int ooo=INF;
 33     for(int i=f(x);i;i=n(i))
 34     if(v(i)!=fa)
 35         if(f[v(i)]==len)tans++;
 36         else{q[++tail]=make_pair(f[v(i)],v(i));}
 37     sort(q+head,q+tail+1);
 38     while(head<tail)
 39     {
 40         if(q[head].first+q[tail].first>len)
 41         {
 42             if(q[tail].first<ooo)
 43             {
 44                 if(ooo!=INF)tans++;
 45                 ooo=q[tail].first;tail--;
 46             }
 47             else tans++,tail--;
 48         }
 49         else head++,tail--,tans++;
 50     }
 51     if(head==tail)
 52     {
 53         if(q[tail].first<ooo)
 54         {
 55             if(ooo!=INF)tans++;
 56             ooo=q[tail].first;
 57         }
 58         else tans++;
 59     }
 60     if(ooo!=INF)f[x]+=ooo;
 61 }
 62 int solve(int len)
 63 {
 64     ma(f);tans=0;    
 65     dfs(1,0,len);
 66     if(f[1])tans++;
 67     return tans;
 68 }    
 69 inline int read();
 70 inline void add(int u,int v);
 71 signed main()
 72 {    
 73 //    freopen("in.txt","r",stdin);    
 74     //freopen("1.out","w",stdout);
 75 
 76     n=read();int ta,tb;
 77     for(int i=1;i<n;i++)
 78         ta=read(),tb=read(),
 79         add(ta,tb),add(tb,ta);
 80     int l=1,r=n,mid,ans=solve(n),now=INF;
 81     int ans1,ans2;
 82     while(l<=r)
 83     {
 84         mid=(l+r)>>1;    
 85         int res=solve(mid);    
 86         if(res>ans)l=mid+1;
 87         else 
 88         {
 89             r=mid-1,ans=res;
 90             ans1=mid;ans2=res;
 91         }
 92     }
 93     printf("%lld %lld\n",ans2,ans1);
 94 }
 95 inline int read()
 96 {    
 97     int s=0,f=1;char a=getchar();
 98     while(a<'0'||a>'9'){if(a=='-')f=-1;a=getchar();}
 99     while(a>='0'&&a<='9'){s=s*10+a-'0';a=getchar();}
100     return s*f;
101 }
102 inline void add(int u,int v)
103 {    
104     ++num_e;
105     u(num_e)=u;
106     v(num_e)=v;
107     n(num_e)=f(u);
108     f(u)=num_e;
109 }
View Code

 

标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!