本来从昨天开始就想记录一下国庆这几天的 奥赛生活（结果因为懒。。
那就从今天开始好了（结果今天才刚开始
------------------------------------------------------8:23
需要从半夜时梦见帅气的学弟开始写起吗？（什么狗
需要从早上泡了一杯掺了奥利奥碎的原味咖啡写起吗？
需要从今天突然降温裹了一件外套来学校写起吗？
算了 老师来了 今天早上又是模拟考。。。
------------------------------------------------------8:31
离考试结束还有大概半个小时。。有点慌。。。
T1 图论题 EMM 先反向求两次最短路，然后再按拓扑序去DP，可是题目并没有保
证图中没有环啊（不知道会有多少分。。
T2 看起来像是树形DP，然而不会，按数据范围打了50分暴力（如果连50都没有。
。）
T3 好，完全懵逼。瞎推了一些神秘的性质，可是还是不知道怎么做啊，写了一个
根本不可能拿到分的骗分。。
最好的情况是拿到150分，也许连100分都没有。。
----------------------------------------------------11:35
老师还在给高一讲课，没空过来给我们评测 12:10
提交成功 12:13
原地爆炸 第一题才四十分，果然图是有环的！ 12:22
来到学校 开始消化 14:51
T1 只要三遍最短路就可以解决啊啊啊啊 ！
T1 完成-------15:13:06
T2 完成-------16:37:03
感觉T2的题解一点都不清晰 自己写一个好了
题目如下
【具有极强迷惑性又毫无用处的样例】
8 10 2 5 9 11 15 19 20 1 4 1 3 1 7 4 6 2 8 2 3 3 5
树形DP 那肯定要写出 DP方程啦 观察到n很小 所以
f[u][cu] 表示 假设此时距离u点最近的中心位于cu点，此时整个子树的费用之和（注意还没有加上在cu建立中心的K值）
然后每次更新的时候去找u 的所有直系儿子 设为 v 好了，可以用 f[v][cu] 直接更新，但如果离v最近的中心不是cu,而是是cv，那么就在v 的子树内部枚举所有的cv 找到一个最小的f[v][cv]，用f[v][cv]+K 来更新，那为什么此时又加上K了呢？因为cv在v 的子树内，u的中心已经不是cv而是cu了，你这个cv被建为中心已经是板上钉钉的事，必须加上K。
然后可以用一遍Floyd处理出任意两点间的距离，代码如下
1 #include<iostream> 2 #include<cstdio> 3 #include<cstdlib> 4 #include<cstring> 5 #include<algorithm> 6 7 #define For(i,a,b) for(register int i=a;i<=b;++i) 8 #define Re register 9 #define inf 0x7f7f7f 10 using namespace std; 11 const int N=210; 12 int f[N][N]; 13 int d[N][N],dx[N]; 14 int head[N],nxt[N*2],v[N*2],cnt=0,nf; 15 int n,m,x,y,z,kx; 16 inline void read(int &v){ 17 v=0; 18 char c=getchar(); 19 while(c<'0'||c>'9')c=getchar(); 20 while(c>='0'&&c<='9')v=v*10+c-'0',c=getchar(); 21 } 22 void add(int ux,int vx){ 23 cnt++; 24 nxt[cnt]=head[ux]; head[ux]=cnt; v[cnt]=vx; 25 cnt++; 26 nxt[cnt]=head[vx]; head[vx]=cnt; v[cnt]=ux; 27 } 28 29 int dfx[N],tot=0,bh[N],sz[N]; 30 void DFS(int x,int fa){ 31 dfx[++tot]=x; 32 sz[x]=1; bh[x]=tot; 33 for(Re int i=head[x];i;i=nxt[i])if(v[i]!=fa){ 34 DFS(v[i],x); 35 sz[x]+=sz[v[i]]; 36 } 37 } 38 39 void TrDP(int x,int fa){ 40 For(i,1,n)f[x][i]=dx[d[x][i]]; 41 for(Re int i=head[x];i;i=nxt[i])if(v[i]!=fa){ 42 int vv=v[i]; 43 TrDP(vv,x); 44 int mn=nf; 45 For(j,bh[vv],bh[vv]+sz[vv]-1){ 46 mn=min(mn,f[vv][dfx[j]]); 47 } 48 For(j,1,n)f[x][j]+=min(f[vv][j],mn+kx); 49 } 50 } 51 52 int main(){ 53 // freopen("ex.in","r",stdin); 54 freopen("logistics.in","r",stdin); 55 freopen("logistics.out","w",stdout); 56 read(n); read(kx); 57 For(i,1,n-1)read(dx[i]); 58 memset(d,inf,sizeof(d)); 59 nf=d[0][0]; 60 For(i,1,n-1){ 61 read(x); read(y); 62 add(x,y); 63 d[x][y]=d[y][x]=1; 64 } 65 66 For(i,1,n)d[i][i]=0; 67 For(k,1,n) For(i,1,n) For(j,1,n) 68 if(i!=j&&j!=k&&k!=i&&d[i][k]!=nf&&d[k][j]!=nf){ 69 d[i][j]=min(d[i][j],d[i][k]+d[k][j]); 70 } 71 72 73 DFS(1,0); 74 TrDP(1,0); 75 int fn=nf; 76 For(i,1,n)fn=min(fn,f[1][i]+kx); 77 cout<<fn<<endl; 78 fclose(stdin); fclose(stdout); 79 return 0; 80 }
17:13:25
18:33:30 回教室给桌上刚发芽的菊花浇了水，然后准备参加洛谷的模拟赛
21:55:23 还有35分钟结束，感觉还不错（？） 希望能拿到满分（？）
毕竟今天洛谷的运势可是大吉
22:34:45
AK了 不过有14个人AK呢 而且我是倒数第二的，不能太高兴，不能太高兴。。。。。
简略写一下思路吧（一定是因为 T2 开了没什么用的的三维DP数组 才慢的和狗一样）
T1 根据排列组合可以推出一个通式 Ans= n×2n-1-(2n-1) 就OK啦 不过草稿纸上的过程惨不忍睹
T2 DP题（开了三维的数组）f[i][j][0/1] 第 i 个Portal，时间到了j，第i个Portal 是否取，然后又开了一个Pa数组来记录是从哪里更新来的，最后从后往前去找出所有方案。
T3 本来以为NOIP不会考网络流的，想了二十几分钟的其它算法，后来才发现可以直接网络流建图（那些没学过网络流的NOIPer就有点可怜了）所以就把每个点拆成两个，建一个类似二分图的模型，流最大流判断是否有方案就OK了。最后还要在残流图上找到所有的方案，输出即可。
————————————23:00:59————————————
好了回家睡觉，今天就这样吧。
You can hear it in the silence~ silence~ you~~
You can feel it on the way home~ way home~ you~~
You can see it when the lights out~ light out~ you~~
——Taylor Swift "You Are in Love"