题解:Codeforces Round #779 (Div. 2) D1、D2(位运算,奇淫技巧)

本文D1,D2做法完全相同,这里以D2为例进行讲解

题目链接:D2. 388535 (Hard Version)

题意:给定区间[l,r],(0≤l≤r<2^17),和长为r-l+的数组a,数组a是由[l,r]全排列组成的某种排列全部异或x后形成.,输出任意一个可能的x

题解:本文题解属于奇淫技巧,仅供参考,首先我们可以发现,某个正解x一定存在于所有的a[i]^l中,我们可以把所有a[i]^l存到一个数组c中作为可选数,如何进一步枚举呢?我们可以发现,如果x是正确的那么\sum_{l}^{r}a[i]^x==\sum_{l}^{r}i,所有我们可以预处理出数组a中每个二进制上1的个数,在枚举所有a[i]^l时可以通过a[i]^l每一位的01情况确定\sum_{l}^{r}a[i]^x,因为能够满足最高条件的情况非常非常少(具体为啥也不知道,靠感觉),所以对于符合的就可以直接暴力判断

代码:D1,D2都适用

#include<iostream>
#include<stack>
#include<list>
#include<set>
#include<vector>
#include<algorithm>
#include<math.h>
#include<numeric>
#include<map>
#include<cstring>
#include<queue>
#include<iomanip>
#include<cmath>
#include<queue>
#include <bitset>
#include<unordered_map>/*
	#ifndef local
	#define endl '\n'
#endif */
#define mkp make_pair
using namespace std;
using std::bitset;
typedef long long ll;
typedef long double ld;
const int inf=0x3f3f3f3f;
const ll MAXN=2e6+10;
const ll N=1e6+100;
const ll mod=998244353;
const ll hash_p1=1610612741;
const ll hash_p2=805306457;
const ll hash_p3=402653189;
//-----------------------------------------------------------------------------------------------------------------*/
// ll head[MAXN],net[MAXN],to[MAXN],edge[MAXN]/*流量*/,cost[MAXN]//费用;
/* 
void add(ll u,ll v,ll w,ll s){
	to[++cnt]=v;net[cnt]=head[u];edge[cnt]=w;cost[cnt]=s;head[u]=cnt;
	to[++cnt]=u;net[cnt]=head[v];edge[cnt]=0;cost[cnt]=-s;head[v]=cnt;
}
struct elemt{
	int p,v;
};
struct comp{
	public:
		bool operator()(elemt v1,elemt v2){
			return v1.v<v2.v;
		}
};
-----------------------------------
求[1,MAXN]组合式和逆元 
ll mi(ll a,ll b){
	ll res=1;
	while(b){
		if(b%2){
			res=res*a%mod;
		}	
		a=a*a%mod;
	}
	return res;
}
ll fac[MAXN],inv[MAXN]
fac[0]=1;inv[0]=1;
for(int i=1;i<=MAXN;i){
	fac[i]=(fac[i-1]*i)%mod;
	inv[i]=mi(fac[i],mod-2);
}
ll C(int m,int n){//组合式C(m,n); 
	if(!n){
		return 1;
	}
	return fac[m]*(inv[n]*inv[m*-n]%mod)%mod;
}
---------------------------------
 unordered_map<int,int>mp;
//优先队列默认小顶堆 , greater<int> --小顶堆  less<int> --大顶堆  
priority_queue<elemt,vector<elemt>,comp>q;
	set<int>::iterator it=st.begin();
*/
// vector<vector<int>>edge; 二维虚拟储存坐标 
//-----------------------------------------------------------------------------------------------------------------*/
  //map<int,bool>mp[N]; 
ll num[30];
ll a[N];
bool vis[N];
ll l,r;
bool check(ll k){//暴力判断k是否可行
	ll num=0;
	for(ll i=l;i<=r;i++){
		vis[i]=0;
	}
	for(ll i=l;i<=r;i++){
		ll tmp=a[i]^k;
		if(tmp<=r&&tmp>=l&&!vis[tmp]){
			num++;
			vis[tmp]=1;
		}
		else{
			return 0;
		}
	}
	return num==r-l+1;
}
int main(){
/*cout<<setiosflags(ios::fixed)<<setprecision(8)<<ans<<endl;//输出ans(float)格式控制为8位小数(不含整数部分)*/
/*cout<<setprecision(8)<<ans<<endl;//输出ans(float)格式控制为8位小数(含整数部分)*/
	ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);//同步流
	ll t;
	cin>>t;
	while(t--){
		memset(num,0,sizeof(num));
		cin>>l>>r;
		vector<ll>c;
		for(ll i=l;i<=r;i++){
			cin>>a[i];
			c.push_back(a[i]^l);//将可能答案存到b中
			for(ll j=0;j<=20;j++){
				if(a[i]&(1ll<<j)){//记录每一位上1的个数
					num[j]++;
				}
			}
		}
		ll ans=0;
		if(check(ans)){
			cout<<ans<<endl;
			continue;
		}
		ll sum=(l+r)*(r-l+1ll)/2;
		for(ll i=c.size()-1;i>=0;i--){//
			ll res=0;
			ll tmp=c[i];
			for(ll j=0;j<=20;j++){
				if(c[i]&(1ll<<j)){//若当前二进制位上位1,那么异或a后a数组上位0的全取1
					res+=(ll)(r-l+1ll-num[j])*(1ll<<j);
				}
				else{//反正取1的全取1
					res+=num[j]*(1ll<<j);
				}
			}
			if(res==sum){
				if(check(c[i])){//暴力判断
					ans=c[i];
					break;
				}
			}
		}
		cout<<ans<<endl;
	}
	return 0;
}

 

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值