Xor Sum(Trie)

problem

Zeus 和 Prometheus 做了一个游戏,Prometheus 给 Zeus 一个集合,集合中包含了N个正整数,随后 Prometheus 将向 Zeus 发起M次询问,每次询问中包含一个正整数 S ,之后 Zeus 需要在集合当中找出一个正整数 K ,使得 K 与 S 的异或结果最大。Prometheus 为了让 Zeus 看到人类的伟大,随即同意 Zeus 可以向人类求助。你能证明人类的智慧么?

Input

输入包含若干组测试数据,每组测试数据包含若干行。
输入的第一行是一个整数T(T < 10),表示共有T组数据。
每组数据的第一行输入两个正整数N,M(<1=N,M<=100000),接下来一行,包含N个正整数,代表 Zeus 的获得的集合,之后M行,每行一个正整数S,代表 Prometheus 询问的正整数。所有正整数均不超过2^32。

Output

对于每组数据,首先需要输出单独一行”Case #?:”,其中问号处应填入当前的数据组数,组数从1开始计算。
对于每个询问,输出一个正整数K,使得K与S异或值最大。

Sample Input

2
3 2
3 4 5
1
5
4 1
4 6 5 6
3

Sample Output

Case #1:
4
3
Case #2:
4

思路

题目给定n个数 和m次询问

求每次询问和n个数中哪个数异或值最大

每次询问暴力找每个数的异或值复杂度是 O(mn) O ( m ∗ n )

现在模拟下这个过程看能否优化

当给我一个数tp时,我要从n中找一个数使得n^tp最大

那肯定看tp的二进制位,从高到低,我希望的串是tp二进制串取反后的串

而且对于二进制串,有这样一个性质就是

高位远强于低位,具体表现就是,高位的1个1比低位所有1之和都大

那现在关键就是寻找这样异或后更高位为1的串

自然想到前缀(高位)

就是说可能有一些数高位是一样的,且都是我想要的(和tp异或后为1)

但是在某一位呢,大家出现了分歧

那肯定选择我想要的那条路径继续看下去

那赤裸裸的是Trie

最坏时间复杂度 O(32mT) O ( 32 ∗ m ∗ T )

注意delete 释放内存(不知道为什么只小了10%的空间) ….

代码示例

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
const int charset=2;

struct Node{
    Node *next[charset];
    int flag;
    Node(){
        for(int i=0;i<charset;++i) next[i]=NULL;
        flag=0;
    }
};

Node *root;

void init()
{
    root=new Node();
}

void ins(char *s){
    int len=strlen(s);
    Node *now =root;
    for(int i=0;i<len;++i){
        int to=s[i]-'0';
        if(now->next[to]==NULL) now->next[to]=new Node();
        now=now->next[to];
    }
    now->flag=1;
}

char temp[40];

ll Find(char *s)
{
    int len=strlen(s);
    Node *now=root;
    for(int i=0;i<len;++i){
        s[i]=((s[i]-'0')^1)+'0';
    }
    ll k=(1LL<<34);
    ll ans=0;
    for(int i=0;i<len;++i){
        int to=s[i]-'0';
        if(now->next[to]==NULL){
            now=now->next[to^1];
        }
        else{
            ans+=k;
            now=now->next[to];
        }
        k>>=1;
    }
    return ans;
}

void del(Node *p)//free
{
    for(int i=0;i<charset;++i){
        if(p->next[i]!=NULL) del(p->next[i]);
    }
    free(p);
}

void trans(ll n)
{
    memset(temp,0,sizeof(temp));
    //35够用  高位补0
    int k=0;
    while(n)
    {
        temp[k++]='0'+(n&1);
        n>>=1;
    }
    for(int i=k;i<35;++i){
        temp[i]='0';
    }
    reverse(temp,temp+35);
}

int main()
{
    //freopen("in.txt","r",stdin);
    ios::sync_with_stdio(false);
    int t,tt;
    cin>>t;
    tt=0;
    while(t--)
    {
        init();
        cout<<"Case #"<<++tt<<":\n";
        int n,m;
        cin>>n>>m;
        ll tp;
        while(n--)
        {
            cin>> tp;
            trans(tp);
            ins(temp);
        }
        while(m--){
            cin>>tp;
            trans(tp);
            cout<<(Find(temp)^tp)<<endl;
            //cout<<strlen(temp)<<endl;
        }
        del(root);
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值