【NOIP2014】解方程 玄学哈希+乱搞

【NOIP2014】解方程 玄学哈希+乱搞_第1张图片
【NOIP2014】解方程 玄学哈希+乱搞_第2张图片

我觉得我考的话最多敲个高精了……

n=0 => n%p=0,显然是对的。

n%p=0 => n=0,显然是错的。

但p多选几个,就极有可能是对的……

所以说我们可以对输入给定的数模好几个质数,然后枚举1~m看看有几个能使方程成立,这样能拿70 。

满分做法:n%p=(n+kp)%p,所以每次只枚举[0,p-1],然后把倍数筛掉即可。

p不能选太大,否则会T。p不能选太小,否则会WA。p不能选太少,否则筛不干净……

总之我这样的脸黑的,做这种题简直丧病……

还有,这个题实测加不加秦九韶好像差不多

代码:

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int SIZE=2000010;
typedef long long LL;

LL p[]={11279,11953,12133,12197,13831};

LL a[SIZE][10];


void add(LL &x,LL mod,char c)
{
    x=(((x<<3)%mod+(x<<1)%mod)%mod+c-'0')%mod;
}

void change(LL &x,LL mod)
{
    x=((-x)%mod+mod)%mod;
}
void scan(LL i)
{
    char c=getchar();
    bool flag=0;
    while(c<'0'||c>'9') {if(c=='-') flag=1;c=getchar();}
    while(c>='0'&&c<='9')
    {
        for(int j=0;j<5;j++)
            add(a[i][j],p[j],c);
        c=getchar();
    }
    if(flag)
    {
        for(int j=0;j<5;j++)
            change(a[i][j],p[j]);   
    }
}
int n,m;
LL f(LL x,LL id)
{
    LL mod=p[id];
    LL ans=a[n][id];
    for(int i=n;i>=1;i--)
    {
        ans=( ( ( (ans%mod)* (x%mod) )%mod +a[i-1][id]%mod) %mod)%mod;
    }
    return ans;
/* LL mod=p[id],ans=0,xx=1; for(int i=0;i<=n;i++) { ans=(ans+((a[i][id]%mod)*(xx%mod))%mod)%mod; xx=((xx%mod)*(x%mod))%mod; } return ans;*/
}
bool use[SIZE];
int ans[SIZE];
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=0;i<=n;i++)
        scan(i);

    for(int j=0;j<5;j++)
    {
        for(int i=0;i<p[j];i++)
        {
            if(f(i,j)!=0)
            {
                for(int k=i;k<=m;k+=p[j]) 
                    use[k]=1;
            } 
        } 
    } 

    for(int i=1;i<=m;i++)
    {
        if(!use[i])
            ans[++ans[0]] = i;  
    }
    for(int i=0;i<=ans[0];i++)
    {
        printf("%d\n",ans[i]);
    }
    return 0;
}


你可能感兴趣的:(noip)