POJ 1060 Modular multiplication of polynomials(多项式的加减乘除，除法转化成减法来求)

题意：给出f(x),g(x),h(x)的 （最高次幂+1）的值，以及它们的各项系数，求f(x)*g(x)/h(x)的余数。

　　这里多项式的系数只有1或0，因为题目要求：这里多项式的加减法是将系数相加/减后再模2，这样其实也就可以用异或运算来代替加减法。

思路：看代码吧，水题一个，主要在于把除法转化成减法，一次一次减就行。

#include <iostream>
#include <stdio.h>
#include <string.h>

using namespace std;
const int maxn=2010;
//f,g,h存储的是多项式的系数,sum存储的是f*g的系数以及最后余数的系数
int f[maxn],g[maxn],h[maxn],sum[maxn];
int lf,lg,lh,ls;//分别为f,g,h,sum的最高次幂

//比较sum表示的多项式与h表示的多项式的大小
int compare() {
    if(ls<lh)
        return -1;
    if(ls>lh)
        return 1;
    for(int i=ls-1; i>=0; i--) {
        if(sum[i]==h[i])
            continue;
        if(sum[i]>h[i])
            return 1;
        if(sum[i]<h[i])
            return -1;
    }
    return 0;
}
int main() {
    int t,d;
    scanf("%d",&t);
    while(t--) {
        memset(h,0,sizeof(h));
        memset(sum,0,sizeof(sum));
        //将f多项式的信息存入f数组
        scanf("%d",&d);
        lf=d-1;
        for(int j=lf; j>=0; j--) {
            scanf("%d",&f[j]);
        }
        //将g多项式的信息存入g数组
        scanf("%d",&d);
        lg=d-1;
        for(int j=lg; j>=0; j--) {
            scanf("%d",&g[j]);
        }
        //将h多项式的信息存入h数组
        scanf("%d",&d);
        lh=d-1;
        for(int j=lh; j>=0; j--) {
            scanf("%d",&h[j]);
        }
        //计算f*g的多项式
        ls=lf+lg;
        for(int i=lf; i>=0; i--) {
            for(int j=lg; j>=0; j--) {
                sum[i+j]=sum[i+j]^(f[i]&g[j]);
            }
        }
        /*
          关键是怎么求余数，这里是先判断sum多项式是否大于h多项式，
          若大于，则sum减一次h，减去后的信息存入sum中。
          再继续判断,直到sum小于h，则此时的sum为余数。
          总之，就是把除法改成较简单的减法操作。
        */
        while(compare()>=0) {
            d=ls-lh;
            for(int i=ls; i-d>=0; i--) {
                sum[i]=sum[i]^h[i-d];
            }
            while(ls && !sum[ls])
                ls--;
            /*
            原先一直WA的代码,在每次更新sum的最高次幂ls时出了错误。
            int mark=0;
            for(int i=ls; i-d>=0; i--) {
                sum[i]=sum[i]^h[i-d];

                下面错误错在这里只判断了i>=d的情况，
                有可能当i>=d的时候sum[i]=0,这样求出mark(也就是结果的最高次幂)为0；
                但是可能有i<d的时候，有sum[i]=1,那么mark就不为0。

                if(sum[i] && !mark) {
                    mark=i;
                }

            }
            ls=mark;
            */

        }

        printf("%d",ls+1);
        for(int i=ls; i>=0; i--) {
            printf(" %d",sum[i]);
        }
        printf("
");
    }
    return 0;
}