[清华集训2014]玛里苟斯

[清华集训2014]玛里苟斯

UOJ
BZOJ
对于k=1,全部或起来除个2(如果某个数该位为1,那么这一位为1的概率为\(\frac{1}{2}\))
对于k=2,枚举两位(如果全有1,概率为\(\frac{1}{2}\)或者\(\frac{1}{4}\),看是否同时出现),最后除个2
对于k>=3,发现\(a_i\leq2^{21}\),暴力枚举答案检验每个答案是否能被异或出来(线性基)

#define ull unsigned long long
#define ll long long
#include<bits/stdc++.h>
using namespace std;
const int _=1e5+5;
ll re(){
    ll x=0,w=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
    while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
    return x*w;
}
int n,k,cnt,s[_];
ll a[_];
ull ans;
void solve1(){
	for(int i=1;i<=n;i++)ans|=a[i];
	printf("%llu",ans>>1);
	if(ans&1)puts(".5");
}
void solve2(){
	for(int i=0;i<=32;i++)
		for(int j=0;j<=32;j++){
			int flag=1,A=0,B=0;
			for(int k=1;k<=n;k++){
				if(a[k]>>i&1)A=1;
				if(a[k]>>j&1)B=1;
				if((a[k]>>i&1)&&(a[k]>>j&1^1))flag=0;
				if((a[k]>>i&1^1)&&(a[k]>>j&1))flag=0;
			}
			if(!A||!B)continue;
			if(flag)ans+=(1ll<<(i+j));
			if(!flag&&i+j>0)ans+=(1ll<<(i+j-1));
		}
	printf("%llu",ans>>1);
	if(ans&1)puts(".5");
}
void insert(int x){
	for(int i=21;~i;i--){
		if(x>>i&1^1)continue;
		if(!s[i]){s[i]=x;return;}
		else x^=s[i];
	}
}
bool check(int x){
	for(int i=21;~i;i--)
		if(x>>i&1)x^=s[i];
	return x==0;
}
__int128 Pow(int x){
	__int128 res=1;for(int i=1;i<=k;i++)res*=x;return res;
}
void solve3(){
	__int128 Ans=0;
	for(int i=1;i<=n;i++)insert(a[i]);
	for(int i=0;i<=(1<<(63/k+1));i++)
		if(check(i))Ans+=Pow(i),cnt++;
	printf("%lld",Ans/cnt);
	if(Ans%cnt)puts(".5");
}
int main(){
    n=re(),k=re();
	for(int i=1;i<=n;i++)a[i]=re();
	if(k==1){solve1();return 0;}
	if(k==2){solve2();return 0;}
	solve3();return 0;
}
posted @ 2018-12-06 16:13  faced  阅读(243)  评论(0编辑  收藏  举报