7-10 多項式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
模擬就完事了
Code:
#include<bits/stdc++.h>
using namespace std;
# define ll long long
# define inf 0x3f3f3f3f
const int maxn = 1e5+100;
double a[maxn],b[maxn],c[maxn];
int main()
{
int n,m,tmp,maxa,maxb;
scanf("%d",&n);
for(int i=0; i<n; i++)
{
scanf("%d",&tmp);
scanf("%lf",&a[tmp]);
if(i==0)
maxa=tmp;
}
scanf("%d",&m);
for(int i=0; i<m; i++)
{
scanf("%d",&tmp);
scanf("%lf",&b[tmp]);
if(i==0)
maxb=tmp;
}
for(int i=maxa; i>=maxb; i--)
{
c[i-maxb]=a[i]/b[maxb];
for(int j=maxb; j>=0; j--)
{
a[i-maxb+j]-=c[i-maxb]*b[j];
}
}
//計算完成
int num1=0,num2=0;
for(int i=maxa-maxb; i>=0; i--)
{
if(fabs(c[i])<0.05)
c[i]=0;
else
num1++;
}
if(num1==0)
{
printf("0 0 0.0\n");
}
else
{
printf("%d",num1);
for(int i=maxa-maxb; i>=0; i--)
{
if(fabs(c[i])>0)
printf(" %d %.1lf",i,c[i]);
}
printf("\n");
}
for(int i=maxb-1; i>=0; i--)
{
// cout<<a[i]<<endl;
if(fabs(a[i])<0.05)
a[i]=0;
else
num2++;
}
if(num2==0)
{
printf("0 0 0.0\n");
}
else
{
printf("%d",num2);
for(int i=maxb-1; i>=0; i--)
{
if(fabs(a[i])>0)
printf(" %d %.1lf",i,a[i]);
}
printf("\n");
}
}