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;
}

猜你喜欢

转载自www.cnblogs.com/IamIron-Man/p/11930165.html