ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

「NOI2014」魔法森林

2019-07-26 19:03:33  阅读:233  来源: 互联网

标签:Index int Max 魔法 Father NOI2014 Tag Child 森林


题目链接

题目分析

实则毫无思路……

由于权值是两维的,首先尝试正常操作排序一维。然后问题就变成了求另外一维的最小。不难想到另外一维是一个最小生成树。每次的答案就是\(1\)到\(n\)路径上的两个最大加起来(不一定是枚举的一维加上另一维的最大)。

然后考虑如何维护。动态加边的最小生成树自然用LCT维护。由于LCT维护的是点值,所以我们要把边拆成一个点加两条边,点上维护原来边上的值。然后直接套LCT维护max即可。

参考程序

由于尝试封装了一下,所以代码巨翔无比……大家可以直接跳过LCT部分。如果要学习LCT的可以戳这里

#include <bits/stdc++.h>
using namespace std;

    struct LCTChild {
        int Child[ 2 ];
    };
    struct LCT {
        int *Father, *Stack, *Tag;
        LCTChild *Child;
        LCT() {}
        LCT( int *_Father, LCTChild *_Child, int *_Stack, int *_Tag ) : Father( _Father ), Child( _Child ), Stack( _Stack ), Tag( _Tag ) {}
        void Clear() { memset( Father, 0, sizeof( Father ) ); memset( Tag, 0, sizeof( Tag ) ); memset( Child, 0, sizeof( Child ) ); return; }
        void Collect( int Index );
        void TagDown( int x );
        bool IsRoot( int x ) { return !( ( Child[ Father[ x ] ].Child[ 0 ] == x ) || ( Child[ Father[ x ] ].Child[ 1 ] == x ) ); }
        void Rotate( int C ) {
            int B = Father[ C ]; int A = Father[ B ]; int Tag = Child[ B ].Child[ 1 ] == C;
            if( !IsRoot( B ) ) Child[ A ].Child[ Child[ A ].Child[ 1 ] == B ] = C; Father[ C ] = A;
            Father[ Child[ C ].Child[ Tag ^ 1 ] ] = B; Child[ B ].Child[ Tag ] = Child[ C ].Child[ Tag ^ 1 ];
            Child[ C ].Child[ Tag ^ 1 ] = B; Father[ B ] = C;
            Collect( B ); Collect( C );
            return;
        }
        void Splay( int x ) {
            int Num = 0; Stack[ ++Num ] = x; int t = x; for( ; !IsRoot( t ); t = Father[ t ] ) Stack[ ++Num ] = Father[ t ]; for( int i = Num; i >= 1; --i ) TagDown( Stack[ i ] );
            while( !IsRoot( x ) ) { 
                int y = Father[ x ]; int z = Father[ y ]; if( !IsRoot( y ) ) if( ( Child[ z ].Child[ 0 ] == y ) ^ ( Child[ y ].Child[ 0 ] == z ) ) Rotate( x ); else Rotate( y ); Rotate( x );
            }
            Collect( x );
            return;
        }
        void Access( int x ) { for( int i = 0; x; i = x, x = Father[ x ] ) { Splay( x ); Child[ x ].Child[ 1 ] = i; Collect( x ); } return; }
        void MakeRoot( int x ) { Access( x ); Splay( x ); Tag[ x ] ^= 1; return; }
        int FindRoot( int x ) { Access( x ); Splay( x ); TagDown( x ); while( Child[ x ].Child[ 0 ] ) { x = Child[ x ].Child[ 0 ]; TagDown( x ); } Splay( x ); return x; }
        void Split( int x, int y ) { MakeRoot( x ); Access( y ); Splay( y ); return; }
        void Link( int x, int y ) { 
            MakeRoot( x ); if( FindRoot( y ) == x ) return; Father[ x ] = y; return; 
        }
        void Cut( int x, int y ) {
            MakeRoot( x ); if( FindRoot( y ) != x || Child[ y ].Child[ 0 ] || Father[ y ] != x ) return; Father[ y ] = Child[ x ].Child[ 1 ] = 0; Collect( x ); return; 
        }
    };
const int Maxn = 50010;
const int Maxm = 100010;
struct edge {
    int X, Y, A, B;
};
edge Edge[ Maxm ];
int N, M;
LCTChild Child[ Maxn + Maxm ];
int Father[ Maxn + Maxm ], Stack[ Maxn + Maxm ], Tag[ Maxn + Maxm ];
LCT Lct( Father, Child, Stack, Tag );
struct rec {
    int A, B;
};
rec Max[ Maxn + Maxm ], Value[ Maxn + Maxm ];
int From[ Maxn + Maxm ];
void LCT::Collect( int Index ) {
    Max[ Index ].A = max( Max[ Child[ Index ].Child[ 0 ] ].A, Max[ Child[ Index ].Child[ 1 ] ].A );
    Max[ Index ].B = max( Max[ Child[ Index ].Child[ 0 ] ].B, Max[ Child[ Index ].Child[ 1 ] ].B );
    if( Max[ Child[ Index ].Child[ 0 ] ].B > Max[ Child[ Index ].Child[ 1 ] ].B ) From[ Index ] = From[ Child[ Index ].Child[ 0 ] ];
    else From[ Index ] = From[ Child[ Index ].Child[ 1 ] ];
    Max[ Index ].A = max( Max[ Index ].A, Value[ Index ].A );
    if( Value[ Index ].B > Max[ Index ].B ) From[ Index ] = Index;
    Max[ Index ].B = max( Max[ Index ].B, Value[ Index ].B );
    return; 
}

void LCT::TagDown( int x ) {
    if( Tag[ x ] ) { 
        Tag[ x ] = 0; 
        std::swap( Child[ x ].Child[ 0 ], Child[ x ].Child[ 1 ] );
        Tag[ Child[ x ].Child[ 0 ] ] ^= 1; Tag[ Child[ x ].Child[ 1 ] ] ^= 1; 
    } 
    return;
}

bool Cmp( edge X, edge Y ) {
    return X.A < Y.A || X.A == Y.A && X.B < Y.B;
}

int main() {
    scanf( "%d%d", &N, &M );
    for( int i = 1; i <= M; ++i )
        scanf( "%d%d%d%d", &Edge[ i ].X, &Edge[ i ].Y, &Edge[ i ].A, &Edge[ i ].B );
    std::sort( Edge + 1, Edge + M + 1, Cmp );
    for( int i = 1; i <= M; ++i ) {
        Value[ i + N ].A = Max[ i + N ].A = Edge[ i ].A;
        Value[ i + N ].B = Max[ i + N ].B = Edge[ i ].B;
        From[ i + N ] = i + N;
    }
    int Ans = -1;
    for( int i = 1; i <= M; ++i ) {
        Lct.MakeRoot( Edge[ i ].X );
        if( Lct.FindRoot( Edge[ i ].Y ) != Edge[ i ].X ) {
            Lct.Link( Edge[ i ].X, i + N );
            Lct.Link( i + N, Edge[ i ].Y );
        } else {
            Lct.Split( Edge[ i ].Y, Edge[ i ].X );
            int T = From[ Edge[ i ].X ] - N;
            if( Edge[ T ].B > Edge[ i ].B ) {
                Lct.Cut( Edge[ T ].X, T + N );
                Lct.Cut( T + N, Edge[ T ].Y );
                Lct.Link( Edge[ i ].X, i + N );
                Lct.Link( i + N, Edge[ i ].Y );
            }
        }
        Lct.MakeRoot( 1 );
        if( Lct.FindRoot( N ) == 1 ) {
            Lct.Split( N, 1 );
            Ans = ( Ans == -1 ) ? ( Max[ 1 ].A + Max[ 1 ].B ) : min( Ans, Max[ 1 ].A + Max[ 1 ].B );
        }
    }
    printf( "%d\n", Ans );
    return 0;
}

标签:Index,int,Max,魔法,Father,NOI2014,Tag,Child,森林
来源: https://www.cnblogs.com/chy-2003/p/11252181.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有