51Nod 算法马拉松21(迎新年)
这次打算法马拉松是在星期五的晚上，发挥还算正常(废话，剩下的题都不会= =)。
讲讲比赛经过吧。
8:00准时发题，拿到之后第一时间开始读。
A配对，看上去像是二分图最大权匹配，一看范围吓傻了，先跳过读后面的题。
B完全二叉树的方差，大概看了一遍，好神的样子，跳过。
C多项式？好吧没学过FFT和NTT的我肯定不会，跳跳跳。
D最大值，哎呦这函数什么破玩意儿，看不懂，跳跳跳。
E B君的射击，卧槽毕克大人您出题就算了出这么一道码农题是要闹那样，跳跳跳。
F那些年，我们一起讲的故事，卧槽这特么简直就是道纯语文题，跳跳跳。
哎等等，跳没了……
好像A比较简单，先想想A算了……
第一思路是二分图最大权匹配，然后费用流暴力，一看数据范围这显然不靠谱，然后就开始思考问题的特殊性质。
不经意间按了几下F5，卧槽怎么你们做的那么快，看来是我想复杂了，这应该是个不太难的题。
想了想匹配的构造方法，贪心似乎不太可行，别的方法一时也想不出来，算了好像可以枚举每条边计算贡献……感觉像是对的，但是又对此深表怀疑……
中间又按了几次F5，卧槽你们做题怎么这么神速……看来这题确实简单，反正WA了也不会有惩罚，看我码一发枚举边计算贡献，随手一交，啊哈居然A了……
自己交的比较早，所以是第18个A的(为什么我记得是第19个……还是说是第19个提交的……)，排名还不错……(机房其他人好多都第几十+才A的……)
看了看其他题，E题貌似之前看少了连通的圆可以使中间那块也变成洞，感觉自己并不会写。F题读了半天还是没读懂……D题的函数看懂了就是十进制位倒过来然后乱搞，然而还是没思路。C题也看了两眼，一看别人的运行时间都是十几ms，感觉像是O(1)的结论题，打了一发print input()(这是Python2)然而过不去样例……
然后开始搞B题，看见方差最小感觉像是一个三分平均数然后二分图最小权匹配，搞了半天搞出来新树直径的结论，然而并没有推出题解的神贪心，写了一发三分平均数+最小费用最大流(二分图最小权匹配)，然后不知道哪儿写挂了，样例都过不去(也可能是函数根本不单峰……)……调了半天调不出来，狠狠心一交，1W19T，算了我弃疗……
脑子昏昏涨涨，就此弃疗。星期六和星期天也没有再来做题，由于第一题交的早，最终成绩rank60，真是便宜我了。(然而排名太低，Rating往下掉了……= =)
随便水水比赛就拿到了排名是20的倍数这个buff，话说运气真是好……
话说今天才收到点头盾，这速度我给差评……
给A题贴个题解：
既然要求匹配点的距离和最大，那么两条路径一定要尽可能相交(如果不相交，交换两条路径的端点不会变劣)。或者说，匹配路径一定要尽可能经过长的边。不过这还是不太好搞，所以放弃构造方法，考虑对每条边计算能出现在几条匹配路径中，也就是说能对答案产生多少贡献。
考虑一条长为w的边，设这条边连接的两个连通块的大小分别为a和b。为了保证距离和最大，我们需要使得尽量多的匹配路径经过这条边，所以应该尽量让边两端的点进行匹配，因此这条边最多经过min{a,b}次。
边与边之间互不影响(我不会证……)，所以一遍dfs/bfs之后枚举每条边计算贡献即可。为了保险(防爆栈)，我用的是bfs。
1 #include<cstdio> 2 #include<cstring> 3 #include<algorithm> 4 #include<vector> 5 using namespace std; 6 const int maxn=100010; 7 struct edge{int u,v,w;}e[maxn]; 8 void bfs(); 9 vector<int>G[maxn]; 10 int n,q[maxn],prt[maxn]={0},size[maxn]={0}; 11 long long ans=0; 12 int main(){ 13 scanf("%d",&n); 14 for(int i=1;i<n;i++){ 15 scanf("%d%d%d",&e[i].u,&e[i].v,&e[i].w); 16 G[e[i].u].push_back(e[i].v); 17 G[e[i].v].push_back(e[i].u); 18 } 19 bfs(); 20 for(int i=1;i<n;i++){ 21 if(prt[e[i].u]==e[i].v)swap(e[i].u,e[i].v); 22 int s=min(size[e[i].v],size[1]-size[e[i].v]); 23 ans+=(long long)s*e[i].w; 24 } 25 printf("%lld",ans); 26 return 0; 27 } 28 void bfs(){ 29 int head=0,tail=0; 30 q[tail++]=1; 31 while(head!=tail){ 32 int x=q[head++]; 33 size[x]=1; 34 for(int i=0;i<(int)G[x].size();i++)if(G[x][i]!=prt[x]){ 35 prt[G[x][i]]=x; 36 q[tail++]=G[x][i]; 37 } 38 } 39 for(int i=n;i;i--){ 40 int x=q[i]; 41 size[prt[x]]+=size[x]; 42 } 43 }
附官方题解：
(话说构造法真的可行嘛……我仍然理解不了贪心构造……只会写O(n2m2)的费用流……囧囧囧)
反思：
这次打比赛的时候有点着急，B题也没有耐心仔细想/写，发挥并不太好。
A题看见别人神速过题之后就有点慌了，一改求稳的作风，凭着直觉交了一份代码，能A估计是运气。今后的比赛或者是考试一定要求稳，毕竟无缘无故掉了几十分上百分可是很要命的。
下一站，UOJ Test Round #2。
加油。
