IOI2008 island

Posted zhou2003

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了IOI2008 island相关的知识,希望对你有一定的参考价值。

题目链接:[IOI2008]Island

题目大意:求基环树直径(由于题目的意思其实是类似于每个点只有一个出度,所以在每个联通块中点数和边数应该是相同的,这就是一棵基环树,所以题目给出的图就是一个基环树森林,又由于乘船的操作,可以知道答案就是将所有的直径相加起来的和)

分析:类似于树的直径,我们可以类推出基环树的直径也是有以下两种情况

? 1、是某个外向树的直径

? 2、两个点处于两课不同的外向树中(中间横跨了一段环上的边)

? 我们首先要把环找到

? 对于1的处理,我们可以对环上的点及其外向树直接进行dfs求解

? 对于2,我们会稍微复杂一点,首先可以明确的是,这两个点肯定是道他们对应的根节点距离最大的点,我们记以(u)为根节点的树中的最大距离是(dis_u),再在环上任选一个点,环上到这个点的距离为(sum_i)(所有的距离应是同向的),环上的每一条边的距离和为(sum_{cnt}),那么答案就应该是这样的

[ans=max(ans,max(dis_i+dis_j+(sum_i-sum_j),dis_i+dis_j+sum_{cnt}-(sum_i-sum_j)))]

将其变形可以得到

[ans=max(ans,max((sum_i+dis[i])-(sum_j-dis[j]),sum_{cnt}-(sum_i-dis[i])+(sum_j+dis[j])))]

注意到所有的(dis)(sum)均已知,且我们保证在处理节点(i)之前已经处理了点(j),并且通过我们的整理,上面的式子其实只剩下了两种形式:(sum_i+dis[i])(sum_i-dis[i])

那么为了保证这个最大值,我们可以记录(mind=min(sum_i-dis[i]))(maxd=max(sum_i+dis[i])),每次寻找(ans)时直接使用这两个值更新答案可以保证最优,同时沿途更新这两个值

#include<iostream>
#include<string.h>
#include<string>
#include<stdio.h>
#include<algorithm>
#include<vector>
#include<queue>
#include<map>
using namespace std;
struct node{
    int to,nxt,cost;
}sq[2000100];
int n,head[1000100],cp[1001000],fa[1001000],
    sq1[1001000][3],all=0,cnt=0;
bool vis[1001000],cir[1001000];
long long dis[1001000],sum[1001000],ans;

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;
}

void add(int u,int v,int w)
{
    all++;sq[all].to=v;sq[all].nxt=head[u];sq[all].cost=w;head[u]=all;
}

void init()
{
    n=read();
    int i;
    memset(head,0,sizeof(head));
    memset(vis,0,sizeof(vis));
    memset(cir,0,sizeof(cir));
    memset(dis,0,sizeof(dis));
    for (i=1;i<=n;i++)
    {
        int v=read(),w=read();
        add(i,v,w);add(v,i,w);
        sq1[i][0]=i;sq1[i][1]=v;sq1[i][2]=w;
    }
}

void findcircle(int u)
{
    vis[u]=1;cnt=0;
    while (1)
    {
        int v=sq1[u][1];
        if (vis[v])
        {
            cp[++cnt]=v;cir[v]=1;int i;sum[cnt]=sq1[v][2];
            for (i=u;i!=v;i=fa[i])
            {
                cir[i]=1;cp[++cnt]=i;sum[cnt]=sq1[i][2];
            }
            for (i=1;i<=cnt;i++) sum[i]+=sum[i-1];
            break;
        }
        else {vis[v]=1;fa[v]=u;}
        u=v;
    }
}

void find_d(int u)
{
    cir[u]=1;
    int i;
    for (i=head[u];i;i=sq[i].nxt)
    {
        int v=sq[i].to,w=sq[i].cost;
        if (cir[v]) continue;
        find_d(v);
        ans=max(ans,dis[u]+dis[v]+w);
        dis[u]=max(dis[u],dis[v]+w);
    }
}

long long findans(int u)
{
    sum[0]=0;
    findcircle(u);
    int i;
    ans=0;long long mind=1e18+7,maxd=-mind;
    for (i=1;i<=cnt;i++) find_d(cp[i]);
    for (i=1;i<=cnt;i++)
    {
        int v=cp[i];
        ans=max(ans,max(sum[i]+dis[v]-mind,sum[cnt]-(sum[i]-dis[v])+maxd));
        mind=min(mind,sum[i]-dis[v]);
        maxd=max(maxd,sum[i]+dis[v]);
    }
    //cout << 1;
    //for (i=1;i<=n;i++) cout << vis[i] << " ";cout << endl;
    return ans;
}
    

void work()
{
    long long sumans=0;
    int i;
    for (i=1;i<=n;i++) if (!cir[i]) sumans+=findans(i);
    printf("%lld",sumans);
}

int main()
{
    init();
    work();
    return 0;
}
/*
7
3 8
7 2
4 2
1 4
1 9
3 4
2 3
*/

以上是关于IOI2008 island的主要内容,如果未能解决你的问题,请参考以下文章

[IOI2008]Island

[题解] LuoguP4381 [IOI2008]Island

BZOJ 1791 [Ioi2008]Island 岛屿

[IOI2008]Island

[IOI2008]Island

[IoI2008]Island 岛屿