题目:
我们定义n位自幂数为:
如果x是n位数,x的各位的n次方的也为x,那么x为n位自幂数
求n位自幂数有哪些,0< n <=19,x > 0,排序后输出
时限3000ms。
解法1:打表
首先,如果拿到一个数列,可以先试试这个网站:OEIS,输入数列串,就会返回通项公式同时给出表,可以以后直接打表用。
解法2:搜索
怎么搜索??
暴力枚举n位?最差情况9^19种情况,会炸。时间复杂度链接
其实我们观察,174和741,其实是一种东西,所以求解自幂数的关键是求一个数x中它每一位数字重复出现的个数,至于这些数字怎么排列其实没关系,我们要做的就是枚举各个数字重复出现的个数即可。
代码如下:
#include<cstdio>
#include<string.h>
#include<algorithm>
using namespace std;
typedef long long ll;
int n, a[20], f[20], cnt;
ll ans[20];
//快速幂
ll poww(int a, int b) {
ll ans = 1, base = a;
while (b != 0) {
if (b & 1 != 0)
ans *= base;
base *= base;
b >>= 1;
}
return ans;
}
int calc(ll sum) {
ll t = sum;
int num = 0;
memset(f, 0, sizeof(f));
while (t) {
ll tmp = t % 10;
t /= 10;
f[tmp] ++;
num++;
}
return num;
}
//数字比较大,在爆ll的边缘试探,所以这里拿double检测一下
//ps是填充到哪个数字了,hs是还能填充几个数字
void dfs(int ps, int hs, ll sum, double dbsum) {
if (ps == 9) {
a[ps] = hs;
sum += (ll)hs*poww(9, n);
dbsum += (double)hs*poww(9, n);
if (dbsum <= 1e19&&sum > 0) {
if (calc(sum) == n) {
bool flag = 1;
for (int i = 0; i <= 9; i++) {
if (f[i] != a[i])flag = 0;
}
if (flag)ans[cnt++] = sum;
}
}
return;
}
for (int i = 0; i <= hs; i++) {
a[ps] = i;
dfs(ps + 1, hs - i, sum + (ll)i*poww(ps, n), dbsum + (double)i*poww(ps, n));
}
}
int main() {
scanf("%d", &n);
memset(ans, 0, sizeof(ans)); cnt = 0;
dfs(0, n, 0, 0);
if (!cnt) { printf("-1
"); }
else {
sort(ans, ans + cnt);
for (int i = 0; i < cnt; i++) {
printf("%lld ", ans[i]);
}
printf("
");
}
return 0;
}