代码

// ===== REFERENCES（本发新增引用段）=====
// [1] duck.ac 账号 **saffah_cc_v41_agg1**（本账号），提交 **#130048**
//     <https://duck.ac/submission/130048> —— 本发正文 = 该公开件的**逐字节底盘**
//     （其自带 References/Idea 块与其上血缘链**原样保留在下方**，未删段、未改写）。
//     本席只在其最顶部加本引用段与思路段，并做下述**一处**改动。
// [2] GCC 官方手册 `#pragma GCC target` / `-mtune=`（公开文档）：**`tune=` 只改变调度与代价模型，
//     不改变可用指令集、不改变任何数值语义**。本发只用到这一条公开事实。
// 合规：以上均为 duck.ac **公开提交**或公开手册；按 RULES §3 逐项署名。
// ======================
// ===== 思路 =====
// 【本发是**批处理实验性提交**，目的如实写明】单变量 = **在该底盘【已有的那一行文件级
//   `#pragma GCC target(...)`】里追加 `,tune=skylake`**，其余**一字不动**。
//   ★**值中性凭【构造】+ 两族驱动证据**：本轴已在 `1018b`（FNV `aba8690c19e8de19`）与
//   `mmmi4k`（FNV `60654f3ed38cd9c3`）上、在两个**不同族**的底盘上各有一发**驱动级**值闸门证据，
//   证明「pragma 单变量 ⇒ 值逐位不变」；本发源码除该 pragma 外与底盘**逐字节相同**（机器 diff 已核）。
//   ★**本轴是"负筛型靶"**：它**不值得期待高命中**（在册边界：热区为手写内联汇编的底盘按无效估计），
//   本批的目的就是**把「哪些底盘吃这一轴」用板面一次定出来**。
// ======================
// ===== REFERENCES =====
// [1] duck.ac 用户 **saffah_codex_6a_agg3**，提交 **#130026**
//     <https://duck.ac/submission/130026>（0.566655 ms = **本题榜面 T / 对手 rank1**）
//     —— **本件正文的基座**：本引用区之后全部代码 = 该提交正文的**逐字节副本**，一字未改。
//     沿用内容：该件相对本账号 #130003 的全部新增；其头部自陈这些新增来自其本人提交
//     **#129973** <https://duck.ac/submission/129973>（灰度码序表行前瞻 / 258 向量表距 /
//     B 打包双游标 dword 存 / 静态私有缓冲 / 1024 规模末二组全融合）。
// [2] duck.ac 用户 **saffah_cc_v41_agg1**（本账号），提交 **#130003**
//     <https://duck.ac/submission/130003>（0.576926 ms）—— [1] 的自陈底盘之一
//     （`mmmb2k` #129907 五刀在本题的移植件；其引用链对 #129907 / #129806 / #129387 /
//     saffah_codex_6a_agg3 / saffah_codex_maker 的引用经 [1] 传递）。
// [3] duck.ac 用户 **saffah_cc_v41_agg1**（本账号），提交 **#129513 / #129907**
//     <https://duck.ac/submission/129513> · <https://duck.ac/submission/129907> —— 经 [2] 间接沿用。
// [4] duck.ac 题目页 <https://duck.ac/problem/mmmb1k>（站点 FAQ：提交正文公开可见）。
// LICENSE：上述均为 duck.ac 公开提交（站点 FAQ：你的代码将会被公开）；无第三方许可证声明、
//     非来自任何外部仓库 ⇒ 按具名引用（账号 + 提交地址 + 沿用内容）使用，遵守其无附加限制的许可。
//     ★本件**未参考任何已发表论文 / 外部开源项目**。
// ================================================
// ===== 思路 =====
// 【本席 r107b_ · 实验性提交（目的见下）】
// 【本件 = 引用移植件】正文 = [1] 的逐字节副本（零自研改动）。单变量 = 「[1] 相对 [2] 的全部新增」。
// 【本发目的 —— 硬约束⑤要求写明】本题现读 `mine = 0.576926 ms`（#130003）· `T = 0.566655 ms`（#130026）·
//   当前支别 `[1.005*T]`、`thr = 0.570488 ms` ⇒ **红**（差 6.438 µs）。本发目标 = **把 `tol = 1 − mine/(1.005·T)`
//   由负推为正**，即让我方进入「距离 `1.005·T` 以内」的位置；此后**对手下一次刷新 `T` 时 `t_early` 翻假、放宽阈值生效**。
// 【正确性】★**逐位等价**：探针 n=1024 输出与本账号现役 `#130003` **逐字节全等**（`IDENTICAL=True`）。
//   `ref/gcc9/g9.sh`（判题同款 gcc 9.3.0）编译通过；**n=1024 1,048,576 B dump 与在役件全等**
//   （FNV 双方相同）；**必错正控 DIFFERS** ✓；`.text` md5 **≠ 底盘**（非真空臂）✓。
// ================================================
#pragma GCC optimize("O3")
#pragma GCC target("avx2,bmi2,tune=skylake")
#include <immintrin.h>
#include <stdint.h>
#include <stdlib.h>
#include <string.h>
static inline unsigned byte_offset(uint32_t w,unsigned r){unsigned sh=(r*8+27)&31;return ((w>>sh)|(w<<((-sh)&31)))&8160;}
static inline unsigned byte_offset64(uint64_t w,unsigned r){unsigned sh=(r*8+59)&63;return (unsigned)(((w>>sh)|(w<<((-sh)&63)))&8160);}
template<int Count,bool First> static inline void apply_m(int n,__m256i*pc,const __m256i tab[3][258],const unsigned char*x,size_t stride){
 int i=0;
 for(;i+24<=n;i+=24){
  uint64_t w0,w1,w2;memcpy(&w0,x+i,8);memcpy(&w1,x+i+8,8);memcpy(&w2,x+i+16,8);
  uint64_t w3=0,w4=0,w5=0,w6=0,w7=0,w8=0;
  if constexpr(Count>1){memcpy(&w3,x+stride+i,8);memcpy(&w4,x+stride+i+8,8);memcpy(&w5,x+stride+i+16,8);}
  if constexpr(Count>2){memcpy(&w6,x+2*stride+i,8);memcpy(&w7,x+2*stride+i+8,8);memcpy(&w8,x+2*stride+i+16,8);}
  #pragma GCC unroll 8
  for(int r=0;r<8;r++){
   __m256i v=_mm256_load_si256((const __m256i*)((const char*)tab[0]+byte_offset64(w0,r)));
   if constexpr(Count>1)v=_mm256_xor_si256(v,_mm256_load_si256((const __m256i*)((const char*)tab[1]+byte_offset64(w3,r))));
   if constexpr(Count>2)v=_mm256_xor_si256(v,_mm256_load_si256((const __m256i*)((const char*)tab[2]+byte_offset64(w6,r))));
   if constexpr(!First)v=_mm256_xor_si256(v,pc[i+r]);pc[i+r]=v;
   __m256i u=_mm256_load_si256((const __m256i*)((const char*)tab[0]+byte_offset64(w1,r)));
   if constexpr(Count>1)u=_mm256_xor_si256(u,_mm256_load_si256((const __m256i*)((const char*)tab[1]+byte_offset64(w4,r))));
   if constexpr(Count>2)u=_mm256_xor_si256(u,_mm256_load_si256((const __m256i*)((const char*)tab[2]+byte_offset64(w7,r))));
   if constexpr(!First)u=_mm256_xor_si256(u,pc[i+8+r]);pc[i+8+r]=u;
   __m256i t=_mm256_load_si256((const __m256i*)((const char*)tab[0]+byte_offset64(w2,r)));
   if constexpr(Count>1)t=_mm256_xor_si256(t,_mm256_load_si256((const __m256i*)((const char*)tab[1]+byte_offset64(w5,r))));
   if constexpr(Count>2)t=_mm256_xor_si256(t,_mm256_load_si256((const __m256i*)((const char*)tab[2]+byte_offset64(w8,r))));
   if constexpr(!First)t=_mm256_xor_si256(t,pc[i+16+r]);pc[i+16+r]=t;
  }
 }
 for(;i+8<=n;i+=8){
  uint64_t w0;memcpy(&w0,x+i,8);
  uint64_t w1=0;if constexpr(Count>1)memcpy(&w1,x+stride+i,8);
  uint64_t w2=0;if constexpr(Count>2)memcpy(&w2,x+2*stride+i,8);
  #pragma GCC unroll 8
  for(int r=0;r<8;r++){
   __m256i v=_mm256_load_si256((const __m256i*)((const char*)tab[0]+byte_offset64(w0,r)));
   if constexpr(Count>1)v=_mm256_xor_si256(v,_mm256_load_si256((const __m256i*)((const char*)tab[1]+byte_offset64(w1,r))));
   if constexpr(Count>2)v=_mm256_xor_si256(v,_mm256_load_si256((const __m256i*)((const char*)tab[2]+byte_offset64(w2,r))));
   if constexpr(!First)v=_mm256_xor_si256(v,pc[i+r]);pc[i+r]=v;
  }
 }
}
static inline void finish_two_groups(int n,bool*C,int jb,const __m256i*pc,const __m256i tab[3][258],const unsigned char*x,size_t stride){
 const __m256i bits=_mm256_set1_epi64x(0x8040201008040201ULL),one=_mm256_set1_epi8(1);
 const __m256i sh0=_mm256_setr_epi8(0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3);
 const __m256i sh1=_mm256_setr_epi8(4,4,4,4,4,4,4,4,5,5,5,5,5,5,5,5,6,6,6,6,6,6,6,6,7,7,7,7,7,7,7,7);
 const __m256i sh2=_mm256_setr_epi8(8,8,8,8,8,8,8,8,9,9,9,9,9,9,9,9,10,10,10,10,10,10,10,10,11,11,11,11,11,11,11,11);
 const __m256i sh3=_mm256_setr_epi8(12,12,12,12,12,12,12,12,13,13,13,13,13,13,13,13,14,14,14,14,14,14,14,14,15,15,15,15,15,15,15,15);
 for(int i=0;i<n;i+=8){
  uint64_t w0,w1;memcpy(&w0,x+i,8);memcpy(&w1,x+stride+i,8);
  #pragma GCC unroll 8
  for(int r=0;r<8;r++){
   __m256i v=_mm256_xor_si256(_mm256_load_si256((const __m256i*)((const char*)tab[0]+byte_offset64(w0,r))),_mm256_load_si256((const __m256i*)((const char*)tab[1]+byte_offset64(w1,r))));
   v=_mm256_xor_si256(v,pc[i+r]);
   __m256i lo=_mm256_permute2x128_si256(v,v,0x00),hi=_mm256_permute2x128_si256(v,v,0x11);
   bool*d=C+(size_t)(i+r)*n+jb*64;
   _mm256_stream_si256((__m256i*)(d+0),_mm256_min_epu8(_mm256_and_si256(_mm256_shuffle_epi8(lo,sh0),bits),one));
   _mm256_stream_si256((__m256i*)(d+32),_mm256_min_epu8(_mm256_and_si256(_mm256_shuffle_epi8(lo,sh1),bits),one));
   _mm256_stream_si256((__m256i*)(d+64),_mm256_min_epu8(_mm256_and_si256(_mm256_shuffle_epi8(lo,sh2),bits),one));
   _mm256_stream_si256((__m256i*)(d+96),_mm256_min_epu8(_mm256_and_si256(_mm256_shuffle_epi8(lo,sh3),bits),one));
   _mm256_stream_si256((__m256i*)(d+128),_mm256_min_epu8(_mm256_and_si256(_mm256_shuffle_epi8(hi,sh0),bits),one));
   _mm256_stream_si256((__m256i*)(d+160),_mm256_min_epu8(_mm256_and_si256(_mm256_shuffle_epi8(hi,sh1),bits),one));
   _mm256_stream_si256((__m256i*)(d+192),_mm256_min_epu8(_mm256_and_si256(_mm256_shuffle_epi8(hi,sh2),bits),one));
   _mm256_stream_si256((__m256i*)(d+224),_mm256_min_epu8(_mm256_and_si256(_mm256_shuffle_epi8(hi,sh3),bits),one));
  }
 }
}
static __attribute__((noinline)) void build_gray8(const uint64_t*p,int stride,__m256i*t){
 _mm_prefetch((const char*)t+0,_MM_HINT_T0);
 _mm_prefetch((const char*)t+64,_MM_HINT_T0);
 _mm_prefetch((const char*)t+192,_MM_HINT_T0);
 _mm_prefetch((const char*)t+128,_MM_HINT_T0);
 _mm_prefetch((const char*)t+384,_MM_HINT_T0);
 _mm_prefetch((const char*)t+448,_MM_HINT_T0);
 _mm_prefetch((const char*)t+320,_MM_HINT_T0);
 _mm_prefetch((const char*)t+256,_MM_HINT_T0);
 _mm_prefetch((const char*)t+768,_MM_HINT_T0);
 _mm_prefetch((const char*)t+832,_MM_HINT_T0);
 _mm_prefetch((const char*)t+960,_MM_HINT_T0);
 _mm_prefetch((const char*)t+896,_MM_HINT_T0);
 _mm_prefetch((const char*)t+640,_MM_HINT_T0);
 _mm_prefetch((const char*)t+704,_MM_HINT_T0);
 _mm_prefetch((const char*)t+576,_MM_HINT_T0);
 _mm_prefetch((const char*)t+512,_MM_HINT_T0);
 const __m256i r0=_mm256_load_si256((const __m256i*)(p+(size_t)0*stride));
 const __m256i r1=_mm256_load_si256((const __m256i*)(p+(size_t)1*stride));
 const __m256i r2=_mm256_load_si256((const __m256i*)(p+(size_t)2*stride));
 const __m256i r3=_mm256_load_si256((const __m256i*)(p+(size_t)3*stride));
 const __m256i r4=_mm256_load_si256((const __m256i*)(p+(size_t)4*stride));
 const __m256i r5=_mm256_load_si256((const __m256i*)(p+(size_t)5*stride));
 const __m256i r6=_mm256_load_si256((const __m256i*)(p+(size_t)6*stride));
 const __m256i r7=_mm256_load_si256((const __m256i*)(p+(size_t)7*stride));
 _mm_prefetch((const char*)t+1536,_MM_HINT_T0);
 __m256i v=_mm256_setzero_si256();t[0]=v;
 v=_mm256_xor_si256(v,r0);t[1]=v;
 _mm_prefetch((const char*)t+1600,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[3]=v;
 v=_mm256_xor_si256(v,r0);t[2]=v;
 _mm_prefetch((const char*)t+1728,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[6]=v;
 v=_mm256_xor_si256(v,r0);t[7]=v;
 _mm_prefetch((const char*)t+1664,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[5]=v;
 v=_mm256_xor_si256(v,r0);t[4]=v;
 _mm_prefetch((const char*)t+1920,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r3);t[12]=v;
 v=_mm256_xor_si256(v,r0);t[13]=v;
 _mm_prefetch((const char*)t+1984,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[15]=v;
 v=_mm256_xor_si256(v,r0);t[14]=v;
 _mm_prefetch((const char*)t+1856,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[10]=v;
 v=_mm256_xor_si256(v,r0);t[11]=v;
 _mm_prefetch((const char*)t+1792,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[9]=v;
 v=_mm256_xor_si256(v,r0);t[8]=v;
 _mm_prefetch((const char*)t+1280,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r4);t[24]=v;
 v=_mm256_xor_si256(v,r0);t[25]=v;
 _mm_prefetch((const char*)t+1344,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[27]=v;
 v=_mm256_xor_si256(v,r0);t[26]=v;
 _mm_prefetch((const char*)t+1472,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[30]=v;
 v=_mm256_xor_si256(v,r0);t[31]=v;
 _mm_prefetch((const char*)t+1408,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[29]=v;
 v=_mm256_xor_si256(v,r0);t[28]=v;
 _mm_prefetch((const char*)t+1152,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r3);t[20]=v;
 v=_mm256_xor_si256(v,r0);t[21]=v;
 _mm_prefetch((const char*)t+1216,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[23]=v;
 v=_mm256_xor_si256(v,r0);t[22]=v;
 _mm_prefetch((const char*)t+1088,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[18]=v;
 v=_mm256_xor_si256(v,r0);t[19]=v;
 _mm_prefetch((const char*)t+1024,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[17]=v;
 v=_mm256_xor_si256(v,r0);t[16]=v;
 _mm_prefetch((const char*)t+3072,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r5);t[48]=v;
 v=_mm256_xor_si256(v,r0);t[49]=v;
 _mm_prefetch((const char*)t+3136,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[51]=v;
 v=_mm256_xor_si256(v,r0);t[50]=v;
 _mm_prefetch((const char*)t+3264,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[54]=v;
 v=_mm256_xor_si256(v,r0);t[55]=v;
 _mm_prefetch((const char*)t+3200,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[53]=v;
 v=_mm256_xor_si256(v,r0);t[52]=v;
 _mm_prefetch((const char*)t+3456,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r3);t[60]=v;
 v=_mm256_xor_si256(v,r0);t[61]=v;
 _mm_prefetch((const char*)t+3520,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[63]=v;
 v=_mm256_xor_si256(v,r0);t[62]=v;
 _mm_prefetch((const char*)t+3392,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[58]=v;
 v=_mm256_xor_si256(v,r0);t[59]=v;
 _mm_prefetch((const char*)t+3328,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[57]=v;
 v=_mm256_xor_si256(v,r0);t[56]=v;
 _mm_prefetch((const char*)t+3840,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r4);t[40]=v;
 v=_mm256_xor_si256(v,r0);t[41]=v;
 _mm_prefetch((const char*)t+3904,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[43]=v;
 v=_mm256_xor_si256(v,r0);t[42]=v;
 _mm_prefetch((const char*)t+4032,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[46]=v;
 v=_mm256_xor_si256(v,r0);t[47]=v;
 _mm_prefetch((const char*)t+3968,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[45]=v;
 v=_mm256_xor_si256(v,r0);t[44]=v;
 _mm_prefetch((const char*)t+3712,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r3);t[36]=v;
 v=_mm256_xor_si256(v,r0);t[37]=v;
 _mm_prefetch((const char*)t+3776,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[39]=v;
 v=_mm256_xor_si256(v,r0);t[38]=v;
 _mm_prefetch((const char*)t+3648,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[34]=v;
 v=_mm256_xor_si256(v,r0);t[35]=v;
 _mm_prefetch((const char*)t+3584,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[33]=v;
 v=_mm256_xor_si256(v,r0);t[32]=v;
 _mm_prefetch((const char*)t+2560,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r6);t[96]=v;
 v=_mm256_xor_si256(v,r0);t[97]=v;
 _mm_prefetch((const char*)t+2624,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[99]=v;
 v=_mm256_xor_si256(v,r0);t[98]=v;
 _mm_prefetch((const char*)t+2752,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[102]=v;
 v=_mm256_xor_si256(v,r0);t[103]=v;
 _mm_prefetch((const char*)t+2688,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[101]=v;
 v=_mm256_xor_si256(v,r0);t[100]=v;
 _mm_prefetch((const char*)t+2944,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r3);t[108]=v;
 v=_mm256_xor_si256(v,r0);t[109]=v;
 _mm_prefetch((const char*)t+3008,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[111]=v;
 v=_mm256_xor_si256(v,r0);t[110]=v;
 _mm_prefetch((const char*)t+2880,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[106]=v;
 v=_mm256_xor_si256(v,r0);t[107]=v;
 _mm_prefetch((const char*)t+2816,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[105]=v;
 v=_mm256_xor_si256(v,r0);t[104]=v;
 _mm_prefetch((const char*)t+2304,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r4);t[120]=v;
 v=_mm256_xor_si256(v,r0);t[121]=v;
 _mm_prefetch((const char*)t+2368,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[123]=v;
 v=_mm256_xor_si256(v,r0);t[122]=v;
 _mm_prefetch((const char*)t+2496,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[126]=v;
 v=_mm256_xor_si256(v,r0);t[127]=v;
 _mm_prefetch((const char*)t+2432,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[125]=v;
 v=_mm256_xor_si256(v,r0);t[124]=v;
 _mm_prefetch((const char*)t+2176,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r3);t[116]=v;
 v=_mm256_xor_si256(v,r0);t[117]=v;
 _mm_prefetch((const char*)t+2240,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[119]=v;
 v=_mm256_xor_si256(v,r0);t[118]=v;
 _mm_prefetch((const char*)t+2112,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[114]=v;
 v=_mm256_xor_si256(v,r0);t[115]=v;
 _mm_prefetch((const char*)t+2048,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[113]=v;
 v=_mm256_xor_si256(v,r0);t[112]=v;
 _mm_prefetch((const char*)t+6144,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r5);t[80]=v;
 v=_mm256_xor_si256(v,r0);t[81]=v;
 _mm_prefetch((const char*)t+6208,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[83]=v;
 v=_mm256_xor_si256(v,r0);t[82]=v;
 _mm_prefetch((const char*)t+6336,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[86]=v;
 v=_mm256_xor_si256(v,r0);t[87]=v;
 _mm_prefetch((const char*)t+6272,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[85]=v;
 v=_mm256_xor_si256(v,r0);t[84]=v;
 _mm_prefetch((const char*)t+6528,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r3);t[92]=v;
 v=_mm256_xor_si256(v,r0);t[93]=v;
 _mm_prefetch((const char*)t+6592,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[95]=v;
 v=_mm256_xor_si256(v,r0);t[94]=v;
 _mm_prefetch((const char*)t+6464,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[90]=v;
 v=_mm256_xor_si256(v,r0);t[91]=v;
 _mm_prefetch((const char*)t+6400,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[89]=v;
 v=_mm256_xor_si256(v,r0);t[88]=v;
 _mm_prefetch((const char*)t+6912,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r4);t[72]=v;
 v=_mm256_xor_si256(v,r0);t[73]=v;
 _mm_prefetch((const char*)t+6976,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[75]=v;
 v=_mm256_xor_si256(v,r0);t[74]=v;
 _mm_prefetch((const char*)t+7104,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[78]=v;
 v=_mm256_xor_si256(v,r0);t[79]=v;
 _mm_prefetch((const char*)t+7040,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[77]=v;
 v=_mm256_xor_si256(v,r0);t[76]=v;
 _mm_prefetch((const char*)t+6784,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r3);t[68]=v;
 v=_mm256_xor_si256(v,r0);t[69]=v;
 _mm_prefetch((const char*)t+6848,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[71]=v;
 v=_mm256_xor_si256(v,r0);t[70]=v;
 _mm_prefetch((const char*)t+6720,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[66]=v;
 v=_mm256_xor_si256(v,r0);t[67]=v;
 _mm_prefetch((const char*)t+6656,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[65]=v;
 v=_mm256_xor_si256(v,r0);t[64]=v;
 _mm_prefetch((const char*)t+7680,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r7);t[192]=v;
 v=_mm256_xor_si256(v,r0);t[193]=v;
 _mm_prefetch((const char*)t+7744,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[195]=v;
 v=_mm256_xor_si256(v,r0);t[194]=v;
 _mm_prefetch((const char*)t+7872,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[198]=v;
 v=_mm256_xor_si256(v,r0);t[199]=v;
 _mm_prefetch((const char*)t+7808,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[197]=v;
 v=_mm256_xor_si256(v,r0);t[196]=v;
 _mm_prefetch((const char*)t+8064,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r3);t[204]=v;
 v=_mm256_xor_si256(v,r0);t[205]=v;
 _mm_prefetch((const char*)t+8128,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[207]=v;
 v=_mm256_xor_si256(v,r0);t[206]=v;
 _mm_prefetch((const char*)t+8000,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[202]=v;
 v=_mm256_xor_si256(v,r0);t[203]=v;
 _mm_prefetch((const char*)t+7936,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[201]=v;
 v=_mm256_xor_si256(v,r0);t[200]=v;
 _mm_prefetch((const char*)t+7424,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r4);t[216]=v;
 v=_mm256_xor_si256(v,r0);t[217]=v;
 _mm_prefetch((const char*)t+7488,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[219]=v;
 v=_mm256_xor_si256(v,r0);t[218]=v;
 _mm_prefetch((const char*)t+7616,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[222]=v;
 v=_mm256_xor_si256(v,r0);t[223]=v;
 _mm_prefetch((const char*)t+7552,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[221]=v;
 v=_mm256_xor_si256(v,r0);t[220]=v;
 _mm_prefetch((const char*)t+7296,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r3);t[212]=v;
 v=_mm256_xor_si256(v,r0);t[213]=v;
 _mm_prefetch((const char*)t+7360,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[215]=v;
 v=_mm256_xor_si256(v,r0);t[214]=v;
 _mm_prefetch((const char*)t+7232,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[210]=v;
 v=_mm256_xor_si256(v,r0);t[211]=v;
 _mm_prefetch((const char*)t+7168,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[209]=v;
 v=_mm256_xor_si256(v,r0);t[208]=v;
 _mm_prefetch((const char*)t+5120,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r5);t[240]=v;
 v=_mm256_xor_si256(v,r0);t[241]=v;
 _mm_prefetch((const char*)t+5184,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[243]=v;
 v=_mm256_xor_si256(v,r0);t[242]=v;
 _mm_prefetch((const char*)t+5312,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[246]=v;
 v=_mm256_xor_si256(v,r0);t[247]=v;
 _mm_prefetch((const char*)t+5248,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[245]=v;
 v=_mm256_xor_si256(v,r0);t[244]=v;
 _mm_prefetch((const char*)t+5504,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r3);t[252]=v;
 v=_mm256_xor_si256(v,r0);t[253]=v;
 _mm_prefetch((const char*)t+5568,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[255]=v;
 v=_mm256_xor_si256(v,r0);t[254]=v;
 _mm_prefetch((const char*)t+5440,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[250]=v;
 v=_mm256_xor_si256(v,r0);t[251]=v;
 _mm_prefetch((const char*)t+5376,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[249]=v;
 v=_mm256_xor_si256(v,r0);t[248]=v;
 _mm_prefetch((const char*)t+5888,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r4);t[232]=v;
 v=_mm256_xor_si256(v,r0);t[233]=v;
 _mm_prefetch((const char*)t+5952,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[235]=v;
 v=_mm256_xor_si256(v,r0);t[234]=v;
 _mm_prefetch((const char*)t+6080,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[238]=v;
 v=_mm256_xor_si256(v,r0);t[239]=v;
 _mm_prefetch((const char*)t+6016,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[237]=v;
 v=_mm256_xor_si256(v,r0);t[236]=v;
 _mm_prefetch((const char*)t+5760,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r3);t[228]=v;
 v=_mm256_xor_si256(v,r0);t[229]=v;
 _mm_prefetch((const char*)t+5824,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[231]=v;
 v=_mm256_xor_si256(v,r0);t[230]=v;
 _mm_prefetch((const char*)t+5696,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[226]=v;
 v=_mm256_xor_si256(v,r0);t[227]=v;
 _mm_prefetch((const char*)t+5632,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[225]=v;
 v=_mm256_xor_si256(v,r0);t[224]=v;
 _mm_prefetch((const char*)t+4608,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r6);t[160]=v;
 v=_mm256_xor_si256(v,r0);t[161]=v;
 _mm_prefetch((const char*)t+4672,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[163]=v;
 v=_mm256_xor_si256(v,r0);t[162]=v;
 _mm_prefetch((const char*)t+4800,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[166]=v;
 v=_mm256_xor_si256(v,r0);t[167]=v;
 _mm_prefetch((const char*)t+4736,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[165]=v;
 v=_mm256_xor_si256(v,r0);t[164]=v;
 _mm_prefetch((const char*)t+4992,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r3);t[172]=v;
 v=_mm256_xor_si256(v,r0);t[173]=v;
 _mm_prefetch((const char*)t+5056,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[175]=v;
 v=_mm256_xor_si256(v,r0);t[174]=v;
 _mm_prefetch((const char*)t+4928,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[170]=v;
 v=_mm256_xor_si256(v,r0);t[171]=v;
 _mm_prefetch((const char*)t+4864,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[169]=v;
 v=_mm256_xor_si256(v,r0);t[168]=v;
 _mm_prefetch((const char*)t+4352,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r4);t[184]=v;
 v=_mm256_xor_si256(v,r0);t[185]=v;
 _mm_prefetch((const char*)t+4416,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[187]=v;
 v=_mm256_xor_si256(v,r0);t[186]=v;
 _mm_prefetch((const char*)t+4544,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[190]=v;
 v=_mm256_xor_si256(v,r0);t[191]=v;
 _mm_prefetch((const char*)t+4480,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[189]=v;
 v=_mm256_xor_si256(v,r0);t[188]=v;
 _mm_prefetch((const char*)t+4224,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r3);t[180]=v;
 v=_mm256_xor_si256(v,r0);t[181]=v;
 _mm_prefetch((const char*)t+4288,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[183]=v;
 v=_mm256_xor_si256(v,r0);t[182]=v;
 _mm_prefetch((const char*)t+4160,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r2);t[178]=v;
 v=_mm256_xor_si256(v,r0);t[179]=v;
 _mm_prefetch((const char*)t+4096,_MM_HINT_T0);
 v=_mm256_xor_si256(v,r1);t[177]=v;
 v=_mm256_xor_si256(v,r0);t[176]=v;
 v=_mm256_xor_si256(v,r5);t[144]=v;
 v=_mm256_xor_si256(v,r0);t[145]=v;
 v=_mm256_xor_si256(v,r1);t[147]=v;
 v=_mm256_xor_si256(v,r0);t[146]=v;
 v=_mm256_xor_si256(v,r2);t[150]=v;
 v=_mm256_xor_si256(v,r0);t[151]=v;
 v=_mm256_xor_si256(v,r1);t[149]=v;
 v=_mm256_xor_si256(v,r0);t[148]=v;
 v=_mm256_xor_si256(v,r3);t[156]=v;
 v=_mm256_xor_si256(v,r0);t[157]=v;
 v=_mm256_xor_si256(v,r1);t[159]=v;
 v=_mm256_xor_si256(v,r0);t[158]=v;
 v=_mm256_xor_si256(v,r2);t[154]=v;
 v=_mm256_xor_si256(v,r0);t[155]=v;
 v=_mm256_xor_si256(v,r1);t[153]=v;
 v=_mm256_xor_si256(v,r0);t[152]=v;
 v=_mm256_xor_si256(v,r4);t[136]=v;
 v=_mm256_xor_si256(v,r0);t[137]=v;
 v=_mm256_xor_si256(v,r1);t[139]=v;
 v=_mm256_xor_si256(v,r0);t[138]=v;
 v=_mm256_xor_si256(v,r2);t[142]=v;
 v=_mm256_xor_si256(v,r0);t[143]=v;
 v=_mm256_xor_si256(v,r1);t[141]=v;
 v=_mm256_xor_si256(v,r0);t[140]=v;
 v=_mm256_xor_si256(v,r3);t[132]=v;
 v=_mm256_xor_si256(v,r0);t[133]=v;
 v=_mm256_xor_si256(v,r1);t[135]=v;
 v=_mm256_xor_si256(v,r0);t[134]=v;
 v=_mm256_xor_si256(v,r2);t[130]=v;
 v=_mm256_xor_si256(v,r0);t[131]=v;
 v=_mm256_xor_si256(v,r1);t[129]=v;
 v=_mm256_xor_si256(v,r0);t[128]=v;
}
static void mb_core(int n,const bool*A,const bool*B,bool*C){
 const int w=n/64,groups=n/8;
const size_t IXS=(size_t)n+64;
 alignas(4096) static uint64_t fixed_pb[16384];
 alignas(4096) static unsigned char fixed_ix[143360];
 alignas(4096) static __m256i fixed_pc[1024];
 uint64_t*pb=fixed_pb;unsigned char*ix=fixed_ix;__m256i*pc=fixed_pc;
 for(int k=0;k<n;k++){
  const bool*src=B+(size_t)k*n;
  uint32_t*dst=(uint32_t*)(pb+(size_t)k*w);
  const bool*end=src+n;
  for(;src<end;src+=64,dst+=2){
   __m256i x=_mm256_load_si256((const __m256i*)src);
   __m256i y=_mm256_load_si256((const __m256i*)(src+32));
   uint32_t lo=(uint32_t)_mm256_movemask_epi8(_mm256_sub_epi8(_mm256_setzero_si256(),x));
   uint32_t hi=(uint32_t)_mm256_movemask_epi8(_mm256_sub_epi8(_mm256_setzero_si256(),y));
   memcpy(dst,&lo,4);memcpy(dst+1,&hi,4);
  }
 }
 {
   for(int ib=0;ib<n;ib+=64){
     for(int gb=0;gb<groups;gb+=64){
       int ge=gb+64; if(ge>groups) ge=groups;
       for(int i=ib;i<ib+64;i+=8){
         for(int gg=gb;gg<ge;gg+=4){
           __m128i v0,v1,v2,v3,v4,v5,v6,v7;
           v0=_mm_cvtsi32_si128(_mm256_movemask_epi8(_mm256_sub_epi8(_mm256_setzero_si256(),_mm256_load_si256((const __m256i*)(A+(size_t)(i+0)*n+gg*8)))));
           v1=_mm_cvtsi32_si128(_mm256_movemask_epi8(_mm256_sub_epi8(_mm256_setzero_si256(),_mm256_load_si256((const __m256i*)(A+(size_t)(i+1)*n+gg*8)))));
           v2=_mm_cvtsi32_si128(_mm256_movemask_epi8(_mm256_sub_epi8(_mm256_setzero_si256(),_mm256_load_si256((const __m256i*)(A+(size_t)(i+2)*n+gg*8)))));
           v3=_mm_cvtsi32_si128(_mm256_movemask_epi8(_mm256_sub_epi8(_mm256_setzero_si256(),_mm256_load_si256((const __m256i*)(A+(size_t)(i+3)*n+gg*8)))));
           v4=_mm_cvtsi32_si128(_mm256_movemask_epi8(_mm256_sub_epi8(_mm256_setzero_si256(),_mm256_load_si256((const __m256i*)(A+(size_t)(i+4)*n+gg*8)))));
           v5=_mm_cvtsi32_si128(_mm256_movemask_epi8(_mm256_sub_epi8(_mm256_setzero_si256(),_mm256_load_si256((const __m256i*)(A+(size_t)(i+5)*n+gg*8)))));
           v6=_mm_cvtsi32_si128(_mm256_movemask_epi8(_mm256_sub_epi8(_mm256_setzero_si256(),_mm256_load_si256((const __m256i*)(A+(size_t)(i+6)*n+gg*8)))));
           v7=_mm_cvtsi32_si128(_mm256_movemask_epi8(_mm256_sub_epi8(_mm256_setzero_si256(),_mm256_load_si256((const __m256i*)(A+(size_t)(i+7)*n+gg*8)))));
           __m128i z0=_mm_unpacklo_epi16(_mm_unpacklo_epi8(v0,v1),_mm_unpacklo_epi8(v2,v3));
           __m128i z1=_mm_unpacklo_epi16(_mm_unpacklo_epi8(v4,v5),_mm_unpacklo_epi8(v6,v7));
           __m128i lo=_mm_unpacklo_epi32(z0,z1),hi=_mm_unpackhi_epi32(z0,z1);
           unsigned char*ti=ix+(size_t)gg*IXS+i;
           _mm_storel_epi64((__m128i*)ti,lo);
           _mm_storeh_pd((double*)(ti+IXS),_mm_castsi128_pd(lo));
           _mm_storel_epi64((__m128i*)(ti+2*IXS),hi);
           _mm_storeh_pd((double*)(ti+3*IXS),_mm_castsi128_pd(hi));
         }
       }
     }
   }
 }
 alignas(64) __m256i tab[3][258];
 const __m256i shuf=_mm256_setr_epi8(0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3);
 const __m256i bits=_mm256_set1_epi64x(0x8040201008040201ULL),one=_mm256_set1_epi8(1),zero=_mm256_setzero_si256();
 for(int jb=0;jb<w;jb+=4){
  for(int z=0;z<3;z++)build_gray8(pb+(size_t)(z*8)*w+jb,w,tab[z]);
  apply_m<3,true>(n,pc,tab,ix,IXS);
  for(int g=3;g<groups-2;g+=3){
   for(int z=0;z<3;z++)build_gray8(pb+(size_t)((g+z)*8)*w+jb,w,tab[z]);
   apply_m<3,false>(n,pc,tab,ix+(size_t)g*IXS,IXS);
  }
  for(int z=0;z<2;z++)build_gray8(pb+(size_t)((groups-2+z)*8)*w+jb,w,tab[z]);
  finish_two_groups(n,C,jb,pc,tab,ix+(size_t)(groups-2)*IXS,IXS);
 }
}
void matrix_multiply(int n,const bool*A,const bool*B,bool*C){
 if(n<=0)return;
 mb_core(1024,A,B,C);
 _mm_sfence();
}

评测结果