原题链接 https://www.luogu.org/problem/P5022



本着快csp了,做点往年的NOIp的题试试水来着,没想到水这么深 难度还挺大的,耗了我一天的时间(可能是我太菜了)
题目大意:
给你 n 个点和 m 条边,问如何遍历每个结点才能使最后的字典序最小,注意只有遍历完一棵字树的所有结点后才能回溯到他的父亲结点;
前 60 pts:
作为 NOIp2018 day2T1 来说,部分分确实给的挺足的,这 60 pts 就是哦;
看到 m = n-1 说明这是一棵树,考虑用搜索:
既然要保证是字典序最小,那么我们第一个点一定要选 1 了,考虑接下来只需要从他的儿子里按照编号从小到大遍历就好了,注意要遍历完一棵子树才能回溯上去;
后 40 pts:
m = n ?嗯,基环树!
所谓基环树,就是说在一个图里,有 n 个结点和 n 条边,那么这个图内有且仅有一个环;
我们仍然可以按照前 60 pts 的做法上去想:
按理来说我们遍历完 n 个结点只需要走 n-1 条边就好了啊,所以一定有一条边是多余的,也就是说我们根本遍历不到它,那么我们可以枚举每一条边,暂时把它删掉,然后和刚才的做法一样跑一遍 dfs,这样能求得一个字典序,我们取所有字典序中最小的一个就是答案了;
优化
显然如果我们像上述做法那样暴力的话,是会有几个点 TLE 的(毕竟NOIp也不会水到这种程度吧qwq),所以我们要考虑优化;
优化一:
题目中说了,我们的目的是遍历完每个点,那么假设我们删掉的边不是环上的边,那么这个图一定会变得不连通,那么就无法完成任务了;
所以我们只要找出基环树上的环,只需删掉环上的边就好了,至于找环嘛,这里我用的 tarjan(其实是只会这个qwq)
优化二:
最优性剪枝:假如我们在 dfs 的过程中,发现求得的字典序不如之前的答案优,那么我们直接可以 return 了;
有了这两个小优化,终于可以愉快的 AC 本题了,时间复杂度 O(n2);
至于更神仙的 O(n log n)甚至 O(n)的做法,咕咕咕~
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
int read()
{
char ch=getchar();
int a=0,x=1;
while(ch<'0'||ch>'9')
{
if(ch=='-') x=-x;
ch=getchar();
}
while(ch>='0'&&ch<='9')
{
a=(a<<1)+(a<<3)+(ch-'0');
ch=getchar();
}
return a*x;
}
const int N=5005;
int n,m,bi,bj,tot,top,tim,start,edge_sum,scc_sum,ans_sum=1;
int u[N],v[N],dfn[N],low[N],scc[N],st[N],vis[N],size[N],head[N],ans[N],f[N][N];
struct node
{
int to,next;
}a[N<<1];
void add(int from,int to) //链表存图
{
edge_sum++;
a[edge_sum].to=to;
a[edge_sum].next=head[from];
head[from]=edge_sum;
}
void tarjan(int u,int fa) //tarjan找环
{
dfn[u]=low[u]=++tim;
st[++top]=u;
vis[u]=1;
for(int i=head[u];i;i=a[i].next)
{
int v=a[i].to;
if(v==fa) continue;
if(!dfn[v])
{
tarjan(v,u);
low[u]=min(low[u],low[v]);
}
else if(vis[v]) low[u]=min(low[u],dfn[v]);
}
if(dfn[u]==low[u])
{
scc_sum++;
if(u==1) start=scc_sum;
while(st[top]!=u)
{
vis[st[top]]=0;
scc[st[top]]=scc_sum;
top--;
}
vis[st[top]]=0;
scc[st[top]]=scc_sum;
top--;
}
}
void dfs(int u,int fa)
{
tot++;
if(u<ans[tot]||bj) //如果比答案更优,或者之前就已经比答案优的话,就更新答案
{
bj=1; //发现更优解了
ans[tot]=u;
}
else if(u>ans[tot]) //最优性剪枝:不如答案优就直接返回
{
bi=1; //没有找到更优解
return ;
}
vis[u]=1;
for(int i=1;i<=n;i++) //从小到大去枚举每个点
{
if(i!=fa&&f[u][i]&&!vis[i]) //判断是否连接
{
dfs(i,u);
if(bi) return ; //剪枝
}
}
}
int main()
{
n=read();m=read();
for(int i=1;i<=m;i++)
{
u[i]=read();v[i]=read();
add(u[i],v[i]);add(v[i],u[i]); //链表存图,tarjan的时候方便
f[u[i]][v[i]]=1;f[v[i]][u[i]]=1; //邻接矩阵存图,选择贪心决策的时候方便从小到大找儿子
}
for(int i=1;i<=n;i++) //tarjan找环
{
if(!dfn[i]) tarjan(i,0);
}
memset(ans,0x3f,sizeof(ans));
if(m==n-1) //前60pts
{
dfs(1,0);
for(int i=1;i<=n;i++) printf("%d ",ans[i]);
return 0;
}
else //后40pts:基环树
{
for(int i=1;i<=m;i++)
{
if(scc[u[i]]==scc[v[i]]) //只需删除同一联通块里的边
{
tot=0;bi=0;bj=0; //bi表示是否一定不会找到更优解,bj表示是否找到了更优解
memset(vis,0,sizeof(vis));
f[u[i]][v[i]]=0; //暂时删边,这里用邻接矩阵就比较方便了
f[v[i]][u[i]]=0;
dfs(1,0); //更新答案
f[u[i]][v[i]]=1; //记得恢复
f[v[i]][u[i]]=1;
}
}
}
for(int i=1;i<=n;i++) printf("%d ",ans[i]);
return 0;
}