PAT 1104 天长地久 (20 分) C语言

        这个题目的关键在于对代码的简化，很容易就超时，所以一定要尽可能简化，我的简化思路放在文章最后，欢迎参考与指正。

题目：

1104 天长地久 (20 分)

“天长地久数”是指一个 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

代码： 

#include<stdio.h>
#include<math.h>
#include<stdlib.h>
typedef struct{        /*因为我没有办法使记录数据时的顺序与题目要求的输出顺序相同*/
    int sum;           /*数字A+1后的各位数字之和*/
    int number;        /*数字A*/
}numb;
int comp1(const void* a,const void* b){    /*排序函数*/
    numb* x = (numb*)a;
    numb* y = (numb*)b;
    return (x->sum != y->sum) ? x->sum - y->sum : x->number - y->number;
}
int gcd(int x,int y);    /*求取最大公因数*/
int check(int n);        /*是否为质数*/
int sum(int n);          /*求数字各位数字之和*/
int main(){
    int n , k , m, number , num_9 , cishu;    /*k是位数，m是指定和，num_9为结尾数字9的个数*/
    numb jilu[100000];                      /*记录输出数据*/
    scanf("%d",&n);
    for(int i = 1 ; i <= n; i++)
    {
        number = 0;
        printf("Case %d\n",i);        /*对case的输出*/
        scanf("%d %d", &k , &m);         /*给定k和a*/
        if(m >= k * 9)                    /*如果每一位都是9都比给定各位之和小的话，则不可能存在*/
        {
            printf("No Solution\n");
            continue;
        }
        for(num_9 = 1 ; m - num_9 * 9 + 1 > 0 ; num_9++)            /*化简的关键处*/
        {
            if( check( gcd(m , m - num_9 * 9 + 1) )==1)
            {
                cishu = 0 ;
                for(int num = pow(10,k-1) + pow(10 , num_9) - 1 ; num < pow(10 , k) ; num += pow(10 , num_9) , cishu++)
                {
                    if(cishu % 10 == 9)         //如果在第一个非9的位置（从右向左）加上9次1则整个数字会多一个九
                        continue;
                    if(sum(num) == m)     //如果符合要求
                    {
                        jilu[number].sum = sum(num+1) ;
                        jilu[number].number = num ;
                        number++ ;
                    }
                }
            }
        }
        if(number == 0)
        {
             printf("No Solution\n");
                continue;
        }
         qsort(jilu , number , sizeof(numb) , comp1);       /*按要求排序*/
         for(int x = 0;x < number ; x++)
             printf("%d %d\n",jilu[x].sum,jilu[x].number);
    }
}
int gcd(int x,int y){        /*求取最大公因数，辗转相除法*/
    int j ;
    while(y != 0)
    {
        j = x ;
        x = y ;
        y = j % x ;
    }
    return x;
}
int check(int n){
    if(n <= 2)
        return 0;
    if(n == 3)
        return 1 ;
    if(n % 6 != 1 && n % 6 != 5)
        return 0 ;
    int i;
    for(i = 6 ; i - 1 <= sqrt(n) ; i += 6)
    {
        if(n % (i + 1) == 0 || n % (i - 1) == 0)
            return 0;
    }
    return 1;
}
int sum(int n){
    int number = 0;
    do{
        number += n % 10;
        n /= 10 ;
    }while(n != 0);
    return number;
}

       化简思路： 

1.首先最后一位数字必定为9，且最后必为1或多个连续的9。

    证明：设 符合要求的数字各位数字之和为a,数字 A%10！=9，A的各位数字之和为sum，则A+1的各位数字之和为sum+1，sum和sum+1只可能有公因数1，必不成立。若a和a-9+1有质数的最大公因数，则存在最后一个数字为9的符合要求的数字，否则不存在。同理，若a和a-9*3+1有质数的最大公因数，则存在最后3个数字为9的符合要求的数字。以此类推。

2.求取最大公因数时采用辗转相除法。

        先用一个数除以另一个数,得第一个余数;再用除数除以余数,得第二个余数;又用第一个余数除以第二个余数,得第三个余数;这样逐次用除数作为被除数除以余数得到新的余数,直到余数是0为止.那么,最后一个除数就是所求的最大公约数

感谢 @herbertyellow 对代码错误的指正。