PTA-多项式A除以B

多项式A除以B

这仍然是一道关于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

知道我看到这个有多高兴吗？！我调试了一晚上呀555，这是纯手打亲自想的一个方法，总结一下，首先，整型是integer型，32767，然后测试的时候关联文件都加上吧，看题给的格式不要错了，然后输入的时候用scanf是快，但是你得想好数据是什么格式就%什么，虽然输的是整数，但是给double输，就得先lf，不然后面一堆奇怪错误我也不知道为啥，然后四舍五入这么写没问题：

    	ans_xishu[i] = (double) ((int)(ans_xishu[i]*10 + (ans_xishu[i]<0?-0.5:0.5)))/10;

然后还有一堆粗心，比如都算完了然后输出没算的那个数，排序的时候排指数而不是系数，等等，设计测试样例！细心！我也不知道这道题花这么多时间值不值得，可能应该在难题上练代码调试，这种题看看正解？明天看看大佬们都怎么写的

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<cmath>
#include<algorithm>
using namespace std;
struct dxs{
	int zhishu;
	double xishu;
}c[33000];

bool cmp(dxs d1,dxs d2)
{
	return d1.zhishu>d2.zhishu;
}	

int main()
{
//	freopen("test.in","r",stdin);
//	freopen("test.out","w",stdout);
	int a_zhishu[33000],b_zhishu[33000],temp_zhishu[33000],shuchu;
	double temp_xishu[33000],b_xishu[33000];
	double ans_xishu[33000],a_xishu[33000];
	int ans_zhishu[33000];
	bool temp_biaoji[33000],pd,vo;
	int n1,n2,i,num=0,start,tot=0;
	cin>>n1;
	for (i=1;i<=n1;i++)
	  scanf("%d%lf",&a_zhishu[i],&a_xishu[i]);
	cin>>n2;
	for (i=1;i<=n2;i++)
	  scanf("%d%lf",&b_zhishu[i],&b_xishu[i]);
	while (n1!=0 and a_zhishu[1]>=b_zhishu[1])
	{
		start=1;
		num++;
		ans_zhishu[num]=a_zhishu[1]-b_zhishu[1];
		ans_xishu[num]=a_xishu[1]/b_xishu[1];
		for (i=1;i<=n2;i++)
		{
			temp_zhishu[i]=ans_zhishu[num]+b_zhishu[i];
			temp_xishu[i]=ans_xishu[num]*b_xishu[i];
			temp_biaoji[i]=0;
		}
		for (i=1;i<=n1;i++)
		{
		  pd=false;
		  for (int j=start;j<=n2 and temp_zhishu[j]>=a_zhishu[i];j++)
		  	if (a_zhishu[i]==temp_zhishu[j])
		  	{
		  		pd=true;
		  		if (a_xishu[i]!=temp_xishu[j])
		  		{
		  			c[++tot].zhishu=a_zhishu[i];
		  			c[tot].xishu=a_xishu[i]-temp_xishu[j];
		  	    }
				temp_biaoji[j]=1;
		  		start++;
		  		break;
		  	}
		  if (pd==false)
		  {
		  		c[++tot].zhishu=a_zhishu[i];
		  		c[tot].xishu=a_xishu[i];
		  }
	    }
		for (i=1;i<=n2;i++)
		  if (temp_biaoji[i]==0)
		  {
		  	c[++tot].zhishu=temp_zhishu[i];
		  	c[tot].xishu=(-1)*temp_xishu[i];
		  }
		sort(c+1,c+tot+1,cmp);
		for (int i=1;i<=tot;i++)
		{
			a_xishu[i]=c[i].xishu;
			a_zhishu[i]=c[i].zhishu;
		}
		n1=tot;
		tot=0;
	}
	for(int i=1;i<=num;i++)
    	ans_xishu[i] = (double) ((int)(ans_xishu[i]*10 + (ans_xishu[i]<0?-0.5:0.5)))/10;
    for(int i=1;i<=n1;i++)
    	a_xishu[i] = (double)((int)(a_xishu[i]*10 + (a_xishu[i]<0?-0.5:0.5)))/10;
	shuchu=num;
	for (i=1;i<=num;i++)
	  if (ans_xishu[i]==0)
	    shuchu--;
	if (shuchu==0)
	  cout<<"0 0 0.0"<<endl;
	else
	{
		vo=true;
	  printf("%d ",shuchu);
	  for (i=1;i<=num;i++)
	    if (ans_xishu[i]!=0)
	      if (vo==true)
	      {
	      	printf("%d %.1lf",ans_zhishu[i],ans_xishu[i]);
			vo=false;
		  }
		  else
	  	  printf(" %d %.1lf",ans_zhishu[i],ans_xishu[i]);
	  printf("\n");
    }
    shuchu=n1;
    for (i=1;i<=n1;i++)
      if (a_xishu[i]==0)
        shuchu--;
	if (shuchu==0)
	  cout<<"0 0 0.0";
	else
	{
	  printf("%d ",shuchu);
	  vo=true;
	  for (i=1;i<=n1;i++)
	    if (a_xishu[i]!=0)
	      if (vo==true)
	      {
	        printf("%d %.1lf",a_zhishu[i],a_xishu[i]);
	        vo=false;
	      }
	      else
	  	    printf(" %d %.1lf",a_zhishu[i],a_xishu[i]);
    }
	return 0;
}

今天看了大佬的题解，是一个浙大的学姐好像是（www.liuchuo.net），果然厉害，看她博客感触也颇深，大学过得真的充实、佩服！不仅算法、比赛贼强，一堆证书，还自己写了几本书，不仅会敲代码还学习处事等等都非常好，拿博客记录了大学生活，超级强。她就一个核心过程就完成了，下标是指数，内容是系数，系数为零也填充，然后A[i]直接减当前算出来的那一项商，a[i]-=a[t1]/b[t2]（系数）*b[i-(t1-t2)]（t1-t2是当前最高次幂指数，就能保证一一对应了），然后从高到低循环减（for），一直到t1<t2为止（while），捧上她的代码：（原文地址：https://www.liuchuo.net/archives/5585）

#include <cstdio>
#include <cmath>
using namespace std;
int nonNegativeNum(double c[], int start) {
    int cnt = 0;
    for (int i = start; i >= 0; i--)
        if (abs(c[i]) + 0.05 >= 0.1) cnt++;
    return cnt;
}
void printPoly(double c[], int start) {
    printf("%d", nonNegativeNum(c, start));
    if (nonNegativeNum(c, start) == 0) printf(" 0 0.0");
    for (int i = start; i >= 0; i--)
        if (abs(c[i]) + 0.05 >= 0.1)
            printf(" %d %.1f", i, c[i]);
}
double c1[3000], c2[3000], c3[3000];
int main() {
    int m = 0, n = 0, t = 0, max1 = -1, max2= -1;
    scanf("%d", &m);
    for (int i = 0; i < m; i++) {
        scanf("%d", &t);
        max1 = max1 > t ? max1 : t;
        scanf("%lf", &c1[t]);
    }
    scanf("%d", &n);
    for (int i = 0; i < n; i++) {
        scanf("%d", &t);
        max2 = max2 > t ? max2 : t;
        scanf("%lf", &c2[t]);
    }
    int t1 = max1, t2 = max2;
    while (t1 >= t2) {
        double c = c1[t1] / c2[t2];
        c3[t1 - t2] = c;
        for (int i = t1, j = t2; j >= 0; j--, i--) c1[i] -= c2[j] * c;
        while (abs(c1[t1]) < 0.000001) t1--;
    }
    printPoly(c3, max1 - max2);
    printf("\n");
    printPoly(c1, t1);
    return 0;
}