Solution -「CF 1290F」Making Shapes

$\mathscr{Description}$

  给定平面向量集 $\newcommand{\vct}[1]{\boldsymbol{#1}}\{\vct v_n\}$，求从 $(0,0)$ 开始用这些向量作为方向画线段，求能画出的本质不同的封闭凸多边形，满足面积非 $0$，且在 $x$ 轴和 $y$ 轴上的投影长度都不超过 $m$。答案对 $998244353$ 取模。

  $n\le 5$，向量集任意分量在 $[-4,4]$，$m\le10^9$，此外所有数据都是整数。

$\mathscr{Solution}$

  出题人去为文惠君解牛吧！

  抽象一下要求，问题等价与求

$$\begin{cases} \sum_i c_ix_i=0\\ \sum_i c_iy_i=0\\ \sum_{x_i>0} c_ix_i\le m\\ \sum_{y_i>0} c_iy_i\le m \end{cases}$$

中 $\{c_n\}$ 的非负整数解数量。这什么玩意儿？类欧？半平面交？同余最短路？直接解离大谱了属于是。

  好，出题人开发了（或者搬运了）一个极其神奇的 trick：你不是要 $\sum_ic_ix_i=0$ 吗？由于 $x_i$ 很小，所以 每次加法对 bit 的改变有局限，那么，我们，数位 DP！

  由于正负混杂不方便，我们记 $\textit{pxs}$ 表示正的 $x$ 之和，$\textit{nxs}$ 表示负的 $x$ 之和；$\textit{pys}$ 和 $\textit{nys}$ 同理。从这个基本事实入手：确定了 $c$ 的所有最低位之后，四个值的奇偶性（最低位）是确定的。这意味这若现在就有 $\textit{pxs}$ 和 $\textit{nxs}$ 的最低位不同，那么不管 $c$ 怎么取，$\textit{pxs}\neq\textit{nxs}$，联系限制，这种情况可以舍弃。

  推而广之，确定了 $c$ 的所有低 $i$ 位之后，四个值的低 $i$ 位也都是确定的。这个时候，对于前两个限制，只需要记录 $\textit{pxs}$ 等值除掉低 $i$ 位的数值即可进行后续判断；后两个限制直接维护低 $i$ 位是否有 $\sum_{x_i>0}c_i^{(0..i-1)}x_i\le m^{(0..i-1)}$ 以及 $y$ 对应的式子即可。因此，我们得到这样一个状态：

  令 $f(i,\textit{pxs},\textit{nxs},\textit{pys},\textit{nys},0/1,0/1)$ 表示考虑了 $c$ 的低 $i$ 位，取到了这四个值，分别是否 $\le m^{(0..i-1)}$。转移时，只需要确定所有 $c$ 的第 $i+1$ 位，本质上是在取 $\{\vct v_n\}$ 的子集，所以预处理出每个子集对四个值的贡献，就能做到 $\mathcal O(\log m\cdot(nx)^4\cdot2\cdot2\cdot2^n)$ 转移。

$\mathscr{Code}$

/*+Rainybunny+*/

#include <bits/stdc++.h>

#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)

const int MAXN = 5, MAXX = 4, MAXS = 1 << MAXN, MOD = 998244353;
int n, m, x[MAXN], y[MAXN];
int px[MAXS], py[MAXS], nx[MAXS], ny[MAXS];
int f[2][MAXN * MAXX + 1][MAXN * MAXX + 1]
  [MAXN * MAXX + 1][MAXN * MAXX + 1][2][2];

inline void addeq(int& u, const int v) { (u += v) >= MOD && (u -= MOD); }

int main() {
    // freopen("shape.in", "r", stdin);
    // freopen("shape.out", "w", stdout);

    scanf("%d %d", &n, &m);
    rep (i, 0, n - 1) scanf("%d %d", &x[i], &y[i]);
    rep (S, 0, (1 << n) - 1) {
        rep (i, 0, n - 1) if (S >> i & 1) {
            if (x[i] > 0) px[S] += x[i];
            else nx[S] -= x[i];
            if (y[i] > 0) py[S] += y[i];
            else ny[S] -= y[i];
        }
    }

    f[0][0][0][0][0][0][0] = 1;
    int sta = 0;
    rep (i, 0, 31 - __builtin_clz(m)) {
        int mi = m >> i & 1;
        rep (spx, 0, n << 2) rep (snx, 0, n << 2)
        rep (spy, 0, n << 2) rep (sny, 0, n << 2) {
            rep (upx, 0, 1) rep (upy, 0, 1) {
                int& cur = f[sta][spx][snx][spy][sny][upx][upy];
                if (!cur) continue;
                rep (S, 0, (1 << n) - 1) {
                    int tpx = spx + px[S], tnx = snx + nx[S];
                    int tpy = spy + py[S], tny = sny + ny[S];
                    if ((tpx ^ tnx) & 1 || (tpy ^ tny) & 1) continue;
                    addeq(f[!sta][tpx >> 1][tnx >> 1][tpy >> 1][tny >> 1]
                      [(tpx & 1) != mi ? tpx & 1 : upx]
                      [(tpy & 1) != mi ? tpy & 1 : upy], cur);
                }
                cur = 0;
            }
        }
        sta ^= 1;
    }
    printf("%d\n", (f[sta][0][0][0][0][0][0] + MOD - 1) % MOD);
    return 0;
}