学步园

UVA_138

    这个题目让我第一次过了一把打表的瘾，自从悠然告诉我真正的打表法该怎么用了以后，这是第一次实践。

    当然后来又学了一种佩尔方程解法，根据已知我们可以化简得到一个二元一次方程(2n+1)^2-8*m^2=1，然后找到基本解2n+1=3和m=1，依据递推公式找到后面的10对m、n打印出来即可。

//打表程序
#include<stdio.h>
#include<string.h>
int main()
{
    freopen("UVA_138(1).cpp", "w", stdout);
long long int i, j, min, mid, max;
    printf("#include<stdio.h>\n");
    printf("int main()\n");
    printf("{\n");
for(i = 6; i <= 65918161; i ++)
    {
        min = 1;
        max = i + 1;
for(;;)
        {
            mid = (min + max) / 2;
if(mid == min)
break;
if(mid * mid * 2 > i * (i + 1))
                max = mid;
else
                min = mid;
        }
if(mid * mid * 2 == i * (i + 1))
            printf("\tprintf(\"%10lld%10lld\\n\");\n", mid, i);
    }
    printf("\treturn 0;\n");
    printf("}\n");
return 0;
}

 

#include<stdio.h>
int main()
{
    printf("         6         8\n");
    printf("        35        49\n");
    printf("       204       288\n");
    printf("      1189      1681\n");
    printf("      6930      9800\n");
    printf("     40391     57121\n");
    printf("    235416    332928\n");
    printf("   1372105   1940449\n");
    printf("   7997214  11309768\n");
    printf("  46611179  65918161\n");
return 0;
}

#include<stdio.h>
#include<string.h>
int main()
{
int i = 0, j, x1, y1, x, y, newx, newy;
    x = x1 = 3, y = y1 = 1;
for(i = 0 ;i < 10; i ++)
    {
        newx = x * x1 + 8 * y * y1;
        newy = y * x1 + x * y1;
        x = newx, y = newy;
        printf("%10d%10d\n", y, (x - 1) / 2);
    }
return 0;
}