CF388D-Fox and Perfect Sets【dp,线性基】

1|0正题

题目链接:https://www.luogu.com.cn/problem/CF388D

---

1|1题目大意

给出$k$，求有多少个集合$S$满足$S\sube [1,k]$且

$$a\in S,b\in S\Rightarrow a\ xor\ b\in S$$

$1\leq k\leq 10^9$

---

1|2解题思路

发现这个性质和线性基的很像，问题可以转换为问有多少个本质不同的线性基最大异或和不超过$k$。

在我们求$k$大异或和的时候我们把所有相同的线性基转换成了一个相同的形式满足对于它的数组$d$：

• $d_i$如果不是零那么其第$i$位是$1$且对于所有$j\neq i$都有$d_j$的第$i$位为$0$。

此时统计这种线性基的数量就不会算重了，并且全部的$d$异或起来一定是最大的。

然后然后就$dp$了，设$f_{i,j,0/1}$表示到二进制第$i$位，前面有$j$个$d$不是$0$，现在的异或和的前$i$位到达/每到达$k$的上限。

转移的时候对于一位$d$，如果填$0$那么就需要考虑前面的$d$有多少个加上这一位。为了方便转移可以预处理一下$g_i$表示$i$个东西选出奇数个的方案就好了。

时间复杂度$O(\log^2 k)$

---

1|3code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll P=1e9+7,N=35;
ll n,ans,f[N][N][2],C[N][N],g[N];
signed main()
{
	C[0][0]=1;
	for(ll i=1;i<N;i++)
		for(ll j=0;j<=i;j++)
			C[i][j]=(C[i-1][j]+(j?C[i-1][j-1]:0))%P;
	for(ll i=1;i<N;i++)
		for(ll j=1;j<=i;j+=2)g[i]+=C[i][j];
	scanf("%lld",&n);f[30][0][1]=1;
	for(ll i=29;i>=0;i--)
		for(ll j=0;j<30-i;j++){
			if((n>>i)&1){
				(f[i][j+1][1]+=f[i+1][j][1]%P)%=P;
				(f[i][j][1]+=f[i+1][j][1]*g[j]%P)%=P;
				(f[i][j+1][0]+=f[i+1][j][0]%P)%=P;
				(f[i][j][0]+=f[i+1][j][0]*(1ll<<j)%P+f[i+1][j][1]*((1ll<<j)-g[j])%P)%=P;
			}
			else{
				(f[i][j][1]+=f[i+1][j][1]*((1ll<<j)-g[j])%P)%=P;
				(f[i][j+1][0]+=f[i+1][j][0]%P)%=P;
				(f[i][j][0]+=f[i+1][j][0]*(1ll<<j)%P)%=P;
			}
		}
	for(ll j=0;j<30;j++)(ans+=f[0][j][0]+f[0][j][1])%=P;
	printf("%lld\n",ans);
	return 0;
}

__EOF__

本文作者：QuantAsk
本文链接：https://www.cnblogs.com/QuantAsk/p/15408312.html
关于博主：退役OIer，GD划水选手
版权声明：本博客所有文章除特别声明外，均采用 BY-NC-SA 许可协议。转载请注明出处！
声援博主：如果您觉得文章对您有帮助，可以点击文章右下角【推荐】一下。您的鼓励是博主的最大动力！