多项式除法及求模的计算复杂度

这里写图片描述

这里写图片描述

void polynomial_division(int n, int m, long long *A, long long *B, long long *D, long long *R)
{
    static long long A0[MaxN], B0[MaxN];

    int p = 1, t = n - m + 1;
    while(p < t << 1) p <<= 1;

    fill(A0, A0 + p, 0);
    reverse_copy(B, B + m, A0);
    polynomial_inverse(t, A0, B0);
    fill(B0 + t, B0 + p, 0);
    transform(p, B0);

    reverse_copy(A, A + n, A0);
    fill(A0 + t, A0 + p, 0);
    transform(p, A0);

    for(int i = 0; i != p; ++i)
        A0[i] = A0[i] * B0[i] % mod_v;
    inverse_transform(p, A0);
    reverse(A0, A0 + t);
    copy(A0, A0 + t, D);

    for(p = 1; p < n; p <<= 1);
    fill(A0 + t, A0 + p, 0);
    transform(p, A0);
    copy(B, B + m, B0);
    fill(B0 + m, B0 + p, 0);
    transform(p, B0);
    for(int i = 0; i != p; ++i)
        A0[i] = A0[i] * B0[i] % mod_v;
    inverse_transform(p, A0);
    for(int i = 0; i != m; ++i)
        R[i] = (A[i] - A0[i]) % mod_v;
    fill(R + m, R + p, 0);
}

http://blog.miskcoo.com/2015/05/polynomial-division

应用
多项式的扩展欧几里德算法
多项式的乘法逆元
多项式的多点求值

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值