求最大公共多项式

定义1
若h(x)既是f(x)的因式,又是g(x)的因式,则称h(x)是f(x)与g(x)的公因式。 因,c|f(x),c|g(x),并且c!=0,所以任意两个多项式都有公因式。

定义2
设d(x)是f(x)与g(x)的一个公因式,如果对于f(x)与g(x)的 任一个公因式h(x),都有h(x)|d(x)则称d(x)是f(x)与g(x)的一个最大公因式。

定理
如果d(x)是f(x)与g(x)的一个最大公因式,那么它们的所有最大公因式都是形如cd(x)(c属于F,c!=0)的多项式

常见结果
1.g(x)|f(x),g(x)和cg(x)都是f(x)与g(x)的最大公因式
2.0多项式是0多项式与0多项式的最大公因式

代码实现

ll qmi(ll m, ll k, ll mod)
{
    ll res = 1 % mod, t = m;
    while (k)
    {
        if (k&1)
            res = res * t % mod;
        t = t * t % mod;
        k >>= 1;
    }
    return res;
}

/*多项式求最大公共项*/
vector<int> poly_gcd(vector<int> a,vector<int> b)
{
    if (b.size() == 0)
        return a;
    int t = a.size() - b.size();
    vector<int> c;
    for (int i=0; i<=t; i++)
    {
        int tmp=a[i]*qmi(b[0],mod-2,mod)%mod;
        for(ll j=0; j<b.size(); j++)
            a[i+j]=(a[i+j]-tmp*b[j]%mod+mod)%mod;
    }
    int p=-1;
    for(int i=0; i<a.size(); i++)
    {
        if(a[i]!=0)
        {
            p=i;
            break;
        }
    }
    if(p>=0)
    {
        for(int i=p; i<a.size(); i++)
            c.push_back(a[i]);
    }
    return poly_gcd(b,c);
}

int main()
{
    vector<int> ans = poly_gcd(G1[0], G2[1]);
    int cnt = qmi(ans[0], mod-2,mod);
    for (int i = 0; i < ans.size(); i++)
        ans[i] = ans[i] * cnt % mod;
    return 0;
}

  • 0
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值