题面
题解
正解做法是线性基。
两个函数用同一个随机数种子进行随机,可以得到随机化的一样的结果,这是关键。
于是我们随机生成
m
m
m 个
n
n
n 位 01 串,程序 encoder
输出的串的每一位就表示是否异或这个 01 串,最终异或结果就得到原来的长度为
n
n
n 的串。
为了得到原串,我们要安排可以自由选择的那 m − k m-k m−k 个位置是否异或。用线性基就好。
不过,由于线性基上的一个串其实是原集合中的多个串异或的结果,所以要建一棵表示依赖关系的树之类的,最后通过线性基上的异或情况倒推出 m − k m-k m−k 个串的异或情况。
时间复杂度 O ( n 3 w ) O(\frac{n^3}{w}) O(wn3) 。
CODE
#include "password.h"
#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<stack>
#include<random>
#include<bitset>
#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
namespace {
#define MAXN 100005
#define LL long long
#define ULL unsigned long long
#define ENDL putchar('\n')
#define DB double
#define lowbit(x) (-(x) & (x))
#define FI first
#define SE second
#define BI bitset<1005>
BI bt[1005],ch[2100];
int ad[2100];
bool f[2100],v[2100];
vector<int> pr[2100];
void ins(BI x,int n,int id) {
for(int i = n-1;i >= 0;i --) {
if(!x[i]) continue;
if(f[i]) {
pr[id].push_back(ad[i]);
x ^= bt[i];
}
else {
bt[i] = x;
f[i] = 1; ad[i] = id;
break;
}
}return ;
}
}
void encoder(int n, int m, int k, const char* a, const char* b, char* ans) {
mt19937 g(124312);
for(int i = 0;i < m;i ++) {
for(int j = 0;j < n;j ++) ch[i][j] = g()&1;
}
BI me;
for(int i = 0;i < n;i ++) {
me[i] = a[i]-'0';
}
for(int i = 0;i < m;i ++) {
if(b[i] == '?') {
ins(ch[i],n,i);
}
else {
if(b[i] == '1') {
me ^= ch[i];
}
v[i] = b[i]-'0';
}
}
for(int i = n-1;i >= 0;i --) {
if(me[i]) v[ad[i]] = 1,me ^= bt[i];
}
for(int i = m-1;i >= 0;i --) {
if(v[i] && b[i] == '?') {
for(int j = 0;j < (int)pr[i].size();j ++) {
v[pr[i][j]] ^= 1;
}
}
ans[i] = '0'+v[i];
}
return ;
}
void decoder(int n, int m, const char* a, char* ans) {
mt19937 g(124312);
for(int i = 0;i < m;i ++) {
for(int j = 0;j < n;j ++) ch[i][j] = g()&1;
}
BI me;
for(int i = 0;i < m;i ++) {
if(a[i] == '1') me ^= ch[i];
}
for(int i = 0;i < n;i ++) ans[i] = me[i]+'0';
return ;
}