sos dp 题解备忘

用于遗忘sos dp后快速回忆和练习

题一

添加链接描述

题目大意

给你1e5个数（可重复），定义两个数的二进制位不包含即为一组，问给定的数中一共有多少组数满足这个条件。二进制位不包含即两个数都至少有一位是对方没有的。

分析

对于每一个数，只考虑比它大的数，就可以避免重复计算。因此只需要考虑有多少个数是包含这个数的，所以dp[i]表示：

(a[x]&i)==i的a[x]的个数

只需要用剩下的数减去这个数就是当前的数能够增加的组数，注意数可能重复。

代码

#include<bits/stdc++.h>
#define rep(i,n,m) for(int i=n;i<=m;i++)
#define repp(i,n,m) for(int i=n;i>=m;i--)
#define mod 1000000007
#define whlie while
#define itn int
#define icn cin
#define inc cin
#define mian main
typedef long long ll;
using namespace std;
ll a[1000005];
ll dp[2500005];
ll vis[2500005];
int main(){
	ios::sync_with_stdio(false);
	ll n,m;
	cin>>n>>m;
	rep(i,1,n){
		cin>>a[i];
		dp[a[i]]++;
		vis[a[i]]++;
	}
	rep(j,0,m-1){
		repp(i,(1<<m)-1,0){
			if(i&(1<<j)){
				dp[i^(1<<j)]=(dp[i]+dp[i^(1<<j)])%mod;
			}
		}	
	}	
	ll ans=0;
	sort(a+1,a+n+1);
	ll t=n;
	a[0]=-1;
	rep(i,1,n-1){
		if(a[i]!=a[i-1])t-=vis[a[i]];
//		cout<<t<<endl;
//		cout<<a[i]<<" "<<dp[a[i]]<<" ";
		ans+=t-dp[a[i]]+vis[a[i]];
//		cout<<ans<<endl;
	}
	cout<<ans<<endl;
	return 0;
}

题二

添加链接描述

题目大意

给出1e6个数（可重复），求每个数在这个数组中是否能找到另一个数与之进行"&"操作的答案是0，若能，则输出任意一个数，否则输出-1。

分析

本题要求寻找一个与原数没有同一位二进制都为1的数，记dp[i]为与i满足这样关系的任意一个数。易知：

与i满足这样关系的数和比i少一位二进制1的数一定也满足这样的关系

即可轻松得出转移方程

代码

#include<bits/stdc++.h>
#define rep(i,n,m) for(int i=n;i<=m;i++)
#define repp(i,n,m) for(int i=n;i>=m;i--)
#define whlie while
#define itn int
#define icn cin
#define inc cin
#define mian main
typedef int ll;
using namespace std;
ll a[1000005];
ll dp[10000005];
int main(){
	ios::sync_with_stdio(false);
	ll n;
	cin>>n;
//	memset(dp,0,sizeof(dp));
//	printf("%lld\n",(1<<23)-1);
	rep(i,1,n){
		cin>>a[i];
		dp[a[i]^((1<<23)-1)]=a[i];
//		printf("%lld\n",dp[a[i]^((1<<23)-1)]);
	}
	repp(i,((1<<23)-1),0){
		if(!dp[i]){
			rep(j,1,23){
				if(dp[i|(1<<j-1)]){
					dp[i]=dp[i|(1<<j-1)];
					break;
				}
			}
		}
	}
	rep(i,1,n){
		if(dp[a[i]])cout<<dp[a[i]]<<" ";
		else cout<<"-1"<<" ";
	}
	cout<<endl;
	return 0;
}

题三

添加链接描述

题目大意

给定1e6个数（可重复），求这些数中有多少个相互"&"操作后结果是0的组合，由于这个数可能非常大，输出答案%1e9+7的结果即可。

分析

明天再分析