POJ1060 Modular multiplication of polynomials

不想多说 全在代码里

#include <iostream>
#include <cstdio>
#include <cstring>

using namespace std;

int g[1001],f[1001],s[2001],h[1001];
int lg,ls,lh,lf;
int i,j;

void mul()  //乘
{

    for( i=0; i<lg; i++)
        for( j=0; j<lf; j++)
            s[i+j] ^= g[i]*f[j];
    ls=lg+lf-1;
    return ;

}

void mod()
{
    int t;
    int res[2001];

    while(ls>=lh)
    {
        memset(res,0,sizeof(res));
        t=ls-lh;;
        for( i=0; i<lh; i++)
        {
            if(h[i]==1)
                res[i+t]=1;
        }
        for( i=ls-1; i>=t; i--)
            s[i]^=res[i];
        for( i=ls-1; i>=0; i--)
        {
            if(s[i])
            {
                ls=i+1;
                break;
            }
            ls=0;
        }
    }
    return ;

}
int main()
{
    //freopen("in.txt","r",stdin);
    int n;
    scanf("%d",&n);
    while(n--)
    {
        memset(s, 0, sizeof(s));
		memset(f, 0, sizeof(f));
        memset(g, 0, sizeof(g));
        memset(h, 0, sizeof(h));

        scanf("%d",&lg);
        for( i=lg-1; i>=0; i--)
            scanf("%d",&g[i]);
        scanf("%d",&lf);
        for( i=lf-1; i>=0; i--)
            scanf("%d",&f[i]);
        scanf("%d",&lh);
        for( i=lh-1; i>=0; i--)
            scanf("%d",&h[i]);
        mul();
        mod();
        printf("%d",ls);
        if(ls==0)
            printf("0");
        else
            for( i=ls-1; i>=0; i--)
            {
                printf(" %d",s[i]);
            }
        printf("\n");
    }
    return 0;
}