NOIP2014/洛谷P2312 解方程

Posted 沐灵_hh

tags:

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

NOIP2014/洛谷P2312 解方程

题目描述

已知多项式方程:

a0+a1x+a2x^2+..+anx^n=0

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

输入输出格式

输入格式:

输入文件名为equation .in。

输入共n + 2 行。

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

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

输出格式:

输出文件名为equation .out 。

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

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

输入输出样例

输入样例#1: 复制
2 10 
1
-2
1
输出样例#1: 复制
1
1
输入样例#2: 复制
2 10
2
-3
1
输出样例#2: 复制
2
1
2
输入样例#3: 复制
2 10 
1  
3  
2  
 
输出样例#3: 复制
0

说明

对于30%的数据:0<n<=2,|ai|<=100,an!=0,m<100

对于50%的数据:0<n<=100,|ai|<=10^100,an!=0,m<100

对于70%的数据:0<n<=100,|ai|<=10^10000,an!=0,m<10000

对于100%的数据:0<n<=100,|ai|<=10^10000,an!=0,m<1000000

分析

对于30%的数据,用公式解就可以了。

对于50%的数据,高精度?

对于100%,m的范围不大,在这个范围内枚举,计算最后的值是否为0。要用到秦九韶算法,这样只需要做n次乘法,n次加法。但是a很大,高精度也很难表示怎么办?如果方程左面的值为0,那么方程两边同时mod一个值还是成立的。如果只mod一个值的话,答案很有可能不是正确的,因为如果这个值要是mod这个值的倍数的话,mod的值也是0。我们可以多mod几个值,那么这个数就必须满足是这几个值的最小公倍数的倍数才能在mod之后都为0,否则的话这个数就是0。但是前一种情况机率还是很小的,因为如果模数都取质数的话最小公倍数就是这几个数的乘积,还是很大的。但是这样的话时间复杂度是O(n*m)的,会超时的。如果一个数x mod p !=0 的话,那么x+p mod p!=0,这样就可以少计算一些无用的数。这样的话应该能通过官方数据。

关于秦九韶算法a[0]是否要乘以x,从公式上来看是不需要的,但是如果到最后结果如果是0的话就算乘以x之后还是0,对结果没有什么太大的影响,而且好记,代码短(lazy!)。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=1000010;
const int mo1=99991;
const int mo2=998244353;
const int mo3=1e9+7;
int n,m,cnt,a[N],b[N],c[N],ans[N]; bool vis[N];
inline void read(int i){
    ll x1=0,x2=0,x3=0,f=1; char ch=getchar();
    while(ch<0||ch>9){if(ch==-)f=-1;ch=getchar();}
    while(ch>=0&&ch<=9){
        x1=(x1*10+ch-0)%mo1;
        x2=(x2*10+ch-0)%mo2;
        x3=(x3*10+ch-0)%mo3;
        ch=getchar();
    }
    a[i]=x1*f;b[i]=x2*f;c[i]=x3*f;
}
inline bool calc1(int x){
    ll sum=0;
    for(int i=n;i>=0;--i)
        sum=((a[i]+sum)*x)%mo1;
    return !sum;
}
inline bool calc2(int x){
    ll sum1=0,sum2=0;
    for(int i=n;i>=0;--i){
        sum1=((b[i]+sum1)*x)%mo2;
        sum2=((c[i]+sum2)*x)%mo3;
    }
    return (!sum1)&&(!sum2);
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=0;i<=n;++i) read(i);
    for(int i=1;i<=mo1;++i)
    if(calc1(i)){
        for(int j=i;j<=m;j+=mo1)
        if(calc2(j)) vis[j]=1;
    }
    for(int i=1;i<=m;++i) if(vis[i]) ans[++cnt]=i;
    printf("%d\n",cnt);
    for(int i=1;i<=cnt;++i) printf("%d\n",ans[i]);
    return 0;
}

 

以上是关于NOIP2014/洛谷P2312 解方程的主要内容,如果未能解决你的问题,请参考以下文章

[noip2014]P2312 解方程

[NOIP2014提高组]解方程

洛谷 P2312 解方程

NOIP2014提高组解方程

P2312 解方程

P2312 解方程