#21. DFA

最新推荐文章于 2022-06-25 15:51:37 发布

wamach

最新推荐文章于 2022-06-25 15:51:37 发布

阅读量252

点赞数

分类专栏： simpleOJ

本文链接：https://blog.csdn.net/wamach/article/details/54882695

版权

simpleOJ 专栏收录该内容

25 篇文章 0 订阅

订阅专栏

首先，观察题目可知，该题不存在无解的情况，可以建一颗深度为n的满二叉树，以根为起点，收到0时向左走，收到1时向右走，再观察题目发现(1<<n+1)-1正好9192≈10000，无需证明，这就是对的。所以接下来只要把能够合并的点都合并起来，就能得到符合题意的自动机。

YYMHL（%A题大爷）说了，只要n^2枚举，两两合并，就可以A掉这题。时间复杂度O（kn^2），k为未知常数，听说很小。然后施大说正解就是这样的。顿时无语到不想鏼这题。

然后考试后的汪大就默默地想出了一个n（logn）^2的算法，怒艹标算！然后我就烧四五个小时把它膜了一遍，然后就没有然后了……

因此，蒟蒻决定重点分析这个算法。

先挂个代码，不然讲不清楚。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
#define rep(i,j,k) for(i=j;i<=k;++i)
#define per(i,j,k) for(i=j;i>=k;--i)
#define ll long long
#define ldb long double
#define pii pair<int,int>
#define mkp make_pair
#define X first
#define Y second
const int N=10005;
int n,nn,dep[N],ha[N],par[N],tot,u,vis[N],pos[N],lst[N],key[N];pii tmp[N];
struct TREE{int lc,rc;bool flg;}T[N];
void rebuild(int x){
	if(vis[x])return;vis[x]=1;
	pos[key[x]=++u]=x;if(T[x].flg)lst[++tot]=u;
	rebuild(T[x].lc=par[T[x].lc]),rebuild(T[x].rc=par[T[x].rc]);
}
int main(){
	int i,x,j,l;char str[14];
	scanf("%d",&n);nn=(1<<n+1)-1;
	rep(i,1,nn){
		scanf("%s",str);
		x=1;l=strlen(str)-1;
		rep(j,1,l)x=x<<1|str[j]=='1';
		if(x<<1<nn)T[x]=(TREE){x<<1,x<<1|1,str[0]=='+'};
		else T[x]=(TREE){x,x,str[0]=='+'};
	}
	rep(i,1,nn)dep[par[i]=i]=dep[i>>1]+1;
	for(l=nn,i=0;i<n;++i,l>>=1){
		rep(j,1,l)tmp[j]=mkp(ha[j]=ha[T[j].lc]+l*ha[T[j].rc]<<1|T[j].flg,j);
		sort(tmp+1,tmp+l+1);
		x=ha[tmp[1].Y]=0;
		rep(j,2,l)ha[tmp[j].Y]=x+=tmp[j].X>tmp[j-1].X;
		x=tmp[1].Y;
		rep(j,2,l)
			if(tmp[j].X==tmp[j-1].X){
				if(dep[tmp[j].Y]>n-i)par[tmp[j].Y]=x;
			}
			else x=tmp[j].Y;
	}
	rebuild(1);
	printf("%d 1\n%d",u,tot);
	rep(i,1,tot)printf(" %d",lst[i]);puts("");
	rep(i,1,u)printf("%d %d\n",key[T[pos[i]].lc],key[T[pos[i]].rc]);
	return 0;
}

另外附一张优先级的图（压行走火入魔了……）。