XTU 1352 Fraction

题目描述 埃及人只能理解分子为1的分数，所以他们把分子不为1的分数理解成若干个分母不相同的分数的累加和。 比如23=12+16。 如果存在多个解，我们希望按分母从小到大的排列，然后取字典序最小的一个。 比如35=12+110=13+14+160。 我们按分母从小到大排列，前者为<2,10>,后者为<3,4,60>,显然前者的字典序更小。 现在给你一个分数，请求出符合要求的埃及人能理解的分数。 输入 第一行是一个整数T(1≤T≤1000)，表示样例的个数， 每个样例是两个整数a,b(1≤a<b≤50)，表示分数ab。 输出 每行输出一个样例的结果，为分母从小到大排列，每个分母之间用一个空格隔开。 样例输入 3 1 2 2 3 3 7 样例输出 2 2 6 3 11 231

 

初始代码：

#include<stdio.h>
int main()
{
    int T;
    scanf("%d", &T);
    while (T--) {
        int a, b;
        scanf("%d %d", &a, &b);
        int cnt = 2;
        int fz = a;
        int fm = b;
        while (fz != 1) {
            int ret = 0;
            if (fz * cnt - fm >= 0) {
                fz = fz * cnt - fm;
                fm *= cnt;
                ret = 1;
            }
            if (ret)printf("%d ", cnt);
            cnt++;
        }
        printf("%d\n", fm);
    }
    return 0;
}

 遇到困难：

1、超时，算法的优化

2、类型的选择

3、整数的干扰（可能会让除数为零）

修改后代码：

#include<stdio.h>
long long S(long long a, long long b);
int main()
{
    int T;
    scanf("%d", &T);
    while (T--) {
        long long a, b, t;
        scanf("%lld %lld", &a, &b);
        t = S(a, b);
        a /= t;
        b /= t;
        long long m = b / a;
        int ret = 1;
        double n = (double)b / a;
        long long fm ;
        if (n == m)fm = m;
        else fm = m + 1;
        while (a != 1) {
            printf("%lld ", fm);
            a = a * fm - b;
            b *= fm;
            t = S(a, b);
            a /= t;
            b /= t;
            if (a != 0)m = b / a;
            else {
                ret = 0;
                break;
            }
            n = (double)b / a;
            if (n == m)fm = m;
            else fm = m + 1;
        }
        if (ret)printf("%lld\n", fm);
        else printf("\n");
    }
    return 0;
}
long long S(long long a, long long b)
{
    long long t;
    while (b != 0) {
        t = a % b;
        a = b;
        b = t;
    }
    return a;
}