[JSOI2016]无界单词

题目

题意：求( m border)长度为(0)的(n)位(0,1)字符串个数，并求字典序第(k)小的那一个。

首先是计数，正向不是很好算，考虑正难则反；设(f_i)表示长度为(i)的( m |border|=0)的串的个数

一个串可能有多个( m border)，我们考虑在其最小的( m border)长度(j)时计算它

则有(f_i=2^i-sum_{j=1}^{lfloor frac{i}{2} floor}f_j2^{i-2 imes j})；不用考虑长度大于(lfloor frac{i}{2} floor)的( m border)是因为当( m border)长度大于(lfloor frac{i}{2} floor)的时候一定存在一个小于(lfloor frac{i}{2} floor)的( m border)。

第二问，考虑大力二分；二分一个串，求所有字典序小于等于这个串中，( m border=0)的有多少个

设(g_{0,i})表示长度为(i)且字典序严格小于当前串的前(i)位的串的个数；(g_{1,i})表示长度为(i)且等于当前串的前(i)位的串的个数；不难发现(g_{1,i}=[i otin m Border])

还是跑上面那个容斥，从(g_{1,j})向(g_{0,i})转移的时候要在套一个二分判断一下

代码，其中有一些诡异的行为都是为了防爆unsigned long long

#include<bits/stdc++.h>
#define re register
#define LL unsigned long long
int n,T;LL m,f[66],v[66];
LL g[2][66];int nxt[66],a[66];
inline void split(LL nw) {
	for(re int i=1;i<=n;i++) a[i]=nw>>(n-i)&1;
}
inline LL Calc(LL val,LL R,int i,int j,LL H) {
	LL L=0,k=0;
	while(L<=R) {
		LL mid=(L>>1ull)+(R>>1ull);
		if((L&1ull)&&(R&1ull)) mid++;
		LL t=(((val<<i)|mid)<<j)|val;
		if(t<H) k=mid+1,L=mid+1;else {if(mid) R=mid-1;else break;}
	}
	return k;
}
inline LL calc(LL nw) {
	split(nw);nxt[1]=0;
	for(re int i=2;i<=n;i++) {
		int p=nxt[i-1];
		while(a[p+1]!=a[i]&&p) p=nxt[p];
		if(a[p+1]==a[i]) nxt[i]=p+1;else nxt[i]=0;
	}
	memset(g,0,sizeof(g));
	LL x=a[1];g[1][1]=1;
	if(x) g[0][1]=1;v[1]=x;
	for(re int i=2;i<=n;i++) {
		x=(x<<1ull)|a[i];v[i]=x;
		g[0][i]=x;
		for(re int j=1;j+j<=i;++j) {
			g[0][i]-=g[0][j]*(1ull<<(i-j-j));
			if(g[1][j]) g[0][i]-=Calc(v[j],(1ull<<(i-j-j))-1,i-j-j,j,x);
		}
		if(!nxt[i]) g[1][i]=1;
	}
	return g[0][n]+g[1][n];
}
int main() {
	scanf("%d",&T);f[1]=2;
	while(T--) {
		std::cin>>n>>m;
		for(re int i=2;i<=n;i++) {
			f[i]=1ull<<(i-1);f[i]--;f[i]+=1ull<<(i-1);
			for(re int j=1;2*j<=i;++j) f[i]-=(1ull<<(i-2*j))*f[j];
			f[i]++;
		}
		std::cout<<f[n]<<std::endl;
		LL l=0,r=(1ull<<(n-1)),ans=0;r+=r-1;
		while(l<=r) {
			LL mid=(l>>1ull)+(r>>1ull);
			if((l&1ull)&&(r&1ull)) mid++;
			if(calc(mid)>=m) {ans=mid;if(mid) r=mid-1;else break;}
			else l=mid+1;
		}
		split(ans);
		for(re int i=1;i<=n;i++)putchar(a[i]+'a');puts("");
	}
}