scanf関数の入力が繰り返される

Question

前提・実現したいこと

c言語でピタゴラス数を求めるプログラムを作っています。
入力された変数c_0に対して、
1<=a<=b<=c_0 かつ a^2+b^2=c_0^2を満たし、
その中でbが最大となるa、bを求めたいです。

発生している問題・エラーメッセージ

c_0の入力が済んでも、
printf関数による出力が行われないこと

該当のソースコード

#include <stdio.h>
#include <math.h>

int main(int argc, char *argv[]){
  int a, b, c_0;
  scanf("%d", &c_0);
  b=c_0;
  while(c_0!=sqrt(a^2+b^2) || a==0){
    b--;
    if(a==0){
      a++;
    }
  }
  printf("%d %d %d\n", a, b, c_0);
  return 0;
}

Comments

[cateye]

”a^2+b^2”は。a xor 2 + b xor 2になりますが合っていますか？→https://www.sejuku.net/blog/24512
べき乗するのはpow()関数です。→http://www.c-tipsref.com/reference/math.html
パラメータも復帰値もdoubleです。
まぁ、この場合はa*a+b*bでも良さそうですが？

[退会済みユーザー]

間違っていました。
ありがとうございます

[cateye]

あと、c_0!=sqrt()は永久ループになりかねない。
→doubleとintの比較→（暗黙の型変換で）double同士の比較に成る→正確に同じものにはならない。
→http://www.cc.kyoto-su.ac.jp/~yamada/programming/float.html
intにキャストすれば良いかも？・・・経験なし。

いっそのこと、(c_0*c_0) != (a*a+b*b)のほうが良いかも？

[dodox86]

大したことではありませんが、「関数型プログラミング」タグは違うと思います。

Answer 1

aが初期化されていなかったり、コマンドライン引数の記述の仕方をしているのにプログラムの中でその入力を受け取っている様子がないのが原因ではないでしょうか？

Answer 2

総当りになってますが・・・

#include <stdio.h>
#include <stdlib.h>

int main(int agc, char *agv[])
{
    if (agc < 2) {
        return (1);
    }
    int max_n = (int)strtol(agv[1], NULL, 10);
    /* Maximum value */
    int a_max = 0, b_max = 0, c_max = 0;
    /* from here */
    for (int c = 1; c < max_n; c++) {
        int cc = c * c;
        for (int b = 1; b < c; b++) {
            int bb = b * b;
            for (int a = 1; a < b; a++) {
                int aa = a * a;
                if (aa + bb == cc) {
//                    printf("%3d %3d %3d\n", a, b, c);
                    if (b > b_max) {
                        a_max = a;
                        b_max = b;
                        c_max = c;
                    }
                }
            }
        }
    }
    printf("a:%3d b:%3d c:%d\n", a_max, b_max, c_max);

    return (0);
}

[結果] 実行時パラメータで数値を指定。

usr ~/Project/test % time ./a.out 1000
a:140 b:975 c:985
real 0.18
user 0.17
sys 0.00
usr ~/Project/test % 

Comments

[cateye]

一万でやると３分ぐらいかかる＾＾：
usr ~/Project/test % time ./a.out 10000
a:141 b:9940 c:9941
real 171.65
user 171.56
sys 0.00