BZOJ 3105: [cqoi2013]新Nim游戏 [高斯消元XOR 线性基]

以后我也要用传送门!


题意:一些数,选择一个权值最大的异或和不为0的集合


 

终于有点明白线性基是什么了...等会再整理

求一个权值最大的线性无关子集

线性无关子集满足拟阵的性质,贪心选择权值最大的,用高斯消元判断是否和已选择的线性相关

每一位记录pivot[i]为i用到的行

枚举要加入的数字的每一个二进制为1的位,如果有pivot[i]那么就异或一下(消元),否则pivot[i]=这个数并退出

如果最后异或成0了就说明线性相关...

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <bitset>
using namespace std;
typedef long long ll;
const int N=105,INF=1e9;
inline int read(){
    char c=getchar();int x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}
int n,a[N],bin[35];
ll ans,sum;
void ini(){
    bin[0]=1;
    for(int i=1;i<=30;i++) bin[i]=bin[i-1]<<1;
}
int pivot[N];
void Gauss(){
    for(int i=1;i<=n;i++){
        int _=a[i];
        for(int j=30;j>=0;j--) if(a[i]&bin[j]){
            if(pivot[j]) a[i]^=a[pivot[j]];
            else {pivot[j]=i;break;}
        }
        if(a[i]) ans+=_;
    }
}
int main(){
    freopen("in","r",stdin);
    ini();
    n=read();
    for(int i=1;i<=n;i++) a[i]=read(),sum+=a[i];
    sort(a+1,a+1+n,greater<int>());
    Gauss();
    if(!ans) puts("-1");
    else printf("%lld",sum-ans);
}

 

bzoj2460 一样的

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <bitset>
using namespace std;
typedef long long ll;
const int N=1005;
inline ll read(){
    char c=getchar();ll x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}
int n,pivot[N];
ll ans,bin[N];
struct data{
    ll a,w;
    bool operator <(const data &r)const{return w>r.w;}
}a[N];
void ini(){
    bin[0]=1;
    for(int i=1;i<=60;i++) bin[i]=bin[i-1]<<1;
}
void Gauss(){
    for(int i=1;i<=n;i++){
        for(int j=60;j>=0;j--) if(a[i].a&bin[j]){
            if(pivot[j]) a[i].a^=a[pivot[j]].a;
            else {pivot[j]=i;break;}
        }
        if(a[i].a) ans+=a[i].w;
    }
}
int main(){
    freopen("in","r",stdin);
    ini();
    n=read();
    for(int i=1;i<=n;i++) a[i].a=read(),a[i].w=read();
    sort(a+1,a+1+n);
    Gauss();
    printf("%lld",ans);
}

 

转载于:https://www.cnblogs.com/candy99/p/6414237.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值