2019浙江理工大学校赛 D LIS（线性基）

最新推荐文章于 2022-04-20 11:44:44 发布

纸箱xyr

最新推荐文章于 2022-04-20 11:44:44 发布

阅读量510

点赞数

分类专栏：数论相关 ACM

本文链接：https://blog.csdn.net/qq_40619297/article/details/88768667

版权

ACM 同时被 2 个专栏收录

24 篇文章 0 订阅

订阅专栏

数论相关

3 篇文章 0 订阅

订阅专栏

题目链接

题目描述

令LIS(S)为序列S的最长递增子序列的长度
给你n个非负整数，a[0],a[1],...,a[n-1]，你可以对这个数组进行零次或多次操作，每次操作选择一个i(0<=i<=n-2),将a[i+1]变成(a[i+1] xor a[i])
你的任务是使得LIS(a)越大越好，输出LIS(a)的最大值

输入

第一行输入一个整数n (1 <= n <= 100)
第二行输入n个整数ai(0 <= ai <= 10^18)

输出

输出一个整数，即LIS(a)的最大值

样例输入

【样例输入1】
3
4 2 1
【样例输入2】
10
0 0 0 0 0 0 0 0 0 0
【样例输入3】
10
1 2 3 4 5 6 7 8 9 10
【样例输入4】
15 	
1008 42 7 5 2 9 75 0 0 12 3 6 81 4 3

样例输出

【样例输出1】
3
【样例输出2】
1
【样例输出3】
10
【样例输出4】
12

稍微推一下易知，变化后的 a_i 是前面i-1个数中任意个、任意位置与最初的 a_i 亦或的结果。

很容易想到线性基
假如当前我们正在确定第i个元素的贡献，我们就把前i个元素加入线性基。
然后查询当前元素排在第j位的最小值是，与原来在这个位置的值比较取min
- 若在开头，即使当前当前元素的最小值是什么
- 若在中间，既是求刚好大于dp[i-1][j-1]的值
  - 这个就像取第k大或者取最大那样二进制拆分一下就好了
访问一下dp数组最长达是多少

#include<algorithm>
#include<vector>
#include<iostream>
#include<math.h>
#include<cstring>
#include<string>
#include<stack>
#include<map>
#include<set>
#include<unordered_map>
#include<queue>
#include<assert.h>
#include<iomanip>
#include<bitset>
#include<stdio.h>

#define qcin; ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define pb push_back
#define mp make_pair
#define clr(x) memset(x,0,sizeof x)
#define fmax(x) memset(x,0x3f,sizeof x)
#define finit(x) memset(x,-1,sizeof x)
#define iio(n,m) io(n),io(m)
#define lc(p) (p<<1)
#define rc(p) ((p<<1)|1)
#define VI vector<int>

#define dis(l,r) r-l+1
#define gstr(str) scanf("%s",str)
#define glen(str) strlen(str)
using namespace std;

typedef long long ll;

typedef pair<int,int>pii;
const int maxn = 150;
const int mod =  1e9+7;
const ll INF = 2147483647;

typedef ll arr[maxn];
typedef char str[maxn];
void file(int x){if(x&&fopen("123.in","r")){freopen("123.in","r",stdin);}}

const long double pi = acos(-1);
const double eps=1e-10;
const double delta=0.993;

arr b,a;
ll dp[2][maxn];
int p=0,q=1;
ll res,n;

inline void ins(ll x){
    for(int i=62;~i;i--){
        if(x>>i&1){
            if(!b[i]){
                b[i]=x;
                for(int j=i-1;~j;j--)if(b[j]&&(b[i]>>j&1)){
                    b[i]^=b[j];
                }
                for(int j=62;j>i;j--)if((b[j]>>i)&1){
                    b[j]^=b[i];
                }
                return;
            }else x^=b[i];
        }
    }
}
inline ll getmin(ll x){
    for(int i=62;~i;i--){
        x=min(x,x^b[i]);
    }
    return x;
}
inline ll getnxt(int k,ll pre){
    ll x=a[k];
    for(int i=62;~i;--i){
        x=max(x,x^b[i]);
    }
    if(x<=pre)return dp[0][0];
    else{
        for(int i=62;~i;--i){
            ll tmp=x^b[i];
            if(tmp>x)continue;
            if(tmp>pre)x=min(x,tmp);
        }
        return x;
    }
}

int main(){
    file(1);
    while(~scanf("%lld",&n)){
        clr(b);fmax(dp[p]);
        for(int i=1;i<=n;i++){
            scanf("%lld",&a[i]);
        }
        dp[p][1]=a[1];
        for(int i=2;i<=n;i++){
            ins(a[i-1]);
            swap(p,q);fmax(dp[p]);
            dp[p][1]=min(dp[q][1],getmin(a[i]));
            for(int j=2;j<=i;j++){
                if(dp[q][j-1]==dp[0][0])break;
                dp[p][j]=min(dp[q][j],getnxt(i,dp[q][j-1]));
            }
        }
        res=1;
        for(int i=1;i<=n;i++){
            if(dp[p][i]!=dp[p][0])res=i;
        }
        printf("%lld\n",res);
    }
}