bzoj3242(基环树DP)

最新推荐文章于 2024-07-30 13:58:21 发布

qkoqhh

最新推荐文章于 2024-07-30 13:58:21 发布

阅读量342

点赞数

分类专栏： DP

本文链接：https://blog.csdn.net/qkoqhh/article/details/81267140

版权

DP 专栏收录该内容

116 篇文章 1 订阅

订阅专栏

先考虑树的情况。。显然是直径的一半（然而窝并没有看出来qaq

然后环上的情况。。。首先处理一下外向树的情况。。跑一遍dfs可以解决。。顺便处理最大深度。。。

然后是最远路径经过环的情况。。逐点对求实在是非常麻烦，考虑的情况很多。。

这需要特殊的姿势（就是题解）。。

考虑已经确定了答案点，那么这个点延伸出去能得到一棵最短路路径树。。显然这棵最短路径树必将不会经过环上的一条边，因此，可以枚举环上的边，依次删除求直径。。这样O(n^2)显然不行。。

我们将外向树缩成点后可以当成序列环DP来做。。（然后这里就暴露出了窝做题太少。。

考虑环上断了之后，直径有2种情况。。一种是不跨过序列，即在序列断边后的两侧，一种是跨过序列，即在左右都有半段。。

所以要对前后缀算出完整的直径和一半的直径，枚举断边分别取最大值即可。。

要注意到断边之后可能会出现新叶子节点。。所以直接把环上节点也当叶子节点看待就可以合并这个特殊情况。。

/**
 *　　　　　　  　　┏┓　　 　┏┓
 * 　　　　　  　　┏┛┗━━━━━━━┛┗━━━┓
 * 　　　　　  　　┃　　　　　　　┃ 　
 * 　　　　　  　　┃　　　━　　 　┃
 * 　　　　　  　　┃　＞　　　＜　┃
 * 　　　　　  　　┃　　　　　　　┃
 * 　　　　　  　　┃...　⌒　... 　┃
 * 　　　　　　  　┃              ┃
 * 　　　　　　  　┗━┓          ┏━┛
 * 　　　　　　　　　┃          ┃　Code is far away from bug with the animal protecting　　　　　　　　　　
 * 　　　　　　　　　┃          ┃   神兽保佑,代码无bug
 * 　　　　　　　　　┃          ┃　　　　　　　　　　　
 * 　　　　　　　　　┃          ┃  　　　　　　
 * 　　　　　　　　　┃          ┃
 * 　　　　　　　　　┃          ┃　　　　　　　　　　　
 * 　　　　　　　　　┃          ┗━━━┓
 * 　　　　　　　　　┃              ┣┓
 * 　　　　　　　　　┃              ┏┛
 * 　　　　　　　　　┗┓┓┏━━━━━━━━┳┓┏┛
 * 　　　　　　　　　　┃┫┫       ┃┫┫
 * 　　　　　　　　　　┗┻┛       ┗┻┛
 */
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
#include<map>
#include<stack>
#include<set>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define eps 1e-8
#define succ(x) (1LL<<x)
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define mid (x+y>>1)
#define NM 100005
#define nm 400000
#define N 1000005
#define M(x,y) x=max(x,y)
const double pi=acos(-1);
const ll inf=1e18;
using namespace std;
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}
 
 
 
struct edge{int t;ll v;edge*next;}e[nm],*h[NM],*o=e;
void add(int x,int y,ll v){o->t=y;o->v=v;o->next=h[x];h[x]=o++;}
int n,c[NM],f[NM],tot,cnt,low[NM];
ll d[NM],a[NM],ans,_t,_x,_y,g[NM],s,pre[NM],nxt[NM],b[NM];
bool v[NM];
 
void tar(int x){
    low[x]=++cnt;
    link(x)if(f[x]!=j->t)if(!low[j->t])f[j->t]=x,tar(j->t);
    else if(low[j->t]>low[x])for(int y=j->t;y!=f[x];y=f[y])c[++tot]=y,v[y]++;
}
 
void dfs(int x){
    link(x)if(f[x]!=j->t){
	f[j->t]=x;dfs(j->t);
        if(d[x])s=max(s,d[j->t]+j->v+d[x]);
	d[x]=max(d[x],d[j->t]+j->v);
    }
}
 
int main(){
    n=read();
    inc(i,1,n){_x=read();_y=read();_t=read();add(_x,_y,_t);add(_y,_x,_t);}
    tar(1);c[tot+1]=c[1];inc(i,1,tot)link(c[i])if(j->t==c[i+1])a[i+1]=j->v+a[i];
    inc(i,1,tot)link(c[i])if(!v[j->t]){f[j->t]=c[i];dfs(j->t);if(d[c[i]])s=max(s,d[j->t]+j->v+d[c[i]]);d[c[i]]=max(d[c[i]],d[j->t]+j->v);}
    _x=0;inc(i,1,tot)pre[i]=max(pre[i-1],d[c[i]]+a[i]+_x),_x=max(_x,d[c[i]]-a[i]),b[i]=max(b[i-1],d[c[i]]+a[i]);
    _x=0;dec(i,tot,1)nxt[i]=max(nxt[i+1],d[c[i]]-a[i]+_x),_x=max(_x,d[c[i]]+a[i]),g[i]=max(g[i+1],d[c[i]]+a[tot+1]-a[i]);
    ans=pre[tot];
    inc(i,1,tot-1)ans=min(ans,max(pre[i],max(nxt[i+1],b[i]+g[i+1])));
    ans=max(ans,s);
    return 0*printf("%.1lf\n",ans/2.0);
}

3242: [Noi2013]快餐店

Time Limit: 20 Sec Memory Limit: 512 MB
Submit: 976 Solved: 545
[Submit][Status][Discuss]

Description

小T打算在城市C开设一家外送快餐店。送餐到某一个地点的时间与外卖店到该地点之间最短路径长度是成正比的，小T希望快餐店的地址选在离最远的顾客距离最近的地方。快餐店的顾客分布在城市C的N 个建筑中，这N 个建筑通过恰好N 条双向道路连接起来，不存在任何两条道路连接了相同的两个建筑。任意两个建筑之间至少存在一条由双向道路连接而成的路径。小T的快餐店可以开设在任一建筑中，也可以开设在任意一条道路的某个位置上（该位置与道路两端的建筑的距离不一定是整数）。现给定城市C的地图（道路分布及其长度），请找出最佳的快餐店选址，输出其与最远的顾客之间的距离。