奶牛异或

农民约翰面临一个奶牛异或问题,需要找出奶牛序列中异或值最大的连续子序列。问题涉及排序和数值范围在0..2^21 - 1之间的奶牛。解题关键在于利用 Trie 树进行高效求解,找到最大异或值及其起始和结束位置。

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

题目描述
农民约翰在喂奶牛的时候被另一个问题卡住了。他的所有N(1 <= N <= 100,000)个奶牛在他面前排成一行(按序号1..N的顺序),按照它们的社会等级排序。奶牛#1由最高的社会等级,奶牛#N最低。每个奶牛同时被赋予了一个唯一的数在0..2^21 - 1的范围内。帮助农民约翰找出应该从那一头奶牛开始喂,使得从它开始的某一个连续的自序列上的奶牛的数的异或最大。如果有多个这样的子序列,选择结尾的奶牛社会等级最高的。如果还不唯一,选择最短的。

输入
• 第1行:一个单独的整数N。 • 第2到N + 1行:N个0..2^21 – 1之间的整数,代表每头奶牛的被赋予的数。第j行描述了社会等级j – 1的奶牛。

输出
• 第 1 行: 3个空格隔开的整数,分别为:最大的异或值,序列的起始位置、终止位置。

样例输入
5
1
0
5
4
2
样例输出
6 4 5

Solution

trie树的应用

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
int n,v,tot,ans,cow;
int a[100005],s[100005],b[40];
struct ty
{
    int Next[2];
    void init()
    {
        memset(Next,-1,sizeof(Next));
    }
}p[3000005];
void work(int x)
{
    int k=0;
    while(x>0) 
    {
        k++;
        b[k]=x%2;
        x=x/2;
    }
    for(int i=k+1;i<=31;i++) b[i]=0;
}
void Insert()
{
    int now,k;
    now=0;
    for(int i=31;i>=1;i--) 
    {
        k=b[i];
        if(p[now].Next[k]==-1) 
        {
            tot++;
            p[tot].init();
            p[now].Next[k]=tot;
        }
        now=p[now].Next[k];
    }
}
void query(int id)
{
    int now,k,val;
    now=0;val=0;
    for(int i=31;i>=1;i--) 
    {
        k=b[i];
        if(p[now].Next[k^1]!=-1) 
        {
            val+=(1<<(i-1));
            now=p[now].Next[k^1];
        }
        else now=p[now].Next[k];
    }
    if(val>ans) 
    {
        ans=val;
        cow=id;
    }
    if(val==ans) 
    {
        if(id<cow) cow=id;
    }
}
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++) 
    {
        scanf("%d",&a[i]);
        s[i]=s[i-1]^a[i];
    }
    ans=-1e9;
    p[0].init();
    work(0);
    Insert();
    for(int i=1;i<=n;i++) 
    {
        work(s[i]);
        Insert();
        query(i);
    }
    for(int i=cow;i>=1;i--) 
    {
        v=v^a[i];
        if(v==ans) 
        {
            cout<<ans<<' '<<i<<' '<<cow;
            return 0;
        }
    }
    return 0;
}
#include<cstdio> #include<queue> #define N 500 #define M 200000 #define INF 0x7fffffff using namespace std; int n,f,d,head[N],to[M],c[M],nxt[M],num[N],minn,s,t,cnt,minflow,maxflow,tot; //用邻接表存边 queue<int>q; void adde(int u,int v,int a){//建边 cnt++; to[cnt]=v; c[cnt]=a; nxt[cnt]=head[u]; head[u]=cnt; cnt++; to[cnt]=u; c[cnt]=0; nxt[cnt]=head[v]; head[v]=cnt; } bool add_num() {//给每个点标号 while(!q.empty()) { q.pop(); } for(int i=s;i<=t+n;i++) { num[i]=-1; } num[s]=1; q.push(s); while(!q.empty()) { int now=q.front(); q.pop(); for(int i=head[now];i;i=nxt[i]) { if(c[i]&&num[to[i]]==-1) { num[to[i]]=num[now]+1; q.push(to[i]); } } } if(num[t]==-1) {//如果num[t]==-1,那么就说明到不了汇点,不能给汇点标号。 return false; } else { return true; } } int dfs(int u,int s){ if(u==t){ return s; } for(int i=head[u];i;i=nxt[i]){ if(c[i]&&num[u]+1==num[to[i]]&&(minflow=dfs(to[i],min(s,c[i])))){ c[i]-=minflow; c[i^1]+=minflow; return minflow; } } return 0; } void dinic(){//dinic模板 while(add_num()){ while((minn=dfs(1,INF))){ maxflow+=minn; } } } int main(){ scanf(“%d%d%d”,&n,&f,&d); cnt=1; s=1; t=1+f+n+d+1; for(int i=1;i<=f;i++){//超源连每种食品 adde(s,1+i,1); } for(int i=1;i<=d;i++){//每种饮料连超汇 adde(1+f+n+i,t,1); } for(int i=1;i<=n;i++){//入点连出点 adde(1+f+i,1+f+n+d+1+i,1); } for(int i=1;i<=n;i++){ int dn,fn; scanf(“%d%d”,&fn,&dn); for(int q=1;q<=fn;q++){//食品连喜欢吃它的奶牛 int fi; scanf(“%d”,&fi); adde(1+fi,1+f+i,1); } for(int q=1;q<=dn;q++){//奶牛连他喜欢吃的饮料 int di; scanf(“%d”,&di); adde(1+f+n+d+1+i,1+f+n+di,1); } } dinic(); printf(“%d\n”,maxflow);//输出 return 0; } 将这段代码转换为邻接矩阵实现
最新发布
03-08
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值