多项式A除以B（PTA）

题目

答案

#include<iostream>
#include<math.h>
using namespace std;
const int maxn=10005;
float a[maxn];
int b1[maxn];float b2[maxn];
float c[maxn];
int main()
{
	int n1,first;
	cin>>n1;
	for(int i=0;i<n1;i++)
	{
		int tmp;
		cin>>tmp;
		cin>>a[tmp];
		if(!i) first=tmp;
	}
	int n2;
	cin>>n2;
	for(int i=0;i<n2;i++)
	cin>>b1[i]>>b2[i];
	
	for(int i=first;i>=b1[0];i--)
	{
		int tmp=i-b1[0];
		c[tmp]=a[i]/b2[0];
		for(int j=0;j<n2;j++)
		a[tmp+b1[j]]-=c[tmp]*b2[j];
	}
	int cnt1=0,cnt2=0;
	for(int i=first;i>=0;i--)
	{
		if(fabs(c[i])>=0.05) cnt1++;
		if(fabs(a[i])>=0.05) cnt2++;
	}
	
	cout<<cnt1;
	for(int i=first;i>=0;i--)
	{
		if(fabs(c[i])>=0.05) printf(" %d %.1f",i,c[i]);
	}
	if(!cnt1) cout<<" 0 0.0";
	cout<<endl;
	cout<<cnt2;
	for(int i=first;i>=0;i--)
	{
		if(fabs(a[i])>=0.05) printf(" %d %.1f",i,a[i]);
	}
	if(!cnt2) cout<<" 0 0.0";
}

参考

本文章代码基本参考了这篇文章——多项式A除以B（模拟多项式除法）

总结

关键代码就是这一段：

for(int i=first;i>=b1[0];i--)
{
	int tmp=i-b1[0];
	c[tmp]=a[i]/b2[0];
	for(int j=0;j<n2;j++)
	a[tmp+b1[j]]-=c[tmp]*b2[j];
}

int tmp=i-b1[0];是求商的对应项的指数，
c[tmp]=a[i]/b2[0];是求商的对应项的系数，
a[tmp+b1[j]]-=c[tmp]*b2[j];则是将原来的被除数进行更新——即被除数减商乘除数。