関数化について
何度も質問しているんですがまた行き詰ってしまいました
以前解答をいただき自分なりに進めていき以下のことができるようになりました
１、多項式の係数を入力し、多項式（Ａ）をつくる。
２、それを微分したもの（Ａ’）で割る。→Ａ÷Ａ’＝商・・・余りＢ
３、出てきた余り（Ｂ）を割る数、ひとつ前の作業で割る数だったもの（Ａ’）を割られる数にする。
４、割る数（Ｂ）の最高次の係数の２乗をしたもの(a^2)を割られる数（Ａ’）にかける。（擬除法）
５、割り算を行う。→(a^2)×Ａ’÷Ｂ
６、余りが０になるまで３～５を繰り返す。（０にならなければ終わり）
７、割り切れたときの割る数がＡとＡ’の最大公約数Ｄとなる（正確には最大公因数？）
とここまでできるようになりました。
で、次にさらにまた
８、Ａ÷Ｄ＝Ｅ
９、ＤとＥでユークリッドの互除法により最大公約数Ｇを求める
１０、Ｅ÷Ｇ＝Ｐ→ＰがＡのi次の平方因子
１１ＤをＡとしＤが無平方になるまで１～１０を繰り返す。
といったことをやらなければいけないんですが、
いい加減関数を使わないと長すぎるので同じ作業を簡略化するため関数化することにしました。
で、微分のプログラムは関数化できたんですが、ユークリッドの互除法のプログラムが難しくて関数化できません。
わかる方お願いしますm(_ _)m
以下プログラム
#include <stdio.h>
#include <math.h>
int main(void){
int m,m2,i,j,k,l,n,p,q;
int a[1000],b[1000],c[1000],d[1000],e[1000];
puts("何次の多項式ですか？");
printf("何次:"); scanf("%d",&m);
puts("多項式の係数を入力してください。");
for(i=m;i>=0;i--){
scanf("%d",&a[i]);
//e[i]=a[i];
//printf("e[i]=%d\n",e[i]);
}
derivative(&m,a,&m2,b);
printf("微分された多項式は:");
for(i=m2;i>=0;i--){
printf("%dx^%d+ ",b[i],i);
}
printf("\n\n");
q = 1;
l = 0;
while(q == 1){
printf("割られる数は：");
for(i=m;i>=0;i--){
a[i]=a[i]*(b[m2]*b[m2]);
printf("%dx^%d + ",a[i],i);
}
q = 0;
printf("\n上のは割る数の最高の次数の係数の２乗をかけたもの：b[m2]=%d\n",b[m2]);
printf("--\n");
printf("%d回目の商は\n",l+1);
for(k=m-m2;k>=0;k--){
c[k] = a[m - (m - m2 - k)] / b[m2];
printf("%dx^%d + ",c[k],k);//商の表示
j = m2;
for(i = m - (m - m2 - k);i>= m - (m - m2 - k) - m2 ;i--){
d[i]=a[i]-c[k]*b[j];
a[i] = d[i];
j = j - 1;
}
}
printf("\n");
for (k = m2;k >=0 ;k--){
a[k] = b[k];
}
printf("余り：");
for (k = m2-(m-m2);k >=0 ;k--){
b[k] = d[k];
printf("%dx^%d + ",b[k],k);
if (d[k] != 0){
q = 1;
}
}
for (k = 0;k <= 1000;k++){
d[k] = 0;
}
printf("q:%d\n",q);
p = m2 - 1;
m = m2;
m2 = p;
l = l + 1;
printf("\n");
while(b[m2] == 0){
m2 = m2 - 1;
}
if (m2 <= 0){
break;
}
}
if (q == 0){
printf("割り切れた");
}else{
printf("割り切れなかった");
}
return(0);
}
int derivative(int *m, int a[],int *m2, int b[])
{
int i;
for(i=*m;i>0;i--){
b[i-1]=i*a[i];
}
*m2=*m-1;
return(0);
}
