C言語で友愛数を求めたいです。

Question

前提・実現したいこと

C言語を用いて友愛数を求めるプログラムを書きたいです。
コンパイルは通りましたが結果が表示されません。

以下、問題です。
友愛数を求めるプログラムを作成し、１００００以下の範囲で友愛数を求めよ。
友愛数：異なる2つの自然数の組で、自分自身を除いた約数の和が、互いに他方と等しくなるような数のこと。

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

[xxxxxxx@xxxxxxxx kadai1]$ ./a.out
[xxxxxxx@xxxxxxxx kadai1]$

該当のソースコード

　#include<stdio.h>

//ぷろとたいぷ宣言
int yuuaisuu(int);

int main(){
  int i,j;
  int n = 0;         //sumの配列番号である変数
  int number;        //他方の数の一時退避領域
  int sum[10000]={}; //nの約数である変数
  int yuuaikentei;   //sum[n]から自分自身を引いた数を入れる変数

  for(i = 1 ;i < 10000; i++){ 
   
    for(j = 1;j < i; j++){
      if( (i % j) == 0){
	sum[n] += j;
      }
    }

    /* 友愛数は約数から自分自身をひいたものが他方と等しくなる数なので
       相互に計算する（もう一つはif条件内で計算している）           */
    yuuaikentei = sum[n]-n;
    number = yuuaisuu(yuuaikentei);
   
    if(yuuaikentei == number){
      printf("%d , %d\n",yuuaikentei,number);
    }
    /* ここで配列の番号を変化させる */
    n++;
  }

  return 0;
}


int yuuaisuu(int a){
  int sum=0;
  int i,j;
  /* 約数を求める */
  for( j = 1; j < a ; j++){
    if( (a % j) == 0){
      sum += j;
   }
  }
  /* 約数から自分自身を引く */
  i = sum-j;
  return i;
}

試したこと

1000など小さい数字でも試してみましたが、同様に表示されませんでした。

Answer 1

アルゴリズムも全然違いますし、作り方が悪いです。

まず自分自身を除いた約数の和を求める関数 int sumOfFactors(int) を作ってください。
そして 2 から 20 程度までの数を引数にその関数を呼び出して出力を確認し、正しい出力が得られることを確認してください。
ちなみにこのような結果になります。
1 1 3 1 6 1 7 4 8 1 16 1 10 9 15 1 21 1 22
これ以外の結果が出たら間違いなので作り直してください。

次に 2 から 10000 までの数を引数にこの関数を呼び出し、得られた結果を引数に再度呼び出して得られた結果と元の引数が同じものを抽出してください。

それが成功した後は、重複して現れた同じ結果を省くよう工夫してください。
結果は次のようになります。

220 284 1184 1210 2620 2924 5020 5564 6232 6368

Answer 2

github　上で次のプログラムをみつけました。

• Nまでの友愛数の組を表示するプログラム

https://gist.github.com/kaiware/7652747

これを少し変更してみました。
num.c

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

#define N (100000)

int div_sum(int x) {
  int m = sqrt(x);
  int sum = 0;
  for (int i = 2; i <= m; i++) {
    if (x % i == 0) {
      sum += (i + x / i);
    }
  }

  if (m * m == x) {
    sum -= m;
  }
  return sum + 1;
}

int main() {
  for (int a = 2; a < N; a++) {
    int sum_a = div_sum(a);
    if (sum_a > a) {
      if (a == div_sum(sum_a)) {
        printf("(%d, %d)\n", a, sum_a);
      }
    }
  }
  return 0;
}

実行例

Comments

[langage]

解決できました。ご回答ありがとうございました

Answer 3

コンパイルが通ったからと言って期待する挙動はしません(書いたとおり動くので結果が出なかったりしたらプログラムが間違ってます)

#include<stdio.h>
#include <string.h>

#define N (10000)

int main(void)
{
	int sum[N + 1];
	int i = 0, j = 0;
	
	memset(sum, 0, sizeof(sum));
	
	for(i = 1; i <= N; i++) {
		for(j = 1; j <= i / 2; j++) {
			if((i % j) == 0) {
				sum[i] += j;
			}
		}
	}
	
	for(i = 1; i <= N; i++) {
		if(i != sum[i]) {
			if(i == sum[sum[i]]) {
				printf("%dと%d\n", i, sum[i]);
			}
		}
	}
	
    return 0;
}

あえて落ちるよう＆ダブってるのも出してますので後は頑張ってください。