【BZOJ 1791】 [Ioi2008]Island 岛屿

Description

你将要游览一个有N个岛屿的公园。从每一个岛i出发，只建造一座桥。桥的长度以Li表示。公园内总共有N座桥。尽管每座桥由一个岛连到另一个岛，但每座桥均可以双向行走。同时，每一对这样的岛屿，都有一艘专用的往来两岛之间的渡船。 相对于乘船而言，你更喜欢步行。你希望所经过的桥的总长度尽可能的长，但受到以下的限制。 • 可以自行挑选一个岛开始游览。 • 任何一个岛都不能游览一次以上。 • 无论任何时间你都可以由你现在所在的岛S去另一个你从未到过的岛D。由S到D可以有以下方法： o 步行：仅当两个岛之间有一座桥时才有可能。对于这种情况，桥的长度会累加到你步行的总距离；或者 o 渡船：你可以选择这种方法，仅当没有任何桥和/或以前使用过的渡船的组合可以由S走到D（当检查是否可到达时，你应该考虑所有的路径，包括经过你曾游览过的那些岛）。 注意，你不必游览所有的岛，也可能无法走完所有的桥。 任务 编写一个程序，给定N座桥以及它们的长度，按照上述的规则，计算你可以走过的桥的最大长度。 限制 2 <= N <= 1,000,000 公园内的岛屿数目。 1<= Li <= 100,000,000 桥i的长度。

Input

• 第一行包含N个整数，即公园内岛屿的数目。岛屿由1到N编号。 • 随后的N行每一行用来表示一个岛。第i 行由两个以单空格分隔的整数，表示由岛i筑的桥。第一个整数表示桥另一端的岛，第二个整数表示该桥的长度Li。你可以假设对於每座桥，其端点总是位于不同的岛上。

Output

你的程序必须向标准输出写出包含一个整数的单一行，即可能的最大步行距离。 注1：对某些测试，答案可能无法放进32-bit整数，你要取得这道题的满分，可能需要用Pascal的int64或C/C++的long long类型。 注2：在比赛环境运行Pascal程序，由标准输入读入64-bit数据比32-bit数据要慢得多，即使被读取的数据可以32-bit表示。我们建议把输入数据读入到32-bit数据类型。 评分 N不会超过4,000。

Sample Input

7
3 8
7 2
4 2
1 4
1 9
3 4
2 3

Sample Output

24

HINT

此题为寻找基环树上的最长链

又是一道基环树DP，此处处理基环的方式拆开复制一下，然后寻找环上的最长链，此处注意2的时候要特判

关于树上的最长路就直接跑DP就好了

我知道我写的很不详细，但是代码很清楚啊

orz lyd

 #include<cstdio>
 #include<iostream>
 using namespace std;
 const int N=;
 struct ee{int to,next,w;}e[N*];
 long long head[N],c[N],f[N],du[N],d[N],b[N*],a[N*],q[*N];
 int n,timer,cnt;
 long long ans;
 bool vis[N];
 void ins(int u,int v,int w){
     e[++cnt].to=v;e[cnt].next=head[u];e[cnt].w=w;head[u]=cnt;du[u]++;
     e[++cnt].to=u;e[cnt].next=head[v];e[cnt].w=w;head[v]=cnt;du[v]++;
 }

 void dfs(int now,int k){
     c[now]=k;
     for (int i=head[now];i;i=e[i].next){
         int v=e[i].to;
         if(!c[v]) dfs(v,k);
     }
 }

 void topsort(){
     int l=,r=;
     for (int i=;i<=n;i++) if(du[i]==) q[++r]=i;
     while(l<r) {
         int now=q[++l];
         for (int i=head[now];i;i=e[i].next){
             int v=e[i].to;
             if(du[v]>){
                 du[v]--;
                 d[c[now]]=max(d[c[now]],f[now]+f[v]+e[i].w);
                 f[v]=max(f[v],f[now]+e[i].w);
                 if(du[v]==)q[++r]=v;
             }
         }
     }
 }

 void dp(int t,int x){
     int m=,y=x,i;
     do{
         a[++m]=f[y];du[y]=;
         for(i=head[y];i;i=e[i].next){
             int v=e[i].to;
             if(du[v]>){
                 b[m+]=b[m]+e[i].w;
                 y=e[i].to;
                 break;
             }
         }
     }while(i);
     if(m==){//
         int l=;
         for (int i=head[y];i;i=e[i].next)
             if(e[i].to==x) l=max(l,e[i].w);
         d[t]=max(d[t],f[x]+f[y]+l);
         return;
     }
     for(int i=head[y];i;i=e[i].next){
         int v=e[i].to;
         if(v==x) {
             b[m+]=b[m]+e[i].w;
             break;
         }
     }
     for (int i=;i<=m;i++){
             a[i+m]=a[i];
             b[m+i]=b[m+]+b[i];
         }
     int l,r;
     q[l=r=]=;
     for (int i=;i<*m;i++){
         while(l<=r&&i-q[l]>=m)l++;
         d[t]=max(b[i]-b[q[l]]+a[i]+a[q[l]],d[t]);
         while(l<=r&&a[q[r]]+b[i]-b[q[r]]<=a[i]) r--;
         q[++r]=i;
     }

 }

 int main(){
     scanf("%d",&n);
     int v,w;
     for (int i=;i<=n;i++){
         scanf("%d%d",&v,&w);
         ins(i,v,w);
     }
     for (int i=;i<=n;i++) if (!c[i]) dfs(i,++timer);
     topsort();
     for (int i=;i<=n;i++){
         if(du[i]>&&!vis[c[i]]) {
             vis[c[i]]=;
             dp(c[i],i);
             ans+=d[c[i]];
         }
     }
     cout<<ans<<endl;
 }