题意
给定一个数N,问所有满足下述性质的集合的个数:
集合中所有数的数码都只出现一次。
例如:{135,27,49}符合,{134,27,49}不符合(因为出现了两个4),{177,23}不符合(因为177出现了两次7)
题解
第一想法是:
f[S]
表示数码集为S的情况下,组出一个
≤N
的数,有多少种方案?
这一步是一个数位DP,枚举每一位统计即可,唯一需要注意的是,同一个数中数码也只能出现一次,统计的时候改一下即可。
第二步比赛时想偏了,先说说偏了的想法:
注意到现在我们只要依次枚举集合中每个数所取到的数码集,统计就能得到答案。
int oriS;
void dfs(int d, int lS, modn ans) {
for (int S = lS+1; S < (1<<10); ++ S) if (!(oriS&S) && f[S].n) {
oriS ^= S;
dfs(d+1, S, ans * f[S]);
oriS ^= S;
}
res = res + ans;
}
复杂度计算即为:10个不同的元素,放入m个相同的盒子中,m不定。。(可能是 ∑Stirling2(n,k) ),然而这是过不了的。。。
QwQ
有了上面“错误”的想法,接下来的想法也挺直观,即一个类似n!转为2^n的问题。
令
g[S][x]
表示目前数码集为
S
,所取数集合有
第二维主要目的是为了判重,比赛时因为没搞清怎么判重而没去想这个角度Orz,提供一个判重方法是:
思考一个3个数的全排列:
1 2 3
1 3 2
2 1 3
2 3 1
3 1 2
3 2 1
重新排序后:
1 2 3
2 1 3
1 3 2
3 1 2
2 3 1
3 2 1
可以发现,如果把最后一列遮住,那么前面组成了少一个元素的全排列。
由此得到启发,我们每次只加入一个元素的话,就可以保证其求排列的情况下不重复了
for (int SS = (S-1)&S; SS; SS = (SS-1)&S)
g[S][i] += g[S^SS][1] * g[SS][i-1];
那么,在统计答案的时候只需要除以 x <script type="math/tex" id="MathJax-Element-28">x</script>的阶乘即可
res = res + g[S][i] / fact[i];
看了标程以后,发现又更“正常”的写法。。。
我们可以发现,如果每次只加入包含最大的数 or 包含最小的数的集合,这样,就可以轻易地保证不会算重复。。。
包含最大的数:
for(int i = 1; i < 1<<10; ++i) {
int k = 0;
for(int j = 0; i >> j; j++) k = 1 << j;
f[i] = c[i];
for(int j = i ^ k; j; j = (j - 1) & i) {
f[i] = (f[i] + 1ll * f[j] * c[i^j] % Mod) % Mod;
}
ans = (ans + f[i]) % Mod;
}
包含最小的数:
for(int i = 1; i < 1<<10; ++i) {
int k = 0;
f[i] = 0;
for(int j = i; j; j = (j - 1) & i) if (j&(i&-i)) {
f[i] = (f[i] + 1ll * f[i^j] * c[j] % Mod) % Mod;
}
ans = (ans + f[i]) % Mod;
}
另外有个小写法: 清洗掉最后的“1”:S = S & (S-1) or S = S ^ (S & (-S))
code
#include <algorithm>
#include <bitset>
#include <cassert>
#include <cmath>
#include <cassert>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <ctime>
#include <deque>
#include <iostream>
#include <map>
#include <queue>
#include <set>
#include <string>
#include <vector>
using namespace std;
typedef int modn;
typedef long long LL;
modn f[1<<10], g[1<<10][10];
char num[20];
modn res;
int n, m;
modn fact[10], fact_[10];
int solve() {
scanf("%d", &n); ++ n;
sprintf(num, "%d", n); m = strlen(num);
if (n == 1) return 0;
for (int S = 0; S < (1<<10); ++ S) {
int s = 0;
for (int i = 0; i < 10; ++ i) if (S&(1<<i)) ++ s;
if (s > m) f[S] = 0;
else if (s < m) {
if (S&(1<<0))
f[S] = (LL)fact[s-1] * (s-1);
else
f[S] = fact[s];
} else {
f[S] = 0;
int _S = S;
for (int i = 0; i < m; ++ i) {
int j = m-1-i, p = num[i] - '0';
int x;
for (x = 0; x < p; ++ x) {
if (!(i == 0 && x == 0) && (_S&(1<<x))) {
f[S] = f[S] + fact[j];
}
}
if (!(_S&(1<<x))) break;
_S ^= (1<<x);
}
}
}
res = 0;
memset(g, 0, sizeof g);
for (int S = 1; S < (1<<10); ++ S) {
g[S][1] = f[S]; res += f[S];
for (int i = 2; i < 10; ++ i) {
for (int SS = (S-1)&S; SS; SS = (SS-1)&S)
g[S][i] += g[S^SS][1] * g[SS][i-1];
res = res + g[S][i] / fact[i];
}
}
printf("%d\n", res);
return 1;
}
int main() {
// freopen("E.in", "r", stdin);
fact[0] = 1;
for (int i = 1; i < 10; ++ i) fact[i] = fact[i-1] * i;
while (solve());
// for(;;);
return 0;
}