hdu 3949 XOR

这里写链接内容

分析:高斯消元的板子题。其实我也不知道算不算。。
大致就是先把这些数字全部变成线性基,把每个数变成二进制形态高斯消元消就行了。需要注意的是在这里需要把矩阵变成最简的形式,为下一步考虑。
试想现在把 N 个数变成了t个数的线性基,把二进制表示的01矩阵换成上三角矩阵。那么对于第 i 个线性基a[i],从第 i+1 到第 t 个的异或最大值还是无法超过a[i]。这个于是现在只要相当于将 K 分解为二进制表示,二进制为1的位i就将 a[i] 异或进去。但是这样的前提条件就是线性基组成的矩阵需要是最简的,给个例子吧。 [1011] 这个矩阵并不是最简的形态,那么对于第三个值,并不能表示成 a[0]a[1] ,但如果是最简的 [1001] ,这个表示方法就没有问题了。
还需要注意的是如果矩阵秩小于矩阵大小,那么意味着最小的异或值可以为0,可以组成的方案数为 2t 而非 2t1 相应的还有一些改动。

/*****************************************************/
//#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <map>
#include <set>
#include <ctime>
#include <stack>
#include <queue>
#include <cmath>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <sstream>
#include <cstdlib>
#include <iostream>
#include <algorithm>
using namespace std;
#define   offcin        ios::sync_with_stdio(false)
#define   sigma_size    26
#define   lson          l,m,v<<1
#define   rson          m+1,r,v<<1|1
#define   slch          v<<1
#define   srch          v<<1|1
#define   sgetmid       int m = (l+r)>>1
#define   ll            long long
#define   ull           unsigned long long
#define   mem(x,v)      memset(x,v,sizeof(x))
#define   lowbit(x)     (x&-x)
#define   bits(a)       __builtin_popcount(a)
#define   mk            make_pair
#define   pb            push_back
#define   fi            first
#define   se            second

const int    INF    = 0x3f3f3f3f;
const ll     INFF   = 1e18;
const double pi     = acos(-1.0);
const double inf    = 1e18;
const double eps    = 1e-9;
const ll     mod    = 1e9+7;
const int    maxmat = 10;
const ull    BASE   = 133333331;

/*****************************************************/
inline void RI(int &x) {
      char c;
      while((c=getchar())<'0' || c>'9');
      x=c-'0';
      while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0';
 }
/*****************************************************/

const int maxm = 10005;

ll a[maxm];
int N, M;

int Gauss() {
    int k = 1;
    for (int i = 63; i >= 0; i --) {
        int t = 0;
        for (int j = k; j <= N; j ++) {
            if ((a[j] >> i) & 1) {
                t = j;
                break;
            }
        }
        if (t) {
            swap(a[t], a[k]);
            for (int j = 1; j <= N; j ++) //从1开始 化成最简矩阵形式
                if (j != k && (a[j] >> i) & 1) a[j] ^= a[k];
            k ++;
        }
    }
    return k - 1; 
}

int main(int argc, char const *argv[]) {
    int T; cin>>T;
    for (int kase = 1; kase <= T; kase ++) {
        printf("Case #%d:\n", kase);
        cin>>N;
        for (int i = 1; i <= N; i ++) scanf("%I64d", a + i);
        int cnt = Gauss();
        // for (int i = 1; i <= N; i ++) cout<<a[i]<<endl;
        ll num = (1ll << cnt);
        if (cnt == N) num --;  //如果矩阵秩==大小,那么方案数需要减一
        cin>>M;
        while (M --) {
            ll tmp; scanf("%I64d", &tmp);
            if (tmp > num) puts("-1");
            else {
                if (cnt < N) tmp --; //如果最小可以为0,那么相应后面的下标要-1
                ll ans = 0;
                for (int i = 0; i < cnt; i ++) {
                    if ((tmp >> i) & 1) ans ^= a[cnt - i];
                }
                printf("%I64d\n", ans);
            }
        }
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值