codevs 3732 NOIP2014 解方程 解题报告

题目描述 Description

输入描述 Input Description

输入文件名为equation.in。

输入共n+2行。

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

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

输出描述 Output Description

输出文件名为equation.out。

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

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

样例输入 Sample Input

equation.in

2 10

1

-2

1

equation.out

1

1
/

equation.in

2 10

2

-3

1

equation.out

2

1

2
/

equation.in

2 10

1

3

2

equation.out

0
//
数据范围及提示 Data Size & Hint

思路

玄学

代码

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<vector>
using namespace std;
const int M=1e6+5;
const int N=1e4+5;
int b[M],n,m,ansn,a[M],ans[M];
int f[105][3],p[3]={10007,10917,30071};
char str[N];
int modit(int x)
{
    int flag=0,len=strlen(str);
    if (str[0]=='-')flag=1,a[len]=0;
    a[0]=0;
    for (int i=flag;i<=len-1;i++)
    a[len-i]=str[i]-48;
    for (int i=len;i>=1;i--)
    a[i-1]+=a[i]%x*10;
    if (!flag) return a[0]/10;
    else return (x-a[0]/10)%x;
}
int judge(int x)
{
    for (int i=0;i<=2;i++)
    {
        long long s=f[n][i];
        for (int j=n-1;j>=0;j--)
        s=(s*x%p[i]+f[j][i])%p[i];
        if (s)
        {
            while(x<m) {x+=p[i];b[x]=1;}
            return 0;
        }
    }
    return 1;
}
int main()
{
    scanf("%d%d",&n,&m);
    for (int i=0;i<=n;i++)
    {
        scanf("%s",str);
        for (int j=0;j<=2;j++)
        f[i][j]=modit(p[j]);
    }
    for (int i=1;i<=m;i++)
    if (!b[i]) if(judge(i)) ans[++ansn]=i;
    printf("%d\n",ansn);
    for (int i=1;i<=ansn;i++)
    printf("%d\n",ans[i]);
}