“天长地久数”是指一个 K 位正整数 A,其满足条件为:A 的各位数字之和为 m,A+1 的各位数字之和为 n,且 m 与 n 的最大公约数是一个大于 2 的素数。本题就请你找出这些天长地久数。
输入格式:
输入在第一行给出正整数 N(≤5),随后 N 行,每行给出一对 K(3<K<10)和 m(1<m<90),其含义如题面所述。
输出格式:
对每一对输入的 K 和 m,首先在一行中输出 Case X,其中 X 是输出的编号(从 1 开始);然后一行输出对应的 n 和 A,数字间以空格分隔。如果解不唯一,则每组解占一行,按 n 的递增序输出;若仍不唯一,则按 A 的递增序输出。若解不存在,则在一行中输出 No Solution。
输入样例:
2
6 45
7 80
输出样例:
Case 1
10 189999
10 279999
10 369999
10 459999
10 549999
10 639999
10 729999
10 819999
10 909999
Case 2
No Solution
代码长度限制16 KB
时间限制3000 ms
内存限制64 MB
3、4测试点超时:
这道题别的没什么好说的,就有一个坑点在于它给的整数长度最大能有K=9,如果从10的8次方到10的9次方都遍历,那就一定会超时
这是我的超时代码:
#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<math.h>
typedef struct
{
int geshu;
int shuzi;
}STU;
int gcd(int a, int b)
{
return b == 0 ? a : gcd(b, a % b);
}
int isprime(int i)
{
int j;
if (i == 1 || i <= 0)
return 0;
else if (i == 2)
return 1;
for (j = 2; j <= sqrt(i); j++)
{
if (i % j == 0)
return 0;
}
return 1;
}
int cmp(const void* a, const void* b)
{
STU* a1 = (STU*)a; STU* b1 = (STU*)b;
if (a1->geshu != b1->geshu)
return a1->geshu - b1->geshu;
else
return a1->shuzi - b1->shuzi;
}
int main()
{
int n;
scanf("%d", &n);
int i, j, k, m, ch, ch1;
for (i = 0; i < n; i++)
{
int q = 0;
STU stu[10000];
scanf("%d %d", &k, &m);
for (j = (int)pow(10, k - 1); j < (int)pow(10, k); j++)//最无脑的遍历方式
{
ch = j;
int sum = 0, sum1 = 0;
int t;
while (ch != 0)
{
sum = sum + ch % 10;
ch = ch / 10;
}
if (sum == m)
{
ch1 = j + 1;
while (ch1 != 0)
{
sum1 = sum1 + ch1 % 10;
ch1 = ch1 / 10;
}
t = gcd(m, sum1);
if (isprime(t) == 1 && t > 2)
{
stu[q++].geshu = sum1;
stu[q - 1].shuzi = j;
}
}
}
if (q == 0)
{
printf("Case %d\n", i + 1);
printf("No Solution\n");
}
else
{
printf("Case %d\n", i + 1);
qsort(stu, q, sizeof(STU), cmp);
int r;
for (r = 0; r < q; r++)
{
printf("%d %d\n", stu[r].geshu, stu[r].shuzi);
}
}
}
return 0;
}
这里我们只需要把这一句
for (j = (int)pow(10, k - 1); j < (int)pow(10, k); j++)
修改成这样:
for (j = (int)99 + pow(10, k - 1); j < (int)pow(10, k); j += 100)
正确的程序代码:
#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<math.h>
typedef struct
{
int geshu;
int shuzi;
}STU;
int gcd(int a, int b)
{
return b == 0 ? a : gcd(b, a % b);
}
int isprime(int i)
{
int j;
if (i == 1 || i <= 0)
return 0;
else if (i == 2)
return 1;
for (j = 2; j <= sqrt(i); j++)
{
if (i % j == 0)
return 0;
}
return 1;
}
int cmp(const void* a, const void* b)
{
STU* a1 = (STU*)a; STU* b1 = (STU*)b;
if (a1->geshu != b1->geshu)
return a1->geshu - b1->geshu;
else
return a1->shuzi - b1->shuzi;
}
int main()
{
int n;
scanf("%d", &n);
int i, j, k, m, ch, ch1;
for (i = 0; i < n; i++)
{
int q = 0;
STU stu[10000];
scanf("%d %d", &k, &m);
for (j = (int)99 + pow(10, k - 1); j < (int)pow(10, k); j += 100)
{
ch = j;
int sum = 0, sum1 = 0;
int t;
while (ch != 0)
{
sum = sum + ch % 10;
ch = ch / 10;
}
if (sum == m)
{
ch1 = j + 1;
while (ch1 != 0)
{
sum1 = sum1 + ch1 % 10;
ch1 = ch1 / 10;
}
t = gcd(m, sum1);
if (isprime(t) == 1 && t > 2)
{
stu[q++].geshu = sum1;
stu[q - 1].shuzi = j;
}
}
}
if (q == 0)
{
printf("Case %d\n", i + 1);
printf("No Solution\n");
}
else
{
printf("Case %d\n", i + 1);
qsort(stu, q, sizeof(STU), cmp);
int r;
for (r = 0; r < q; r++)
{
printf("%d %d\n", stu[r].geshu, stu[r].shuzi);
}
}
}
return 0;
}