标签：想法 ring namespace next 细节 字典序 void mem algo
CSP-S2 2019 D1T3
考场上写了2h还是爆零……思维题还是写不来啊
思路分析
最开始可以想到最简单的贪心，从小到大枚举每个数字将其移动到最小的节点。但是通过分析样例后可以发现，一个数字在移动的过程中也可能有无关的边的删除，很难处理。显然直接贪心是不可能的。
分析删边对图的影响。可以发现，一条边删去之后，边两端的部分将不会产生任何影响。也就是说，两边的关系只有这一条边。于是还是之前那个贪心的想法，将边的问题转化为点的问题。现在来分析怎么求解。
具体实现
分析样例后可以发现以下性质：
因此，对于每个节点，关于与其相连的边，有以下三种约束条件：
我们可以根据上面提到的性质，通过判断每个节点的约束条件是否有冲突来求解方案。
可以发现，对于每个节点，其所有出边构成一条偏序链，我们可以通过维护一个链表来保证约束条件没有冲突。接下来分析一个数字经过一条边会产生什么冲突。假设某数字从$x$到$y$经过了边$(u,v)$。
偏序链可以用链表O(1)维护。另外，可以发现，每次数字的转移要维护的是一些连续的节点的关系，因此可以用一遍dfsO(n)维护。加上枚举数字的O(n)，总的时间复杂度O(n)。
细节比较多，注意不要漏点错点。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=2e3+100;
int T,n,tot;
int head[N],ver[2*N],Next[2*N];
int rev[N],pre[N],cnt1[N],cnt2[N],cnt3[N],from[N],to[N],p[N][N],header[N][N],tailer[N][N];
bool pd[N];
void add(int x,int y)
{
ver[++tot]=y,Next[tot]=head[x],head[x]=tot;
ver[++tot]=x,Next[tot]=head[y],head[y]=tot;
}
void check(int x,int root)//root就是数字的起点
{
for(int i=head[x],y=ver[i];i;i=Next[i],y=ver[i])
if(y!=pre[x])
{
pre[y]=x,pd[y]=1;
if(x!=root)
{
if(p[x][y]==x || !p[x][y])
pd[y]=0;
if(tailer[x][y]==from[x] && header[x][pre[x]]==to[x] && cnt1[x]*2+cnt2[x]+cnt3[x]-2>0)
pd[y]=0;
if(tailer[x][y]==pre[x])
pd[y]=0;
}//中间边
else
{
if(p[x][y]==x || !p[x][y])
pd[y]=0;
if(from[x] && tailer[x][y]==from[x] && cnt1[x]*2+cnt2[x]+cnt3[x]-1>0)
pd[y]=0;
}//起始边
pd[y]&=pd[x];//x不行y也不行
check(y,root);
}
if(x==root)
pd[x]=0;//起终点相同也不行
else
if(from[x] ||(to[x] && tailer[x][to[x]]==pre[x] && cnt1[x]*2+cnt2[x]+cnt3[x]-1>0))
pd[x]=0;//x作为终点
}
void clear()
{
memset(head,0,sizeof(head));
memset(Next,0,sizeof(Next));
memset(from,0,sizeof(from));//该节点上的数字从哪条边进来
memset(to,0,sizeof(to));//该节点上的数字从哪条边出去
memset(cnt1,0,sizeof(cnt1));//该节点还剩下几条双向没走过的边
memset(cnt2,0,sizeof(cnt2));//该节点还剩下几条单向出边
memset(cnt3,0,sizeof(cnt3));//该节点还剩下几条双向出边
memset(pd,0,sizeof(pd));//可行性
memset(p,0,sizeof(p));//-1表示没走过，0表示双向都走过，x表示以x为起点单向走过这条边，这里的初始化好像没什么用
memset(header,0,sizeof(header));//偏序链起始边
memset(tailer,0,sizeof(tailer));//偏序链终边
tot=0;
}
int main()
{
scanf("%d",&T);
while(T--)
{
clear();
scanf("%d",&n);
for(int i=1;i<=n;i++)
scanf("%d",&rev[i]);
for(int i=1,x,y;i<n;i++)
{
scanf("%d%d",&x,&y);
add(x,y);
cnt1[x]++,cnt1[y]++;
p[x][y]=p[y][x]=-1;
header[x][y]=tailer[x][y]=y,header[y][x]=tailer[y][x]=x;//初始值
}
for(int i=1,now;i<=n;i++)
{
for(int j=1;j<=n;j++)
pre[j]=0;
pd[rev[i]]=1;//初始值
check(rev[i],rev[i]);//dfs判断可行性
for(int j=1;j<=n;j++)
if(pd[j])
{
now=j;
break;
}//找到字典序最小的可行终点
printf("%d ",now);
from[now]=pre[now];
while(pre[now]!=rev[i])
{
if(p[pre[now]][now]==-1)
{
p[pre[now]][now]=p[now][pre[now]]=pre[now];
cnt1[now]--,cnt1[pre[now]]--,cnt3[now]++,cnt2[pre[now]]++;
}//双向没走过
else
{
p[pre[now]][now]=p[now][pre[now]]=0;
cnt2[now]--,cnt3[pre[now]]--;
}//反向走过
header[pre[now]][tailer[pre[now]][now]]=header[pre[now]][pre[pre[now]]];
tailer[pre[now]][header[pre[now]][pre[pre[now]]]]=tailer[pre[now]][now];//链表插入
now=pre[now];
}
if(p[pre[now]][now]==-1)
{
p[pre[now]][now]=p[now][pre[now]]=pre[now];
cnt1[now]--,cnt1[rev[i]]--,cnt3[now]++,cnt2[rev[i]]++;
}
else
{
p[pre[now]][now]=p[now][pre[now]]=0;
cnt2[now]--,cnt3[rev[i]]--;
}
to[rev[i]]=now;
}
puts("");
}
}
标签：想法 ring namespace next 细节 字典序 void mem algo
原文地址：https://www.cnblogs.com/TEoS/p/11969712.html
