BZOJ4762 最小集合（动态规划+容斥原理）

最新推荐文章于 2021-06-14 22:54:33 发布

weixin_30278311

最新推荐文章于 2021-06-14 22:54:33 发布

阅读量139

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/Gloid/p/10697754.html

本文分享了一段复杂的C++代码，该代码用于解决算法竞赛中的难题，涉及到数学运算、位操作、动态规划等高级编程技巧。代码中使用了多种C++标准库，如<iostream>、<cstdio>、<cmath>等，并定义了一系列实用的宏指令，如gcd、read等，以提高代码的效率和可读性。通过对代码的深入解析，读者可以学习到如何在有限的时间内解决复杂的算法问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

　　https://www.cnblogs.com/AwD-/p/6600650.html

#include<iostream> 
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
#define ll long long
#define P 1000000007
#define N 1024
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<'0'||c>'9')) c=getchar();return c;}
int gcd(int n,int m){return m==0?n:gcd(m,n%m);}
int read()
{
	int x=0,f=1;char c=getchar();
	while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
	while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x*f;
}
int n,a[N],f[2][N][N];
void inc(int &x,int y){x+=y;if (x>=P) x-=P;}
int main()
{
#ifndef ONLINE_JUDGE
	freopen("a.in","r",stdin);
	freopen("a.out","w",stdout);
	const char LL[]="%I64d\n";
#else
	const char LL[]="%lld\n";
#endif
	n=read();
	for (int i=1;i<=n;i++) a[i]=read();
	f[0][1023][1023]=1;
	for (int i=0;i<n;i++)
	{
		for (int j=0;j<N;j++)
			for (int k=N-1^j;k>=0;k=k==0?-1:(k-1&(N-1^j)))
			f[i&1^1][j][k|j]=f[i&1][j][k|j];
		for (int j=0;j<N;j++)
			for (int k=N-1^j;k>=0;k=k==0?-1:(k-1&(N-1^j)))
			inc(f[i&1^1][j&a[i+1]][(k|j)&a[i+1]],f[i&1][j][k|j]),
			inc(f[i&1^1][j&a[i+1]][((k|j)&a[i+1])|j],P-f[i&1][j][k|j]);
	}
	cout<<f[n&1][0][0]<<endl;
	return 0;
}

转载于:https://www.cnblogs.com/Gloid/p/10697754.html