L2-018 多项式A除以B （25 分)(两种算法的比较)

这仍然是一道关于A/B的题，只不过A和B都换成了多项式。你需要计算两个多项式相除的商Q和余R，其中R的阶数必须小于B的阶数。

输入格式：

输入分两行，每行给出一个非零多项式，先给出A，再给出B。每行的格式如下：

N e[1] c[1] ... e[N] c[N]

其中N是该多项式非零项的个数，e[i]是第i个非零项的指数，c[i]是第i个非零项的系数。各项按照指数递减的顺序给出，保证所有指数是各不相同的非负整数，所有系数是非零整数，所有整数在整型范围内。

输出格式：

分两行先后输出商和余，输出格式与输入格式相同，输出的系数保留小数点后1位。同行数字间以1个空格分隔，行首尾不得有多余空格。注意：零多项式是一个特殊多项式，对应输出为0 0 0.0。但非零多项式不能输出零系数（包括舍入后为0.0）的项。在样例中，余多项式其实有常数项-1/27，但因其舍入后为0.0，故不输出。

输入样例：

4 4 1 2 -3 1 -1 0 -1
3 2 3 1 -2 0 1

 

输出样例：

3 2 0.3 1 0.2 0 -1.0
1 1 -3.1

 

 

 

 

 

 

 

 

 

 

 

#include<stdio.h>
#include<math.h>
#define N 100005
double a[N],b[N],c[N];
int Ma,Mb;
int main()
{
	int m,n,t,i,j;
	int numa=0,numc=0;
	scanf("%d",&m);
	for(i=0;i<m;i++)
	{
		scanf("%d",&t);
		scanf("%lf",&a[t]);
		if(i==0)
		{
			Ma=t;
		}
	}
	scanf("%d",&n);
	for(i=0;i<n;i++)
	{
		scanf("%d",&t);
		scanf("%lf",&b[t]);
		if(i==0)
		{
			Mb=t;
		}
	}
	
	for(i=Ma;i>=Mb;i--)   //多项式a/b的计算过程 
	{
		c[i-Mb]=a[i]/b[Mb];
		for(j=Mb;j>=0;j--)
		{
			a[i+j-Mb]-=c[i-Mb]*b[j];
		}
	}
	
	for(i=Ma;i>=0;i--)  //计算余数的有效个数 
	{
		if(fabs(a[i])>0.05)
		{
			numa++;
		}
	}
	for(i=Ma-Mb;i>=0;i--)  //计算商的有效个数 
	{
		if(fabs(c[i])>0.05)
		{
			numc++;
		}
	}
	
// 以下为输出部分 
  if(numc==0)
  {
    printf("0 0 0.0\n"); 
  }
  else
  {
  	printf("%d",numc);
  	for(i=Ma-Mb;i>=0;i--)
  	{
  		if(fabs(c[i])>0.05) printf(" %d %.1lf",i,c[i]);
  	}
  	printf("\n");
  }
	
  if(numa==0)
  {
    printf("0 0 0.0\n"); 
  }
  else
  {
  	printf("%d",numa);
  	for(i=Ma;i>=0;i--)
  	{
  		if(fabs(a[i])>0.05) printf(" %d %.1lf",i,a[i]);
  	}
  	printf("\n");
  }
}

 细心的人会发现，多项式运算部分仅几行代码便可完成，主要是数据的存储和输出。存储必须巧妙，下标得表示指数，否则计算会很困难。