質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

90.33%

  • C

    3991questions

    C言語は、1972年にAT&Tベル研究所の、デニス・リッチーが主体となって作成したプログラミング言語です。 B言語の後継言語として開発されたことからC言語と命名。そのため、表記法などはB言語やALGOLに近いとされています。 Cの拡張版であるC++言語とともに、現在世界中でもっとも普及されているプログラミング言語です。

x=a^3+b^3を満たす異なる2つのa,bの組を探す

解決済

回答 3

投稿

  • 評価
  • クリップ 0
  • VIEW 213

vain_taka

score 1

/*
x=a^3+b^3を満たす異なる2つのa,bの組を持つxを探す
*/

#include <stdio.h>
int main(void)
{
    int x, a1, b1, y, a2, b2;

    for(x = 1; x <= 10000; x++) {    
        for(a1 = 1; a1 <= 22; a1++) {
            for(b1 = 1; b1<= 22; b1++) {
                for(y = 1;y <= 10000; y++) {
                    for(a2 = 1; a2 <= 22; a2++) {
                        for(b2 = 1; b2 <= 22; b2++) {
                            x = a1*a1*a1 + b1*b1*b1;
                            y = a2*a2*a2 + b2*b2*b2;
                            if(x == y && a1 != a2 && b1 != b2)
                            printf("%d\n",x);

                        }
                    }
                }    
            }
        }
    }

    return 0;
}


1<=x<=10000の整数xについて、x=a^3+b^3を満たす異なる整数a,bの組を2つ持つxを表示するプログラムを作ろうとしています。
例えば、1729=1^3+12^3=9^3+10^3のように異なる2つのa,bの組を持つxを探します。
じぶんの考えとしては、x=a*a*a+b*b*b, y=c*c*c+d*d*dと置いて、xとyが等しく、aとcが異なり、bとdが異なる時にそのxを表示するようにしたのですが、うまくいきません。a~dの繰り返しが22までなのは3乗して10000を超えないのが22までだからです。どのように修正すればうまくできるでしょうか。言語はCです。

  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 3

+2

元のループでは、ループ内でxyを書き換えてしまっているため、うまくループしません。というより、計算結果しか入れないのであればxyでループさせる意味が全くありません。

そして、効率のことを考えれば、

  1. まず1~22までの3乗にした値の表を作っておいて、それを2つ組み合わせる形で計算する
  2. そして、a^3 + b^3 = c^3 + d^3となる以上、a < c ≦ d < bに限定して探せばいい

というような絞り込みをやれば、探索範囲はぐっと狭くなります。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2018/05/07 14:37

    なお、試しに計算させてみたところ、条件を満たすのは1729と4104だけのようです。

    キャンセル

  • 2018/05/07 21:34

    なるほど、xとyを使う必要はなかったんですね。
    a<c<=b<dにするという発想は浮かびませんでした。非常に参考になります。

    キャンセル

+1

アルゴリズムをすべて検証したわけではありませんが、forループにxとyを使ったものがあるのに、

                            x = a1*a1*a1 + b1*b1*b1;
                            y = a2*a2*a2 + b2*b2*b2;


としてループ内の処理で値を変えてしまうと、うまくループが回らなくなるのではないですか?
xとyは計算で求める値なので、ループ変数にする必要はないと思います。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

checkベストアンサー

0

こんにちは、maisumakunさんがおっしゃるとおり

a^3 + b^3 = c^3 + d^3となる以上、a < c ≦ d < bに限定して探せばいい 

ので、for文でループする際a,b,c,dのループするスタート地点は

それぞれa = 1c = a + 1,d = cb = d + 1とすればいいでしょう。

こうすることで、必然的にa != c , b != dとなります。

ゴール地点ですが基本的に3乗の値が10000未満までという条件式を立てれればいいでしょう。
ただ、効率化を図るためにb < 10000 - (現在のaの3乗)とすれば処理速度は速くなります。dも同様です。
私は、3乗を求めるところを関数にしていますが、関数にしなくてもいいです。

#include <stdio.h>

#define MAX_NUMBER 10000

int Cubed(int);

int main(void){

    int a,b,c,d;
    int sum1,sum2;
    for(a = 1; Cubed(a) < MAX_NUMBER; a++){
       for(c = a + 1; Cubed(c) < MAX_NUMBER; c++){
           for(d = c; Cubed(d) < MAX_NUMBER - Cubed(c) ; d++){
              for(b = d + 1; Cubed(b) < MAX_NUMBER - Cubed(a); b++){
                 sum1 = Cubed(a) + Cubed(b);
                 sum2 = Cubed(c) + Cubed(d);
                 if(sum1 == sum2){
                   printf("%d ^ 3 + %d ^ 3 = %d ^ 3 + %d ^ 3 = %d\n",a,b,c,d,sum1);
                 }
              }
           }
        }
    }
    return 0;
}

int Cubed(int a){
    return a * a * a;
}


<実行結果>
1 ^ 3 + 12 ^ 3 = 9 ^ 3 + 10 ^ 3 = 1729
2 ^ 3 + 16 ^ 3 = 9 ^ 3 + 15 ^ 3 = 4104

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2018/05/07 21:38

    maisumakunさんおっしゃったことをさらに詳しく説明していただきとてもわかりやすかったです。効率化のくだりは参考になりました。

    キャンセル

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 90.33%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

同じタグがついた質問を見る

  • C

    3991questions

    C言語は、1972年にAT&Tベル研究所の、デニス・リッチーが主体となって作成したプログラミング言語です。 B言語の後継言語として開発されたことからC言語と命名。そのため、表記法などはB言語やALGOLに近いとされています。 Cの拡張版であるC++言語とともに、現在世界中でもっとも普及されているプログラミング言語です。