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 }
