事實證明,人類的腦殘是沒有極限的 他的意思就是簡單的要死的動規題 即使你已經搞出了動態轉移方程也會只拿70分
原題:Bessie正在計劃一年一度的奶牛大集會,來自全國各地的奶牛將來參加這一次集會。當然,她會選擇最方便的地點來舉辦這次集會。每個奶牛居住在 N(1<=N<=100,000) 個農場中的一個,這些農場由N-1條道路連接,並且從任意一個農場都能夠到達另外一個農場。道路i連接農場A_i和B_i(1 <= A_i <=N; 1 <= B_i <= N),長度爲L_i(1 <= L_i <= 1,000)。集會可以在N個農場中的任意一個舉行。另外,每個牛棚中居住者C_i(0 <= C_i <= 1,000)只奶牛。在選擇集會的地點的時候,Bessie希望最大化方便的程度(也就是最小化不方便程度)。比如選擇第X個農場作爲集會地點,它的不方便程度是其它牛棚中每隻奶牛去參加集會所走的路程之和,(比如,農場i到達農場X的距離是20,那麼總路程就是C_i*20)。幫助Bessie找出最方便的地點來舉行大集會。
首先 這一看就是一道簡單的樹形DP
對於百分之50的數據
我們可以枚舉根,進行轉移
這就是正常的50分做法。
對於100分算法,我們可以發現,其實真正的耗時的地方在於對於每一個點都要進行O(n)的驗證,顯然這是十分緩慢的,容易發現,從一個節點到其他的兒子進行跟的轉移,只需要將中間的邊進行處理。
所以兩次DFS就可以了
下面是代碼
#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <stdlib.h>
#include <time.h>
#define min(a,b) ((a)<(b)?(a):(b))
#define N 100000+5
#define M 200000+5
using namespace std;
long long ans=0x3f3f3f3f3fll;
int head[N],n,a[N],size[N];
long long f[N];
bool vis[N];
long long g[N];
struct graph
{
int next,to;
long long val;
graph (){}
graph (int _next,int _to,int _val)
:next(_next),to(_to),val(_val){}
}edge[M];
struct point
{
int x;
int val;
bool operator<(const point &z)const
{
return val<z.val;
}
}P[N];
inline void add(int x,int y,int z)
{
static int cnt = 0;
edge[++cnt]=graph(head[x],y,z);
head[x]=cnt;
edge[++cnt]=graph(head[y],x,z);
head[y]=cnt;
}
inline int read()
{
int x = 0, f = 1; char ch = getchar();
while (ch < '0' || ch > '9') { if (ch == '-') f = -1; ch = getchar(); }
while (ch >= '0' && ch <= '9') { x = x * 10 + ch - '0'; ch = getchar(); }
return x * f;
}
bool flag=false;
long long sum;
void DFS(int x,int fa)
{
size[x]=a[x];
f[x]=0;
for(int i=head[x];i;i=edge[i].next)
if(edge[i].to^fa)
{
DFS(edge[i].to,x);
size[x]+=size[edge[i].to];
f[x]+=(long long)size[edge[i].to]*(long long)edge[i].val+f[edge[i].to];
}
}
void reDFS(int x,int fa,int val)
{
if(x!=1)
g[x]=g[fa]+(sum-2*size[x])*(val);
else
g[x]=f[x];
for(int i=head[x];i;i=edge[i].next)
if(edge[i].to^fa)
reDFS(edge[i].to,x,edge[i].val);
}
int main()
{
freopen("meet.in", "r", stdin);
freopen("meet.out", "w", stdout);
cin>>n;
for(int i=1;i<=n;++i)
a[i]=read(),P[i].x=i,P[i].val=a[i];
for(int i=1;i<n;++i)
{
int tmp1=read(),tmp2=read(),tmp3=read();
add(tmp1,tmp2,tmp3);
}
DFS(1,-1);
ans=f[1];
sum=size[1];
reDFS(1,-1,0);
sort(g+1,g+n+1);
cout<<g[1];
}