bzoj4446[Scoi2015]小凸玩密室

4446: [Scoi2015]小凸玩密室

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 419  Solved: 173
[Submit][Status][Discuss]

Description

小凸和小方相約玩密室逃脫，這個密室是一棵有n個節點的徹底二叉樹，每一個節點有一個燈泡。點亮全部燈

泡便可逃出密室。每一個燈泡有個權值Ai，每條邊也有個權值bi。點亮第1個燈泡不須要花費，以後每點亮4

個新的燈泡V的花費，等於上一個被點亮的燈泡U到這個點V的距離Du,v，乘以這個點的權值Av。在點燈

的過程當中，要保證任意時刻全部被點亮的燈泡必須連通，在點亮一個燈泡後必須先點亮其子樹全部燈泡才能點亮其餘燈泡。請告訴他們，逃出密室的最少花費是多少。

 

Input

第1行包含1個數n，表明節點的個數

第2行包含n個數，表明每一個節點的權值ai。(i=l，2，…，n)

第3行包含n-l個數，表明每條邊的權值bi，第i號邊是由第(i+1)/2號點連向第i+l號點的邊。

(i=l，2...N-1)

 

Output

輸出包含1個數，表明最少的花費。

 

Sample Input

3
5 1 2
2 1

Sample Output

5

 

HINT

 

對於100%的數據，1≤N≤2×105，1<Ai,Bi≤10^5

 

神題。
https://www.cnblogs.com/CXCXCXC/p/5312237.html

 1 #include<bits/stdc++.h>
 2 #define ll long long
 3 #define N 200010
 4 using namespace std;
 5 ll f[N][20],g[N][20];
 6 int n,v[N],w[N],dep[N],d[N];
 7 
 8 int main(){
 9     scanf("%d",&n);dep[1]=1;
10     for(int i=1;i<=n;i++)
11     scanf("%d",&v[i]);
12     for(int i=2;i<=n;i++){
13         scanf("%d",&w[i]);
14         dep[i]=dep[i>>1]+1;
15         d[i]=d[i>>1]+w[i];
16     }
17     for(int i=n;i;i--){
18         if((i<<1)>n){
19             for(int j=0;j<=dep[i];j++){
20                 int fa=i>>(dep[i]-j+1),to=(i>>(dep[i]-j))^1;
21                 f[i][j]=1ll*(d[i]+d[to]-2*d[fa])*v[to];
22             }
23         }
24         else if((i<<1)==n){
25             for(int j=0;j<=dep[i];j++){
26                 int to=i<<1;
27                 f[i][j]=f[to][j]+1ll*v[to]*w[to];
28             }
29         }
30         else{
31             int ls=i<<1,rs=ls|1;
32             for(int j=0;j<=dep[i];j++)
33             f[i][j]=min(1ll*v[ls]*w[ls]+f[ls][dep[ls]]+f[rs][j],1ll*v[rs]*w[rs]+f[rs][dep[rs]]+f[ls][j]);
34         }
35     }
36     for(int i=n;i;i--){
37         if((i<<1)>n){
38             for(int j=0;j<=dep[i];j++){
39                 int fa=i>>(dep[i]-j);
40                 g[i][j]=1ll*(d[i]-d[fa])*v[fa];
41             }
42         }
43         else if((i<<1)==n){
44             int s=i<<1;
45             for(int j=0;j<=dep[i];j++)
46             g[i][j]=1ll*v[s]*w[s]+g[s][j];
47         }
48         else{
49             int ls=i<<1,rs=ls|1;
50             for(int j=0;j<=dep[i];j++)
51             g[i][j]=min(1ll*v[ls]*w[ls]+f[ls][dep[ls]]+g[rs][j],1ll*v[rs]*w[rs]+f[rs][dep[rs]]+g[ls][j]);
52         }
53     }
54     ll ans=g[1][0];                     
55     for(int i=2;i<=n;i++){
56         ll ret=g[i][dep[i]-1];
57         for(int x=i;x>1;x>>=1){
58             int y=x^1;
59             if(y>n)ret+=1ll*v[x>>2]*w[x>>1];
60             else ret+=1ll*v[y]*w[y]+g[y][dep[y]-2];
61         }
62         ans=min(ans,ret);
63     }
64     printf("%lld\n",ans);
65     return 0;
66 }