day2就蓝题起步不是要我命吗
目测数据:不都是基环树和树嘛?基环树我不会但是这个树的贪心不是裸的吗?
大手一挥:这部分分我打定了!
int outp[5100],len;
void dfs(int x,int fa)
{
outp[++len]=x;
list<int>q;//链表……不要问我为什么要用这玩意儿(因为没用过想试一试
for(int k=last[x];k;k=a[k].next)
if(a[k].y!=fa)
q.push_back(a[k].y);
q.sort();
for(list<int>::iterator i=q.begin();i!=q.end();i++)
dfs(*i,x);
}
void Main()
{
dfs(1,0);
for(int i=1;i<=n;i++)
printf("%d ",outp[i]);
}
光荣地水到了60分
对于剩下的分,我们来好好观察一下题面:
n
≤
5000
n\le 5000
n≤5000
常数小的
n
2
log
n
n^2\log n
n2logn都可以过啊!
容 易 证 明整棵基环树的环中一定会有一条边是不需要走过的
博主只记了结论所以以后再来填坑啦
所以我们通过删边可以把基环树转化为一棵树,然后直接跑第一种做法
首先我们要找出这个环
通过拓扑排序解决
首先,由于是双向连通图,所以每个点的入度至少为1
我们就从入度为1的点开始,进行宽搜
观察一下样例:
可以看到,3和4这两个“环的入口”度数都为3
所以,任意一个点触碰到环上的点时,都无法将环上的点入度减为1
最后,我们就可以将这些入度不为零的点计入环中。
拓扑排序:
void topsort()
{
queue<int>q;
for(int i=1;i<=n;i++)
if(in_degree[i]==1)q.push(i);
while(q.size())
{
int x=q.front();
for(int k=last[x];k;k=a[k].next)
if(--in_degree[a[k].y]==1)
q.push(a[k].y);
q.pop();
}
}
还真就音译topsort呗
主函数中部分:
topsort();
for(int i=1;i<=n;i++)
on_cycle[i]=in_degree[i]>1;
然后,我们就按照这些点排出环上的边:
for(int i=1;i<=n;i++)
if(on_cycle[i])
{
get_edges(i,i,0);
break;
}
get_edges
:
bool get_edges(int x,int ed,int fa)
{
for(int k=last[x],y=a[k].y;k;k=a[k].next,y=a[k].y)
if(on_cycle[y]&&y!=fa)
{
cycle[++cycle_len]=k;
if(y==ed)return 1;
if(get_edges(y,ed,x))return 1;//这个地方是为了防止st从两边同时访问这个环,造成重复计数
}
}
然后就一根一根地断,将cut
标记置为1
不要忘记双向边
for(int i=1;i<=cycle_len;i++)
{
cut[cycle[i]]=cut[cycle[i]^1]=1;//^1是对向边,为了这句话要将len提前设为1
top=len=0;dfs(1,0);
if(outp[1]==0||comp(tmp,outp))memcpy(outp+1,tmp+1,n<<2);
cut[cycle[i]]=cut[cycle[i]^1]=0;//别忘了回溯
}
dfs
:
void dfs(int x,int fa)
{
const int last_top=top;
tmp[++len]=x;
for(int k=last[x];k;k=a[k].next)
if(a[k].y!=fa&&!cut[k])
sta[++top]=a[k].y;//sta作为一个排序用的数组,防止爆栈又不怕STL常数大
const int now_top=top;
sort(sta+last_top+1,sta+now_top+1);//贪心
for(int i=last_top+1;i<=now_top;i++)
dfs(sta[i],x);
}
comp
:
bool comp(int a[],int b[])
{
for(int i=1;i<=n;i++)
if(a[i]!=b[i])
if(a[i]<b[i])return 1;
else return 0;
return 0;
}
(本来为了省掉comp
用了basic_string<int>
,结果被其他地方的RE莫名劝退了)
最后outp
数组就是答案了
代码:
#include<cstdio>
#include<cstring>
#include<queue>
#include<list>
#include<algorithm>
using namespace std;
inline void read(int &x)
{
x=0;register char c=getchar();
while(c<48||57<c) c=getchar();
for(;48<=c&&c<=57;c=getchar()) x=x*10+(c&15);
}
int n,m;
struct edge
{
int y,next;
}a[11000];
int last[5100],Len;
int in_degree[5100];
inline void ins(int x,int y)
{
a[++Len]=(edge){y,last[x]};
last[x]=Len;
in_degree[y]++;
}
int outp[5100];
namespace solve1
{
int len;
void dfs(int x,int fa)
{
outp[++len]=x;
list<int>q;
for(int k=last[x];k;k=a[k].next)
if(a[k].y!=fa)
q.push_back(a[k].y);
q.sort();
for(list<int>::iterator i=q.begin();i!=q.end();i++)
dfs(*i,x);
}
void Main()
{
dfs(1,0);
for(int i=1;i<=n;i++)
printf("%d ",outp[i]);
}
};
namespace solve2
{
bool on_cycle[5100];
bool cut[11000];
int cycle[5100],cycle_len;
int tmp[5100],len;
int sta[5100],top;
void topsort()
{
queue<int>q;
for(int i=1;i<=n;i++)
if(in_degree[i]==1)q.push(i);
while(q.size())
{
int x=q.front();
for(int k=last[x];k;k=a[k].next)
if(--in_degree[a[k].y]==1)
q.push(a[k].y);
q.pop();
}
}
bool get_edges(int x,int ed,int fa)
{
for(int k=last[x],y=a[k].y;k;k=a[k].next,y=a[k].y)
if(on_cycle[y]&&y!=fa)
{
cycle[++cycle_len]=k;
if(y==ed)return 1;
if(get_edges(y,ed,x))return 1;
}
}
void dfs(int x,int fa)
{
const int last_top=top;
tmp[++len]=x;
for(int k=last[x];k;k=a[k].next)
if(a[k].y!=fa&&!cut[k])
sta[++top]=a[k].y;
const int now_top=top;
sort(sta+last_top+1,sta+now_top+1);
for(int i=last_top+1;i<=now_top;i++)
dfs(sta[i],x);
}
bool comp(int a[],int b[])
{
for(int i=1;i<=n;i++)
if(a[i]!=b[i])
if(a[i]<b[i])return 1;
else return 0;
return 0;
}
void Main()
{
topsort();
for(int i=1;i<=n;i++)
on_cycle[i]=in_degree[i]>1;
for(int i=1;i<=n;i++)
if(on_cycle[i])
{
get_edges(i,i,0);
break;
}
for(int i=1;i<=cycle_len;i++)
{
cut[cycle[i]]=cut[cycle[i]^1]=1;
top=len=0;dfs(1,0);
if(outp[1]==0||comp(tmp,outp))memcpy(outp+1,tmp+1,n<<2);
cut[cycle[i]]=cut[cycle[i]^1]=0;
}
for(int i=1;i<=n;i++)
printf("%d ",outp[i]);
}
};
int main()
{
read(n);read(m);Len=1;
int x,y;
for(int i=1;i<=m;i++)
{
scanf("%d%d",&x,&y);
ins(x,y);
ins(y,x);
}
if(n==m)
solve2::Main();
else solve1::Main();
return 0;
}