Morris' Blog 我在摸索我存在的意義、生命的意義、內涵及價值。 沒有真正神手般的厲害，套用模式是大部分的情況。 沒那麼厲害的我存在的意義為何，襯托出神的存在？ 有時候不能怪別人不給我機會，我已經開始省思， 為什麼我沒有給別人機會的原因，而是全要靠別人。 有時候我認為我與眾不同，那是因為我做不到很多人能辦到的事情。 我聽不到、玩不起來、說得不夠明瞭、理解得不夠多， 我缺少的語文能力，就是一切一切無能的判定。 為此，我不曉得能做些對別人有幫助的事情。 我知道我自己是個不服輸的性格，但現在我不得不認輸。 我沒有能力，因此我必須認輸。 我到底來這世界作什麼？處處充滿不確定性， 而我總是無法出類拔萃，想當個一般人似乎也回不了頭了。 越強大的挑戰需要越強力的支持，而我的支持是什麼？ 最近嘗試與人談話，但都是失敗的結局，沒人能接受， 不需要言語，不需要溝通，一切所想都表達不出來。

49愛的鼓勵 6訂閱站台

a247. Fermat vs. Pythagoras

內容 ：

輸入說明 ：

輸出說明 ：

範例輸入 ：

範例輸出 ：

提示 ：

出處 ：

---

作法 : math + 建表
/**********************************************************************************/
/*  Problem: a247 "Fermat vs. Pythagoras" from ACM 106                            */
/*  Language: C (1313 Bytes)                                                      */
/*  Result: AC(101ms, 12MB) judge by this@ZeroJudge                               */
/*  Author: morris1028 at 2011-10-16 08:53:28                                     */
/**********************************************************************************/


#include<stdio.h>      
#include<stdlib.h>
#include<string.h>
#define MaxN 1000000
#define oo 2147483647
int gcd(int x, int y) {
    int tmp;
    while(x%y) {
        tmp = x, x = y, y = tmp%y;
    }
    return y;
}
int mask1[MaxN+1], mask2[MaxN+1], Ans[MaxN+1];
int min(int x, int y) {
    return x < y ? x : y;
}
int Build() {

    int n, m, i, j, k, a, b, c;    
    for(i = 0; i <= MaxN; i++)
        mask1[i] = oo;
    memset(mask2, 0, sizeof(mask2));
    memset(Ans, 0, sizeof(Ans));
    for(i = 1; i < 1000; i++)
        for(j = i+1; j < 1000; j += 2) {
            if(gcd(i, j) > 1)    continue;
            a = j*j - i*i, b = 2*i*j, c = j*j + i*i;
            if(c > MaxN)    break;
            mask2[c]++;
            for(k = 1; k*c <= MaxN; k++) {
                mask1[k*a] = min(k*c, mask1[k*a]);
                mask1[k*b] = min(k*c, mask1[k*b]);
                mask1[k*c] = min(k*c, mask1[k*c]);
            }
        }
    for(i = 1; i <= MaxN; i++)
        if(mask1[i] <= MaxN)
            Ans[mask1[i]]++;
    for(i = 1; i <= MaxN; i++) {
        mask2[i] += mask2[i-1];
        Ans[i] += Ans[i-1];
    }
}
int main() {   
    int N;
    Build();
    while(scanf("%d", &N) == 1&& N)    {
        printf("%d %d\n", mask2[N], N-Ans[N]);
    }   
   return 0;   
}   
/*利用此一性質，計算所有斜邊小於10000000的素勾股數，*/  
/* abc為三邊  m > n 、 m 和 n 均是正整*/  
/*a = m^2 -n^2, */  
/*b = 2mn, */  
/*c = m^2+n^2*/