第二章 高精度 AcWing 1481. 多项式乘积

第二章 高精度 AcWing 1481. 多项式乘积

原题链接

AcWing 1481. 多项式乘积

算法标签

模拟 高精度

思路

根据多项式乘积手算过程， 进行模拟
多项式表示
a[i]表示$x^i$的系数

多项式计算
根据多项式乘积手算过程， 得知对于多项式b中每个b[i]， 需要与多项式a中每个a[j]进行计算， 计算后系数为b[i]*a[j], x的次方数为i+j

将计算结果存储至数组c中（数组c的长度最少应为 数组a最大长度 + 数组b最大长度）

代码

#pragma GCC optimize(2)
#pragma GCC optimize(3)
#include<bits/stdc++.h>
#define int long long
#define x first
#define y second
#define ump unordered_map
#define pq priority_queue
#define rep(i, a, b) for(int i=a;i<b;++i)
#define Rep(i, a, b) for(int i=a;i>=b;--i)
using namespace std;
typedef pair<int, int> PII;
const int N = 1005;
//int t, n, m, cnt, ans; 
double a[N], b[N], c[N*2];
inline int rd(){
   int s=0,w=1;
   char ch=getchar();
   while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
   while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();
   return s*w;
}
void input(double a[]){
    int k;
    scanf("%lld", &k);
    while(k--){
        int n;
        double v;
        scanf("%lld %lf", &n, &v);
        a[n]=v;
    }
}
void put(int x) {
    if(x<0) putchar('-'),x=-x;
    if(x>=10) put(x/10);
    putchar(x%10^48);
}
signed main(){
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	input(a),input(b);
	rep(i, 0, N){
	    rep(j, 0, N){
	        c[i+j]+=b[i]*a[j];
	    }
	}
	int cnt=0;
	rep(i, 0, 2*N-1){
	    if(c[i]){
	        cnt++;
	    }
	}
	printf("%lld ", cnt);
	Rep(i, 2*N-1, 0){
	    if(c[i]){
	        printf("%lld %.1lf ", i, c[i]);
	    }
	}
	return 0;
}

参考文献

AcWing 1481. 多项式乘积（PAT甲级辅导课）

原创不易
转载请标明出处
如果对你有所帮助 别忘啦点赞支持哈