Solution -「JOISC 2021」古老的机器
\(\mathcal{Description}\)
Link.
这是一道通信题。
对于长度为一个 \(n\)，仅包含字符 X, Y, Z 的字符串 \(s\)，将其中 \(n\) 个字符按任意顺序删去，定义删除方案的权值为在子串 XYZ 中删除 Y 的次数。实现两个函数：
void Anna(int N, std::vector<char> S)，获取字符串信息，传递不超过 \(7\times10^4\) 个01位用于通信；void Bruno(int N, int L, std::vector<int> A)，获取通信信息，给出权值最大的删除方案。
\(n\le10^5\)。
\(\mathcal{Solution}\)
先看看怎么求最大权值。考虑当前最左侧 X 的位置 \(i\) 与最左侧 Z 的位置 \(j\)，若 \(j<i\)，显然删掉 \(j\) 无任何影响；否则考虑 \(i<k<j\)，\(s_k\) 必然取 X 或 Y。那么从 \(j-1\) 逆序删除到 \(i+1\)，最后删除掉 \(j\)，继续迭代。容易（真的容易）看出这就是最大化权值的方案。
那么 Anna 需要告诉 Bruno 哪些信息呢？——唯一的一个 \(i\)，和若干 \(j\)。一个小小的压缩方式是，对于连续的 Z，仅保留最右端的。特别地，在标记 X 的 \(1\) 的后面强制补一个 \(0\) 站位，我们就把信息串转化为长度为 \(n+1\)，不存在连续 \(1\) 的数字串。每 \(W\) 为一段，利用 Fibonacci 数列求出每段的字典序编号，再转化为二进制输出，就能传递信息啦。
观察标称可知，取 \(W=63\)，此时一段的总情况数接近 \(2^{44}\)，损失较小。运算过程当然是 \(\mathcal O(n)\) 的，信息长度也就比 \(7\times10^4\) 少一百来次。
\(\mathcal{Code}\)
Anna.cpp
/* Clearink */
#include "Anna.h"
#include <vector>
#ifndef MY_REP_DEFINED
#define MY_REP_DEFINED
#define rep( i, l, r ) for ( int i = l, rep##i = r; i <= rep##i; ++i )
#define per( i, r, l ) for ( int i = r, per##i = l; i >= per##i; --i )
#endif
namespace {
typedef unsigned long long ULL;
const int MAXN = 1e5, W = 63, B = 44;
int n;
std::vector<char> s;
bool key[MAXN + W + 5];
ULL fib[W + 5];
inline void mark() {
int i = 0;
while ( i < n && s[i] != 'X' ) ++i;
key[i] = true;
for ( int j = i + 1; j < n; ++j ) {
while ( j < n && s[j] != 'Z' ) ++j;
while ( j + 1 < n && s[j + 1] == 'Z' ) ++j;
if ( i >= n || j >= n ) break;
key[j + 1] = true;
}
}
inline void encode() {
fib[0] = 1, fib[1] = 2;
rep ( i, 2, W ) fib[i] = fib[i - 1] + fib[i - 2];
for ( int l = 0; l <= n; l += W ) {
ULL msg = 0;
rep ( i, l, l + W - 1 ) if ( key[i] ) {
msg += fib[repi - i];
}
rep ( i, 0, B - 1 ) Send( msg >> i & 1 );
}
}
} // namespace.
void Anna( int n, std::vector<char> s ) {
::n = n, ::s = s;
mark(), encode();
}
Bruno.cpp
/* Clearink */
#include "Bruno.h"
#include <vector>
#ifndef MY_REP_DEFINED
#define MY_REP_DEFINED
#define rep( i, l, r ) for ( int i = l, rep##i = r; i <= rep##i; ++i )
#define per( i, r, l ) for ( int i = r, per##i = l; i >= per##i; --i )
#endif
namespace {
typedef long long ULL;
const int MAXN = 1e5, W = 63, B = 44;
ULL fib[W + 5];
int n;
std::vector<int> msg;
bool key[MAXN + W + 5];
inline void decode() {
fib[0] = 1, fib[1] = 2;
rep ( i, 2, W ) fib[i] = fib[i - 1] + fib[i - 2];
for ( int l = 0; l < msg.size(); l += B ) {
ULL val = 0;
per ( i, l + B - 1, l ) val = val << 1 | msg[i];
rep ( i, l / B * W, i + W - 1 ) {
if ( val >= fib[repi - i] ) {
key[i] = true, val -= fib[repi - i];
}
}
}
}
inline void solve() {
int i = 0;
while ( i < n && !key[i] ) Remove( i++ );
if ( i == n ) return ;
for ( int j = i + 1, las = i; j < n; ) {
while ( j < n && !key[j + 1] ) ++j;
per ( k, j - 1, las + 1 ) Remove( k );
if ( j < n ) Remove( las = j++ );
}
Remove( i );
}
} // namespace.
void Bruno( int n, int l, std::vector<int> msg ) {
::n = n, ::msg = msg;
decode(), solve();
}
