[LOJ#2328]「清华集训 2017」避难所
试题描述
“B君啊，你当年的伙伴都不在北京了，为什么你还在北京呢？”
“大概是因为出了一些事故吧，否则这道题就不叫避难所了。”
“唔，那你之后会去哪呢？”
“去一个没有冬天的地方。”
对于一个正整数 (n)，我们定义他在 (b) 进制下，各个位上的数的乘积为 (p = F(n, b))。
比如 (F(3338, 10) = 216)。
考虑这样一个问题，已知 (p) 和 (b)，求最小的 (n) 满足 (p = F(n, b))。
这是一个非常有趣的问题，对于一些 (b) 来说，我们可以贪心来做，比如如果 (b=10, p=216)。
我们可以从 (b-1) 到 (2) 试除，直到 (p) 为 (1) 为止，答案是 (389)，可以验证 (389) 是满足 (p = F(n, b)) 最小的 (n)。
但是对于一些进制 (b)，是不能用贪心做的，比如 (b = 9, p = 216)。使用贪心得到的解是 (3338)，而最优解是 (666)。(均为 (9) 进制下的。)
本题便是在给定进制 (b) 的情况下，举出一个这样的反例，或指出这样的反例不存在。
由于计算资源所限，反例中所有数字的乘积不能超过 (10^{18})。如果最小的反例中所有数字的乘积超过了 (10^{18})，那么也应该输出 (-1)。
输入
从标准输入读入数据。
第一行一个整数 (t)，表示一共有 (t) 组数据。
接下来每行一个整数 (b)，表示进制。
输出
输出到标准输出。
如果不存在反例，输出一行一个整数 (-1)。
如果存在反例，首先输出一个整数 (k)，表示反例 (n) 的位数，接下来在同一行输出 (k) 个十进制整数，表示任意一个反例的最优解。
输入示例
3
8
9
10
输出示例
-1
3 6 6 6
-1
数据规模及约定
对于第 (1) 个测试点，分值为 (30)，(1 leq b leq 32)；
对于第 (2) 个测试点，分值为 (40)，(1 leq b leq 100)；
对于第 (3) 个测试点，分值为 (30)，(1 leq t leq 200, 1 leq n leq 100000)。
题解
这题。。。woc居然是个打表。。。
当然只有 (b le 130) 的时候需要打表，剩下的可以构造。
当 (b > 130) 时，实验证明都可以用三个相同的数字举出反例，令这相同的三个数字值为 (p_1p_2)（这是两个数的乘积），那么需要让他的贪心翻车，就要想方设法让贪心搞出 (4) 位数来，即 (p_2, p_2, p_2, p_1^3) 这样的一个四位数。
那么只需要满足几个条件：
-
(p_1, p_2) 都是素数
-
(p_1p_2 < b)
-
(p_1^3 < b)
-
(p_2^2 ge b)
-
(p_1p_2 < p_1^3)
下面是打表程序：
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <algorithm>
using namespace std;
#define rep(i, s, t) for(int i = (s); i <= (t); i++)
#define dwn(i, s, t) for(int i = (s); i >= (t); i--)
int read() {
int x = 0, f = 1; char c = getchar();
while(!isdigit(c)){ if(c == '-') f = -1; c = getchar(); }
while(isdigit(c)){ x = x * 10 + c - '0'; c = getchar(); }
return x * f;
}
#define maxn 100010
#define LL long long
int prime[maxn], cp;
bool vis[maxn];
void init() {
rep(i, 2, maxn - 1) {
if(!vis[i]) prime[++cp] = i;
for(int j = 1; j <= cp && i * prime[j] < maxn; j++) {
vis[i*prime[j]] = 1;
if(i % prime[j] == 0) break;
}
}
return ;
}
int ans, now[maxn], sol[maxn];
void dfs(int up, int x, int cur) {
if(ans <= cur) return ;
if(x == 1){ ans = cur; rep(i, 1, cur) sol[i] = now[i]; return ; }
dwn(i, up, 2) if(x % i == 0) now[cur+1] = i, dfs(i, x / i, cur + 1);
return ;
}
int main() {
freopen("list.txt", "w", stdout);
init();
int reach = 140;
printf("char Ans[%d][50] = {", reach);
rep(b, 1, reach) {
int x = 1; bool ok = 0;
while(x <= 1000000) {
int tmp = ++x;
// greedy
int cnt = 0;
dwn(i, b - 1, 2) if(tmp % i == 0) tmp /= i, i++, cnt++;
// force
ans = cnt;
dfs(b - 1, x, 0);
if(ans < cnt) {
sort(sol + 1, sol + ans + 1);
printf(""%d ", ans);
rep(i, 1, ans) printf("%d%c", sol[i], i < ans ? ' ' : '"');
ok = 1;
break;
}
}
if(!ok) printf(""-1"");
printf("%s", b < reach ? ", " : "};
");
}
return 0;
}
下面是 AC 程序：
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <algorithm>
using namespace std;
#define rep(i, s, t) for(int i = (s); i <= (t); i++)
#define dwn(i, s, t) for(int i = (s); i >= (t); i--)
int read() {
int x = 0, f = 1; char c = getchar();
while(!isdigit(c)){ if(c == '-') f = -1; c = getchar(); }
while(isdigit(c)){ x = x * 10 + c - '0'; c = getchar(); }
return x * f;
}
#define maxn 100010
#define LL long long
char Ans[140][50] = {"-1", "-1", "-1", "-1", "-1", "-1", "-1", "-1", "3 6 6 6", "-1", "-1", "-1", "3 9 10 10", "3 9 10 10", "3 9 10 10", "-1", "3 12 14 14", "3 12 14 14", "3 10 15 15", "3 10 15 15", "3 10 15 15", "3 10 15 15", "3 10 15 15", "3 10 15 15", "3 10 15 15", "3 15 18 18", "3 15 18 18", "3 14 20 21", "3 15 21 21", "3 15 21 21", "3 15 21 21", "3 15 21 21", "3 15 21 21", "3 15 21 21", "3 15 21 21", "3 21 21 21", "3 21 21 21", "3 21 21 21", "3 21 21 21", "3 21 21 21", "3 21 21 21", "3 21 21 21", "3 21 21 21", "3 21 21 21", "3 21 21 21", "3 21 21 21", "3 21 21 21", "3 21 21 21", "3 21 21 21", "3 22 33 35", "3 22 33 35", "3 22 33 35", "3 22 33 35", "3 22 33 35", "3 22 33 35", "3 26 35 39", "3 26 35 39", "3 26 35 39", "3 26 35 39", "3 26 35 39", "3 26 35 39", "3 26 35 39", "3 26 35 39", "3 26 35 39", "3 26 35 39", "3 22 44 49", "3 22 44 49", "3 22 44 49", "3 22 44 49", "3 22 44 49", "3 22 44 49", "3 22 44 49", "3 22 44 49", "3 22 44 49", "3 22 44 49", "3 22 44 49", "3 22 44 49", "3 26 49 52", "3 26 49 52", "3 26 49 52", "3 26 49 52", "3 26 49 52", "3 26 49 52", "3 26 49 52", "3 26 49 52", "3 26 49 52", "3 26 49 52", "3 26 49 52", "3 26 49 52", "3 26 49 52", "3 26 49 52", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 55 55", "3 33 65 65", "3 33 65 65", "3 33 65 65", "3 33 65 65", "3 33 65 65", "3 33 65 65", "3 33 65 65", "3 33 65 65", "3 33 65 65", "3 33 65 65", "3 33 65 65", "3 33 65 65", "3 33 65 65", "3 33 65 65", "3 33 65 65", "3 33 65 65", "3 33 65 65", "3 33 65 65", "3 33 65 65"};
int prime[maxn], cp;
bool vis[maxn];
void init() {
rep(i, 2, maxn - 1) {
if(!vis[i]) prime[++cp] = i;
for(int j = 1; j <= cp && i * prime[j] < maxn; j++) {
vis[i*prime[j]] = 1;
if(i % prime[j] == 0) break;
}
}
return ;
}
int main() {
init();
int T = read();
while(T--) {
int b = read();
if(b > 140) {
bool ok = 0;
for(int i = 1; i <= cp && prime[i] <= b; i++) {
for(int j = 1; j <= cp && prime[j] <= b; j++) {
int p1 = prime[i], p2 = prime[j];
if((LL)p1 * p2 < b && (LL)p1 * p1 * p1 < b && (LL)p2 * p2 >= b && (LL)p1 * p1 > p2) {
printf("3 %lld %lld %lld
", (LL)p1 * p2, (LL)p1 * p2, (LL)p1 * p2);
ok = 1; break;
}
}
if(ok) break;
}
if(!ok) puts("-1");
}
else puts(Ans[b-1]);
}
return 0;
}