NOIp提高组2014 解方程————数论综合

本文探讨了一元多次方程的求解方法,利用秦九韶定理将高次方程转换为一系列一次方程,进而求解在指定区间内的整数解。通过模大质数的方法处理大数据输入,避免整数溢出问题,并提供了高效的代码实现。

题解:本题主要考查数论综合。
简要题意:a0​+a1​x+a2​x2+⋯+an​xn=0a_0​+a_1​x+a_2​x^2+⋯+a_n​x^n=0a0+a1x+a2x2++anxn=0求这个方程在[1,m][1,m][1,m]内的整数解(nnnmmm均为正整数)。
1.秦九韶定理:解一元多次方程定理:a0​+a1​x+a2​x2+⋯+an​xna_0​+a_1​x+a_2​x^2+⋯+a_n​x^na0+a1x+a2x2++anxn
=a0​+x∗(a1+a2x+a3x2+⋯+anxn−1)=a_0​+x*(a_1+a_2x+a_3x^2+⋯+a_nx^{n-1})=a0+x(a1+a2x+a3x2++anxn1)
=a0​​+x∗(a1+x∗(a2+a3x+⋯+anxn−2))=a_0​​+x*(a_1+x*(a_2+a_3x+⋯+a_nx^{n-2}))=a0+x(a1+x(a2+a3x++anxn2))
=⋯=⋯=
=a0​​+x∗(a1+⋯+x∗(an−2+x∗(an−1+anx+)))=a_0​​+x*(a_1+⋯+x*(a_{n-2}+x*(a_{n-1}+a_nx+) ) )=a0+x(a1++x(an2+x(an1+anx+)))
这样,求nnn次多项式f(x)f(x)f(x)的值就转化为求n个一次多项式的值。
2.一系列的优化:100%100\%100% 的数据,0<n≤100,∣ai∣≤1010000,m≤1060<n≤100,|a_i|≤10^{10000},m≤10^60<n100,ai1010000,m106,数据异常的大,已经爆longlonglong longlonglong了,所以防爆intintint常用方法:模大质数!质数一定要大,最好多模几个大质数。要写读写优化,否则超时
代码如下:

#include<cstdio>
#include<iostream>
#define mod 1000000009
long long n,m,sum;
long long a[105],Ans[1000005];
long long Read()
{
    char ch;int fh=1;
    scanf("%c",&ch);
    while (ch<'0'||ch>'9')
    {
        if(ch=='-') fh=-1;
        scanf("%c",&ch);
    }
    long long num=0;
    while (ch>='0'&&ch<='9')
    {
        num=num*10+ch-'0';
        num%=mod;
        scanf("%c",&ch);
    }
    num*=fh;
    return num;
}
int main()
{
    scanf("%lld%lld",&n,&m);
    for(long long i=0;i<=n;i++)
    a[i]=Read();
    for(long long x=1;x<=m;x++)
    {
        long long ans=0;
        for(long long i=n;i>=1;i--)
        ans=((ans+a[i])*x)%mod;
        if((ans+a[0])%mod==0)
        Ans[++sum]=x;
    }
    printf("%lld\n",sum);
    for(long long i=1;i<=sum;i++)
    printf("%lld\n",Ans[i]);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值