解方程(hash,秦九韶算法)

题目描述

已知多项式方程:

a0+a1x+a2x2+⋯+anxn=0

求这个方程在 [1,m]内的整数解(nm 均为正整数)。

输入输出格式

输入格式:

n+2 行。

第一行包含 2个整数 n,m,每两个整数之间用一个空格隔开。

接下来的 n+1n+1n+1 行每行包含一个整数,依次为 a0,a1,a2…an

输出格式:

第一行输出方程在 [1,m] 内的整数解的个数。

接下来每行一个整数,按照从小到大的顺序依次输出方程在 [1,m]内的一个整数解。

思路:

有一个奇妙的东西叫秦九韶算法

(不知道的请左转百度)

然后我们看到a都很大

于是我们想到了hash

将 a hash掉

之后验证时取膜即可

代码:

// luogu-judger-enable-o2
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define rii register int i
#define rij register int j
#define p1 19260817
#define p2 23333
using namespace std;
char ai[10005];
long long a1[105],a2[105];
int ans[1000005],n,m,cnt;
inline void c1(int wz,int len)
{
    int pd=1;
    int fi=0;
    if(ai[0]=='-')
    {
        pd=-1;
        fi=1;
    }
    long long ans=ai[fi]-'0';
    for(rii=fi+1;i<=len-1;i++)
    {
        ans*=10;
        ans+=ai[i]-'0';
        ans%=p1;
    }
    a1[wz]=ans*pd;
}
inline void c2(int wz,int len)
{
    int pd=1;
    int fi=0;
    if(ai[0]=='-')
    {
        pd=-1;
        fi=1;
    }
    long long ans=ai[fi]-'0';
    for(rii=fi+1;i<=len-1;i++)
    {
        ans*=10;
        ans+=ai[i]-'0';
        ans%=p2;
    }
    a2[wz]=ans*pd;;
}
inline bool check(int v)
{
    long long ans=a1[n];
    for(rii=n;i>=1;i--)
    {
        ans*=v;
        ans%=p1;
        ans+=a1[i-1];
        ans%=p1;
    }
    if(ans!=0)
    {
        return 0;
    }
    ans=a2[n];
    for(rii=n;i>=1;i--)
    {
        ans*=v;
        ans%=p2;
        ans+=a2[i-1];
        ans%=p2;
    }
    if(ans!=0)
    {
        return 0;
    }
    return 1;
}
inline void sr()
{
    int wz=0;
    char l=getchar();
    while(l!=10)
    {
        ai[wz]=l;
        l=getchar();
        wz++;
    }
}
int main()
{
    scanf("%d%d\n",&n,&m);
    for(rii=0;i<=n;i++)
    {
        memset(ai,0,sizeof(ai));
//        scanf("%s",ai);
        sr();
        int ltt=strlen(ai);
        c1(i,ltt);
        c2(i,ltt);
    }
//    for(rii=0;i<=n;i++)
//    {
//        cout<<a2[i]<<" ";
//    }
    for(rii=1;i<=m;i++)
    {
        if(check(i)==1)
        {
            cnt++;
            ans[cnt]=i;
        }
    }
    cout<<cnt<<endl;
    for(rii=1;i<=cnt;i++)
    {
        printf("%d\n",ans[i]);
    }
}

 

转载于:https://www.cnblogs.com/ztz11/p/9636785.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值