[BZOJ 3751][NOIP2014]解方程(哈希)

Posted Zars19

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[BZOJ 3751][NOIP2014]解方程(哈希)相关的知识,希望对你有一定的参考价值。

Description

已知多项式方程:

a0+a1*x+a2*x^2+...+an*x^n=0

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

Solution

一道很久很久以前就应该做的noip的题

一定要放上来是要见证我人品崩坏的一下午

生无可恋…QAQ

题解其实也很简单啦

随便找几个素数取模验证是不是等于0就好了

随便

随便

随便

随便找几个…素数

在WA\\TLE\\OLE间切换,最后还是抄了别人的几个素数

怀疑人生【望天

(BZOJ上的数据是加强了的,如果是ccf的数据那当然就随便找几个素数啦…)

加23333保平安

(震惊…233和23333竟然都是素数)

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
using namespace std;
typedef long long LL;
int prime[5]={839,233,23333,22877,19997},a[105][5];
int n,m,cnt=0,t[1000005];
char s[10005];
int judge[1000005];
int read(char* s,int p)
{
    int x=0,f=1,i=0;
    while(s[i]<\'0\'||s[i]>\'9\'){if(s[i]==\'-\')f=-1;i++;}
    while(s[i]>=\'0\'&&s[i]<=\'9\')x=(x*10+s[i]-\'0\')%p,i++;
    return x*f;
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=0;i<=n;i++)
    {
        scanf("%s",s);
        for(int j=0;j<5;j++)
        a[i][j]=read(s,prime[j]);
    }
    for(int i=0;i<5;i++)
    {
        for(int j=0;j<prime[i]&&j<=m;j++)
        {
            int pow=1;LL ans=0;
            for(int k=0;k<=n;k++)
            {
                ans+=(1LL*pow*a[k][i])%prime[i];
                ans%=prime[i];
                pow=(1LL*pow*j)%prime[i];
            }
            if(!ans)
            {
                for(int k=0;k*prime[i]+j<=m;k++)
                judge[k*prime[i]+j]++;
            }
        }
    }
    for(int i=1;i<=m;i++)
    if(judge[i]==5)t[++cnt]=i;
    printf("%d\\n",cnt);
    for(int i=1;i<=cnt;i++)
    printf("%d\\n",t[i]);
    return 0;
}

 

以上是关于[BZOJ 3751][NOIP2014]解方程(哈希)的主要内容,如果未能解决你的问题,请参考以下文章

[BZOJ3751] [NOIP2014] 解方程 (数学)

[BZOJ 3751][NOIP2014]解方程(哈希)

[BZOJ3751][NOIP2014]解方程(数学相关+乱搞)

hashBZOJ3751-[NOIP2014]解方程

noip2014解方程解题报告

[NOIP2014提高组]解方程