CF1290F Making Shapes——数位背包DP

CF1290F Making Shapes

题目描述

一位 怀好意的拉题人的翻译：

题解

前置知识：计算几何，背包DP，数位DP。

首先容易发现，如果你选的每个向量的系数 $c_i$（可能为0）确定了，那么把它们按极角排序后的凸包是唯一确定的，所以我们把问题转换为求不同的系数有多少组。

如果我们把两坐标按正负分开考虑，就会发现
$$\begin{gathered} \sum _{x_i>0}{c}_i{x}_i=\sum _{x_j<0}{c}_j(-x_j) \\ \sum _{y_i>0}{c}_i{y}_i=\sum _{y_j<0}{c}_j(-y_j) \end{gathered}$$这是因为最终所有线段必须围城闭合的凸多边形，所以向量和为0。而这个 $m$ 的限制呢？显然两坐标的极差就等于所有同号的 $x$ 或 $y$ 的和，也就是说
$$\sum _{x_i>0}{c}_i{x}_i\le m且\sum _{y_i>0}{c}_i{y}_i\le m$$这就是为什么上面的向量和为零要把正负分开写而不是直接求和。

显然，如果 $m$ 很小的话，这是可以做背包的。

但是这题 $m$ 很大，不能用常规的背包。那怎么办呢，这时一般要放弃背包而想想其它办法…

不，我们可以想办法优化背包！

一般背包的优化有哪些呢？最著名的莫过于倍增优化了。传统的倍增优化是先分开考虑每种物品，然后用倍增拆开它的数量。这里我们把它颠倒一下，先统一用倍增拆开向量的系数，拆成二进制的30位，然后再考虑每个向量的系数在每个二进制位上是否有值。
这时显然是会产生进位的，所以就要从低位到高位考虑。每次转移需要考虑 $2^n$ 种情况，进位最多有 $\sum x_i\le 4n$，但是这题 $n$ 和 $x_i$ 非常小，所以这么DP居然是可行的。

这就是更高级的倍增优化：考虑进位的数位背包DP！

具体地，我们需要记录7个状态（7维DP）：
$dp[w][a][b][c][d][p=0/1][q=0/1]$，表示
当前考虑第 $w$ 个二进制位，
正的 $x_i$ 产生的进位为 $a$，
负的 $x_i$ 产生的进位为 $b$，
正的 $y_i$ 产生的进位为 $c$，
负的 $y_i$ 产生的进位为 $d$，
$p$：${\sum }_{x_i>0}{c}_i{x}_i$ 的前 $w-1$ 位是否 $\le m$，
$q$：${\sum }_{y_i>0}{c}_i{y}_i$ 的前 $w-1$ 位是否 $\le m$。

转移的时候要保证每个二进制位上${\sum }_{x_i>0}{c}_i{x}_i={\sum }_{x_j<0}{c}_j(-x_j)$，${\sum }_{y_i>0}{c}_i{y}_i={\sum }_{y_j<0}{c}_j(-y_j)$，且保证最高位没有进位，$p$ 和 $q$ 都为1。

最后要减去系数全为0的1种情况，因为要求凸包面积非0。

总复杂度为 $O(2^n(4n{)}^4\log m)$，写成记忆化搜索会方便一些。

代码

#include<bits/stdc++.h>//JZM yyds!!
#define ll long long
#define uns unsigned
#define IF (it->first)
#define IS (it->second)
#define END putchar('\n')
using namespace std;
const int MAXN=-1;
const ll INF=1e18;
inline ll read(){
	ll x=0;bool f=1;char s=getchar();
	while((s<'0'||s>'9')&&s>0){if(s=='-')f^=1;s=getchar();}
	while(s>='0'&&s<='9')x=(x<<1)+(x<<3)+(s^48),s=getchar();
	return f?x:-x;
}
int ptf[30],lpt;
inline void print(ll x,char c='\n'){
	if(x<0)putchar('-'),x=-x;
	ptf[lpt=1]=x%10;
	while(x>9)x/=10,ptf[++lpt]=x%10;
	while(lpt)putchar(ptf[lpt--]^48);
	if(c>0)putchar(c);
}
inline ll lowbit(ll x){return x&-x;}

const ll MOD=998244353;
int n,k,x[10],y[10];
ll m,ans;
inline ll ads(ll x){return x>0?x:-x;}
inline void ad(int&a,int b){a+=b;if(a>=MOD)a-=MOD;}
int dp[32][22][22][22][22][2][2];
inline int DP(int w,int a,int b,int c,int d,bool p,bool q){
	if(w==30)return (!a&&!b&&!c&&!d&&p&&q);
	if(dp[w][a][b][c][d][p][q]>0)return dp[w][a][b][c][d][p][q]-1;
	int&res=dp[w][a][b][c][d][p][q],e=(m>>w)&1;
	for(int s=0;s<(1<<n);s++){
		int aa=a,bb=b,cc=c,dd=d;
		for(int i=1;i<=n;i++)if((s>>(i-1))&1){
			if(x[i]>=0)aa+=x[i];else bb-=x[i];
			if(y[i]>=0)cc+=y[i];else dd-=y[i];
		}if((aa&1)==(bb&1)&&(cc&1)==(dd&1))
			ad(res,DP(w+1,aa>>1,bb>>1,cc>>1,dd>>1,
			(aa&1)<e||((aa&1)==e&&p),(cc&1)<e||((cc&1)==e&&q)));
	}return res++;
}
signed main()
{
	n=read(),m=read();
	for(int i=1;i<=n;i++)x[i]=read(),y[i]=read();
	print((DP(0,0,0,0,0,1,1)-1+MOD)%MOD);
	return 0;
}