多倍長の乗算について。

Question

こんにちは。現在、多倍長演算として、乗算(掛け算)を作っているのですが、うまくいかないです。プログラムは以下にありますが、a*b=cにしたいです。
構造体の定義はこちらです。
struct NUMBER
{
int n[KETA];/KETA個の桁/
int sign;/符号、正なら1、負なら-1/
};
clearByZero(a)は、配列の中身を全て0にします。
mulBy10(b,c)は、bを10倍したものをcに代入します。
add(d,e,f)は、dとeを足したものをfに代入します。
これらは正しく作動するものとします。(というかしています。)

---

int multiple(struct NUMBER *a,struct NUMBER *b,struct NUMBER *c)
{
struct NUMBER d;
int asign=a->sign;
int bsign=b->sign;
int i,j;
int h=0;
int x;

clearByZero(c);

a->sign=1;
b->sign=1;

for(i=0;i<KETA-1;i++)
{
for(j=0;j<KETA;j++)
{
d.n[j]=((a->n[j]*b->n[i])+h)%10;
h=((a->n[j]*b->n[i])+h)/10;
}
for(x=i;x>0;x--)
{
mulBy10(&d,&d);
}
add(c,&d,c);
}

if(asign==bsign)
{
c->sign=1;
}
else
{
c->sign=-1;
}
return (0);
}

実行結果としては、どうしようと0になります。
正負の扱いに関しては上手くいっています。

添削の程よろしくお願いいたします。

Comments

[退会済みユーザー]

他のユーザから「意図的に内容が抹消された質問」という指摘を受けました
解決後に編集機能を用いて質問内容を改変し関係のない内容にしたり、内容を削除する行為は禁止しています。
投稿していただいた質問は、後に他の誰かが困ったときに助けになる情報資産になると考えるからです。
「質問を編集する」ボタンから編集を行い、他のユーザにも質問内容が見えるように修正してください。

[shanxia]

他のユーザから「意図的に内容が抹消された質問」という指摘を受けました
解決後に編集機能を用いて質問内容を改変し関係のない内容にしたり、内容を削除する行為は禁止しています。
投稿していただいた質問は、後に他の誰かが困ったときに助けになる情報資産になると考えるからです。
「質問を編集する」ボタンから編集を行い、他のユーザにも質問内容が見えるように修正してください。

[coco_bauer]

プログラムの最後から２行目が、return (0); となっていますから、返り値は当然０になるでしょう。引数として渡した c の値も、実行後に常に０になるのですか？

Answer 1

解凍してくださった方々ありがとうございました。

Answer 2

多倍長演算のルーチンを「使うこと」が目的なら、boostライブラリに多倍長演算のライブラリがあります。
boost多倍長整数
boost多倍長浮動小数点数

ルーチンを「自作すること」が目的ならすいません。

Comments

[Gogetsu_azusa]

解答ありがとうございます。
そうですね、自作したいです。

Answer 3

こんにちは。

KETA=3、a=12, b=2でやってみましたら、c=24と適切な結果を得ることができました。

これらは正しく作動するものとします。(というかしています。)

本当に正しく動作していますでしょうか？
下記は入力と出力に同じ領域を指定していますが、そのような使い方でも正しく動作することを確認されてますか？

mulBy10(&d,&d);
add(c,&d,c);

しかし、インプレース処理を考慮してなくても、結果が0になることはあまりないです。
signについては確認されているようなので、coco_bauerさんの指摘が当たっていると言うことはないとは思うのですが．．．

---

【追記】
動作確認したソースを上げておきます。
NUMBERとmultipleの内容は修正していません。

#include <iostream>

#define KETA 3

struct NUMBER 
{ 
  int n[KETA];/*KETA個の桁*/ 
  int sign;/*符号、正なら1、負なら-1*/ 
}; 

void clearByZero(struct NUMBER *c)
{
    for (int i=0; i < KETA; ++i)
        c->n[i]=0;
    c->sign=0;
}
void mulBy10(struct NUMBER *a, struct NUMBER *b)
{
    for (int i=KETA-1; 0 < i; --i)
        b->n[i]=a->n[i-1];
}
void add(struct NUMBER *a, struct NUMBER *b, struct NUMBER *c)
{
    int h=0;
    for (int i=0; i < KETA; ++i)
    {
        int x=a->n[i]+b->n[i]+h;
        c->n[i]=x%10;
        h=x/10;
    }
}

void print(char const* t, struct NUMBER *a)
{
	std::cout << t;
    if (a->sign > 0) {
        std::cout << "+ ";
    } else if (a->sign < 0) {
        std::cout << "- ";
    } else {
        std::cout << "  ";
	}
    for (int i=KETA-1; 0 <= i; --i)
        std::cout << a->n[i] << " ";
    std::cout << "\n";
}

int multiple(struct NUMBER *a,struct NUMBER *b,struct NUMBER *c) 
{ 
  struct NUMBER d; 
  int asign=a->sign; 
  int bsign=b->sign; 
  int i,j; 
  int h=0; 
  int x;

  clearByZero(c);

  a->sign=1; 
  b->sign=1;

  for(i=0;i<KETA-1;i++) 
  { 
    for(j=0;j<KETA;j++) 
    { 
      d.n[j]=((a->n[j]*b->n[i])+h)%10; 
      h=((a->n[j]*b->n[i])+h)/10; 
    } 
    for(x=i;x>0;x--) 
    { 
      mulBy10(&d,&d); 
    } 
    add(c,&d,c); 
  }

  if(asign==bsign) 
  { 
    c->sign=1; 
  } 
  else 
  { 
    c->sign=-1; 
  } 
  return (0); 
}

int main()
{
    NUMBER a={{2,1,0},1};
    NUMBER b={{2,0,0},1};
    NUMBER c;
    clearByZero(&c);

    print("a=", &a);
    print("b=", &b);
    print("c=", &c);
    multiple(&a, &b, &c);
	std::cout << "\n";
    print("c=", &c);

    return 0;
}

MinGW 5.2.0とmsvc 2015で下記のように出力されます。

a=+ 0 1 2
b=+ 0 0 2
c=  0 0 0

c=+ 0 2 4

コピペ＋やっつけなので、変数名等は適当です。

Comments

[Gogetsu_azusa]

解答ありがとうございました！
おっしゃる通りで、add(足し算)に問題がありました。
a=a+bにできませんでした。
直したら上手くいきました！