旅行 Noip 2018 Day2T1

梦想与她 提交于 2019-11-30 04:30:13

突然发现好久没写博客了,那就讲一讲Noip的题吧

-----------------------------------我是分割线---------------------------------------

直接上题:

题目描述

小 Y 是一个爱好旅行的 OIer。她来到 X 国,打算将各个城市都玩一遍。

小Y了解到, X国的 n 个城市之间有 m 条双向道路。每条双向道路连接两个城市。 不存在两条连接同一对城市的道路,也不存在一条连接一个城市和它本身的道路。并且, 从任意一个城市出发,通过这些道路都可以到达任意一个其他城市。小 Y 只能通过这些 道路从一个城市前往另一个城市。

小 Y 的旅行方案是这样的:任意选定一个城市作为起点,然后从起点开始,每次可 以选择一条与当前城市相连的道路,走向一个没有去过的城市,或者沿着第一次访问该 城市时经过的道路后退到上一个城市。当小 Y 回到起点时,她可以选择结束这次旅行或 继续旅行。需要注意的是,小 Y 要求在旅行方案中,每个城市都被访问到。

为了让自己的旅行更有意义,小 Y 决定在每到达一个新的城市(包括起点)时,将 它的编号记录下来。她知道这样会形成一个长度为 nn 的序列。她希望这个序列的字典序 最小,你能帮帮她吗? 对于两个长度均为 nn 的序列 AA 和 BB,当且仅当存在一个正整数 xx,满足以下条件时, 我们说序列 AA 的字典序小于 BB。

  • 对于任意正整数 1 ≤ i < x1i<x,序列 A 的第 i 个元素 Ai 和序列 B 的第 i 个元素 B_iBi 相同。
  • 序列 A 的第 x 个元素的值小于序列 B 的第 x 个元素的值。

输入格式

输入文件共 m + 1m+1 行。第一行包含两个整数 n,m(m ≤ n)n,m(mn),中间用一个空格分隔。

接下来 m 行,每行包含两个整数 u,v (1 ≤ u,v ≤ n)u,v(1u,vn) ,表示编号为 uu 和 vv 的城市之 间有一条道路,两个整数之间用一个空格分隔。

输出格式

输出文件包含一行,nn 个整数,表示字典序最小的序列。相邻两个整数之间用一个 空格分隔。

输入输出样例

输入 #1
6 5 
1 3 
2 3 
2 5 
3 4 
4 6
输出 #1
1 3 2 5 4 6
输入 #2
6 6 
1 3 
2 3 
2 5 
3 4 
4 5 
4 6
输出 #2
1 3 2 4 5 6
考场上的心理:看起来还行,不难,不断dfs找编号小的点。long time later,代码实现:好像是错的,这什么玩意儿,凉凉回归正题,这道题目该怎么做:回到题目,看到提示和说明:m=n-1|m=n;看起来好像没有什么用,但仔细想想啊当m=n时,图不就退化成为一棵树了吗(暴力dfs,相信看到人都能够熟练的运用\手动滑稽\)但是,第二种情况:m=n-1,这看起来就很麻烦了,因为有环,且只有一个环那怎么办呢,好了,这是这道题目最难的一点,也是最难想到的(反正我是这样)既然有且只有一个环,那么妥妥的奇环树我么有一个专业名称:断环法,听起来好像特别高级,但其实通俗的讲就是暴力删边重要的事情说三遍,不是每个边,只有环上的边*3\滑稽\那么找环,就要拿出Tarjan,上一篇有详细讲过第二种就是拓扑排序 topsort(应该也不用说)那么找环代码如下:
void topsort(){
    int hd=0,tl=0;
    for(int i=1;i<=number;i++)
       if(in[i]==1)que[++tl]=i;
    while(hd<tl){
        int now=que[++hd];
        for(int i=head[now];i;i=nxt[i]){
            int y=ver[i];
            if(in[y]>1){
                in[y]--;
                if(in[y]==1)que[++tl]=y;
            }
        }
    }
}

que是不在环上的点,in表示入度,那么环上的点就是入度>=2的点

找到环以后,就开始删边了

将环上任意两个相邻点标记,mapp表示,代表这两个点之间的边被强行删去

那么删除之后,发现不又回到了m=n-1的情况吗

那么,最后只有判断当前的序列是否更优

代码如下

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

const int N=30002;
int number,magic,ans[N],res[N],sum=0,need[N];
int ver[N],nxt[N],tot=0,head[N],size[N];
int in[N],que[N],mapp[5002][5002];

int read(){
	int s=0,w=1;char ch=getchar();
	while(ch<'0'||ch>'9')w=(ch=='-')?-1:1,ch=getchar();
	while(ch>='0'&&ch<='9')s=s*10+ch-'0',ch=getchar();
	return s*w;
}

void add(int x,int y){
	ver[++tot]=y;nxt[tot]=head[x];head[x]=tot;
	in[x]++;in[y]++;
}

void dfs(int now,int fa){
	ans[++sum]=now;
	int q[3002]={0},vis=0;
	for(int i=head[now];i;i=nxt[i]){
		int y=ver[i];
		if(y==fa||mapp[now][y])continue;
		q[++vis]=y;
	}
	sort(q+1,q+vis+1);
	for(int i=1;i<=vis;i++)dfs(q[i],now);
}

void topsort(){
	int hd=0,tl=0;
	for(int i=1;i<=number;i++)
	   if(in[i]==1)que[++tl]=i;
	while(hd<tl){
		int now=que[++hd];
		for(int i=head[now];i;i=nxt[i]){
			int y=ver[i];
			if(in[y]>1){
				in[y]--;
				if(in[y]==1)que[++tl]=y;
			}
		}
	}
}

void check(){
	for(int i=1;i<=number;i++)
	   if(ans[i]>need[i])return;
	for(int i=1;i<=number;i++)
	   need[i]=ans[i];
}

void find(int now){
	int x=now,vis=0,xl[3002],last=0;
	do{
		in[x]=1;xl[++vis]=x;last=0;
		for(int i=head[x];i;i=nxt[i]){
			int y=ver[i];
			if(in[i]>1){
				x=y;last=i;break;
			}
		}
	}while(last);
	for(int i=1;i<vis;i++){
		mapp[xl[i]][xl[i+1]]=mapp[xl[i+1]][xl[i]]=1;
		dfs(1,0);
		check();
		mapp[xl[i]][xl[i+1]]=mapp[xl[i+1]][xl[i]]=0;
	}
}

int main(){
	number=read();magic=read();
	for(int i=1;i<=number;i++)need[i]=1e9;
	for(int i=1;i<=magic;i++){
		int x=read(),y=read();
		add(x,y);add(y,x);
	}
	if(number-1==magic)dfs(1,0);
	else {
		topsort();
		for(int i=1;i<=number;i++)
		   if(in[i]>1)find(i);
	}
	for(int i=1;i<=number;i++)
	   cout<<ans[i]<<" ";
	return 0;
}

 

  

 

当然,还没有结束,当你把当前的代码交上去时发现如下:T到起飞,(自动屏蔽wa,输出错了,应该是need,我是ans)那么原因是什么呢,看到
void dfs(int now,int fa){
    ans[++sum]=now;
    int q[3002]={0},vis=0;
    for(int i=head[now];i;i=nxt[i]){
        int y=ver[i];
        if(y==fa||mapp[now][y])continue;
        q[++vis]=y;
    }
    sort(q+1,q+vis+1);
    for(int i=1;i<=vis;i++)dfs(q[i],now);
}

用了sort,可以卡到 O(n^3*log2n)

看看数据范围:n<=5000  接近崩溃的边缘

不过别慌,要用优化来打败TLE————沃兹基硕德

我们发现,超时的主要原因是不知道一个点连接边点的大小顺序,所以要先存起来

那么,我们可以预处理

现将连接着的点排序(程序开始存,不然没有效果)

注意:相信很多人排序都会从小到大排序,但是前项星存是倒过来的,所以要从大到小排

所以AC代码如下:

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

const int M=5050,N=M<<1;
int number,edge_n,ans[M],res[M],sum=0,need[M];
int ver[N],nxt[N],tot=0,head[N],size[N];
int in[M],que[M],mapp[M][M];
struct node{
    int x,y;
}edge[N];

int read(){
    int s=0,w=1;char ch=getchar();
    while(ch<'0'||ch>'9')w=(ch=='-')?-1:1,ch=getchar();
    while(ch>='0'&&ch<='9')s=s*10+ch-'0',ch=getchar();
    return s*w;
}

void add(int x,int y){
    ver[++tot]=y;nxt[tot]=head[x];head[x]=tot;
    in[y]++;
}

void dfs(int now){
    ans[++sum]=now;res[now]=1;
    for(int i=head[now];i;i=nxt[i]){
        int y=ver[i];
        if(res[y]||mapp[now][y])continue;
        dfs(y);
    }
}

void topsort(){
    int hd=0,tl=0;
    for(int i=1;i<=number;i++)
       if(in[i]==1)que[++tl]=i;
    while(hd<tl){
        int now=que[++hd];
        for(int i=head[now];i;i=nxt[i]){
            int y=ver[i];
            if(in[y]>1){
                in[y]--;
                if(in[y]==1)que[++tl]=y;
            }
        }
    }
}

void check(){
    int flag=0;
    for(int i=1;i<=number;i++){
        if(ans[i]<need[i]){
            flag=1;break;
        }
        if(ans[i]>need[i])return;
    }
    if(!flag)return;
    for(int i=1;i<=number;i++)
       need[i]=ans[i];
}

void find(int now){
    int x=now,vis=0,xl[5002],last=0;
    do{
        in[x]=1;xl[++vis]=x;last=0;
        for(int i=head[x];i;i=nxt[i]){
            int y=ver[i];
            if(in[y]>1){
                x=y;last=y;break;
            }
        }
    }while(last);
    xl[++vis]=now;
    for(int i=1;i<vis;i++){
        mapp[xl[i]][xl[i+1]]=mapp[xl[i+1]][xl[i]]=1;
        memset(res,0,sizeof(res));sum=0;
        dfs(1);
        check();
        mapp[xl[i]][xl[i+1]]=mapp[xl[i+1]][xl[i]]=0;
    }
}

bool cmp(node a,node b){
    return a.y>b.y;
}

int main(){
    number=read();edge_n=read();
    for(int i=1;i<=number;i++)need[i]=1e9;
    for(int i=1;i<=edge_n;i++){
        int x=read(),y=read();
        edge[i].x=x;edge[i].y=y;
        edge[i+edge_n].x=y;edge[i+edge_n].y=x;
    }
    sort(edge+1,edge+edge_n*2+1,cmp);
    for(int i=1;i<=edge_n*2;i++)
       add(edge[i].x,edge[i].y);
    if(number-1==edge_n){
        dfs(1);check();
    }
    if(number==edge_n){
        topsort();
        for(int i=1;i<=number;i++)
           if(in[i]>1){
               find(i);break;
           }
    }
    for(int i=1;i<=number;i++)
       cout<<need[i]<<" ";
    return 0;
}

CSP 2019 rp++

 

标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!