[JZOJ3362] 【NOI2013模拟】数数

题目

题目大意

求区间$[A,B]$有多少个数是“完美的”。
一个数是“完美的”，当且仅当这个数的各位能分成两个集合，使得两个集合中数字的和相等。
$B\le 10^9$

---

思考历程

其实思考这题的时候已经没有什么时间了。
但我还是努力地去想正解。
看到的第一眼就会想到数位DP，然后脑中就弹出了一个DP状态。
很快就被自己证伪了，因为显然有重复。
想其它做法，却一直搞不出来。
看时间真的不多了，于是就去打暴力。
思想特别简单，枚举$[A,B]$中的数，然后判断它们是否可行。
最容易想到的是背包，但是我认为背包会爆炸，于是我打了折半搜索。
然后又搞了各种优化……
最终成功水到了50分。

---

水法

大把人是打表过的，我真是醉了……
50分的$10^7$数据范围其实有另一种更简单的方法：还是背包，但由于背包的值为$0$或$1$，范围就在$[0,45]$，就可以用一个long long将DP数组压起来，转移的时候直接位运算搞定……
YMQ说，数位DP一般都会开到$10^{18}$，这题只开到$10^9$，不是可以打表吗？
接下来介绍一下分段打表大法：
将$[1,10^{6}x]$（其中$1\le x\le 10^3$）的答案全部跑出来。
这样，在询问中一大段的区间就不用处理了，只需要暴力处理剩下的$10^6$个数……
时间还快得很……

---

正解

这题用传统的数位DP似乎不太好做，尤其是判断重复，更加繁琐。
于是我们就试着换一个角度，逃避重复的问题。
可以枚举组成数字的各位数分别有哪些。如果这些数能组成“完美的数”，就进行排列，将范围之内的数统计如答案。
这样显然就不用怕重复了，因为如果是平常的做法，各位数可能有多种组成和相等的两个集合的方法。但是这样不需要管有多少中方法，只要它是成功的就对了。
有重复元素的排列有个公式：$\frac{(\sum nu{m}_i)!}{\prod nu{m}_i!}$
不会证……NOIP初赛的资料上的……
枚举出来的东西不会太多，具体是多少，其实我也没有统计……郑姐说是$40000$到$80000$。
判断它们是否可以组成“完美的数”就用上面说的那个二进制优化背包的方法。
然后统计答案的时候像正常的数位DP一样卡着它的上限，如果当前位选的数小于边界，后面的就可以随便选，用公式算就好了。
所以这题还算是个数位DP吧……这是数位DP的特征啊……

---

代码

using namespace std;
#include <cstdio>
#include <cstring>
#include <algorithm>
int pow10[10],fac[10];
int low,up,w;
int lim;
int num[10];
int ans;
inline bool ok(){
	unsigned long long f=1;
	int s=0;
	for (int i=0;i<=9;++i){
		for (int j=1;j<=num[i];++j)
			f|=f<<i;
		s+=i*num[i];
	}
	return ((s&1^1) && (f>>(s>>1)&1));
}
inline int calc(int m){
	int fm=1;
	for (int i=0;i<=9;++i)
		fm*=fac[num[i]];
	return fac[m]/fm;
}
inline void dfs(int p){
	int th=lim/pow10[p]%10;
	for (int i=0;i<th;++i)
		if (num[i]){
			num[i]--;
			ans+=calc(p);
			num[i]++;
		}
	if (num[th] && p){
		num[th]--;
		dfs(p-1);
		num[th]++;
	}
}
void find(int k,int s){
	if (s==w){
		if (ok())
			dfs(w-1);
		return;
	}
	if (k>9)
		return;
	for (int i=0;i<=w-s;++i){
		num[k]=i;
		find(k+1,s+i);
	}
	num[k]=0;
}
int get_ans(int x){
	lim=++x;
	w=0;
	for (int tmp=lim;tmp;tmp/=10,++w);
	ans=0;
	find(0,0);
	return ans;
}
int main(){
	pow10[0]=1;
	for (int i=1;i<=9;++i)
		pow10[i]=pow10[i-1]*10;
	fac[0]=1;
	for (int i=1;i<=9;++i)
		fac[i]=fac[i-1]*i;
	scanf("%d%d",&low,&up);
	printf("%d\n",get_ans(up)-get_ans(low-1));
	return 0;
}

---

总结

在遇到会出现重复的问题是，可以试着换个角度，避免重复。