其他分享
首页 > 其他分享> > NOIP2019 旅行

NOIP2019 旅行

作者:互联网

注意!注意!前方高能!本题卡常!!!

我们发现,所有的狗血剧情都在告诉我们,树的话直接dfs就出来了

那么基环树呢?

其实只要暴力删边,理论上的复杂度是可以过的qwq

但是删哪条边呢?

这里要引出一个基环树的常用操作:拓扑排序求环。具体方法是:在基环树上拓扑排序,然后拓扑序列中不存在的节点就是环中的节点了。

最后要用到环中的边的时候有一个小技巧,就是存边的时候(我用的是邻接表存双向边)按


input(x,y,z);
if(x>y) swap(x,y);
add(x,y,z);add(y,x,z);

的顺序存。

这样找边去重的时候就比较好找。。。(我也表述不太明白,具体看代码吧~)

Code


#include<iostream>
#include<cstdio>
#include<queue>

using namespace std;

const int N = 1e5+1;
int n,m,r[N];
int h[N],cnt,delu,delv;//假装 del_u-v 这条边已经被删去了qwq 
struct edge
{
    int nxt;
    int to;
}e[N];
int ans[N],ans1[N];
int qh=0,qt=1,qtp[N],tp[N],vis[N],visu[N],visv[N],tpcnt;//巨丑的代码qwq

inline void readx(int &x)
{
    x=0;char ch=getchar();
    while(ch<'0'||ch>'9') ch=getchar();
    while(ch>='0'&&ch<='9') {x=(x<<1)+(x<<3)+ch-48;ch=getchar();}
}

inline void add(register int u,register int v)
{
    e[++cnt].nxt=h[u];
    e[cnt].to=v;
    h[u]=cnt;
}

inline void topo()
{
    for(register int i=1;i<=n;++i) vis[i]=1;
    for(register int u=1;u<=n;++u)
    {
        for(register int i=h[u];i;i=e[i].nxt)
        {
            register int v=e[i].to;
            ++tp[u];++tp[v];
        }
    }
    for(register int i=1;i<=n;++i) if(tp[i]==2) qtp[qt++]=i;

//    for(int i=1;i<=n;++i) printf("%d ",qtp[i]);printf("\n");

    while(qt>qh)
    {
        register int u=qtp[++qh];vis[u]=0;
        for(register int i=h[u];i;i=e[i].nxt)
        {
            register int v=e[i].to;
            tp[u]-=2;tp[v]-=2;
            if(tp[v]==2&&vis[v]) qtp[qt++]=v;
        }
    }
    for(register int u=1;u<=n;++u)
        if(vis[u]==1)
        {
            for(register int i=h[u];i;i=e[i].nxt)
            {
                register int v=e[i].to;
                if(vis[v]==1&&u<v)
                {
                    visu[++tpcnt]=u;visv[tpcnt]=v;
                }
            }
        }
}

inline void update()
{
    for(register int i=1;i<=n;++i) 
    {
        if(ans1[i]==ans[i]) continue;
        if(ans[i]==0) {for(register int j=1;j<=n;++j) ans[j]=ans1[j];return;}
        if(ans1[i]>ans[i]) return;
        //ans1[i]<ans[i] 
        for(register int j=i;j<=n;++j) ans[j]=ans1[j];return;
    }
}

inline void dfs(register int u,register int fa)
{
    ans1[++cnt]=u;
    priority_queue <int,vector<int>,greater<int> > q;//小根堆 
    for(register int i=h[u];i;i=e[i].nxt)
    {
        register int v=e[i].to;
        if(v==fa||(u==delu&&v==delv)||(u==delv&&v==delu)) continue;
        q.push(v);
    }
    while(!q.empty())
    {
        register int v=q.top();q.pop();
        dfs(v,u);
    }
}

signed main()
{
    readx(n);readx(m);
    for(register int i=1;i<=m;++i)
    {
        int u,v;readx(u);readx(v);
        if(u>v) swap(u,v);
        add(u,v);add(v,u);
    }
    cnt=0;
    if(m==n-1)
    {
        dfs(1,0);
        for(register int i=1;i<=n;++i) ans[i]=ans1[i];
    }
    else
    {
        topo();
        for(register int i=1;i<=tpcnt;++i)
        {
            cnt=0;
            delu=visu[i];delv=visv[i];
            dfs(1,0);
            update();
        }
    }
    for(register int i=1;i<=n;++i)
        printf("%d ",ans[i]);
    return 0;
}

其实代码也就是看106行到107行那里,别的地方写得太丑了没法看,而且自己基本能写出来,加油哦qwq!

标签:旅行,ch,int,register,NOIP2019,add,&&,qwq
来源: https://www.cnblogs.com/oierwyh/p/11825256.html