题目链接

https://www.luogu.org/problemnew/show/P2312

分析

这道题很毒啊,这么大的数。

但是如果多项式\(\sum_{i=0}^N a[i]*X^i=0\)则\(\sum_{i=0}^N a[i]*X^i \mod P=0\)

于是我们可以暴力膜一模,然后在\([1,m]\)中枚举就好了。但是呢,万一这个多项式的值是\(P\)的倍数,也会变成0,所以保险起见搞几个又大又质的数膜一膜就好了。

但是\(Exciting\)的是呢,我在洛谷上开O2能过,而BZOJ就不那么友好。

然后luogu题解提供一种减少枚举冗杂的方Fa。我们不是选多个数膜一模吗,如果在膜\(P_i\)的意义下已经不是\(0\)了,枚举其他的就没意义了。于是呢,我们先可以选出一个小点的模数\(P_x\),在\([1,P_x]\)中先枚举一遍,记录多项式值为0的是哪些。最后再枚举\([1,m]\),由于先前的限制,就会减少许多无用选择

然后多项式求值有个叫秦九韶算法的\(O(N)\)方法,不了解的可以看一看

https://www.cnblogs.com/Rye-Catcher/p/9260599.html

代码

我选择了两个数来做模数,较小的是23333,较大的是19260817

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cctype>
#define ll long long
#define ri register int
using namespace std;
const int maxn=105;
const int inf=0x7fffffff;
const int p1=19260817,p2=71806291,p3=23333;
template <class T>inline void read(T &x){
x=0;int ne=0;char c;
while(!isdigit(c=getchar()))ne=c=='-';
x=c-48;
while(isdigit(c=getchar()))x=(x<<3)+(x<<1)+c-48;
x=ne?-x:x;return ;
}
int n,m;
ll a[maxn],c[maxn];
inline void input(int id){
a[id]=c[id]=0;int ne=0;char ch;
while(!isdigit(ch=getchar()))ne=ch=='-';
a[id]=c[id]=ch-48;
while(isdigit(ch=getchar())){
a[id]=((a[id]<<3)%p1+(a[id]<<1)+ch-48)%p1;
c[id]=((c[id]<<3)%p3+(c[id]<<1)+ch-48)%p3;
}
a[id]=ne?-a[id]:a[id];
c[id]=ne?-c[id]:c[id];return ;
}
int ans[1000005],tot=0;
bool ok[1000005];
inline bool pre_calc(ll u){
ll x=0;
for(ri i=n;i>=0;i--){
x=(x*u+c[i])%p3;
}
return x==0?1:0;
}
inline bool calc(ll u){
ll x=0;
for(ri i=n;i>=0;i--){
x=(x*u+a[i])%p1;
}
return x==0?1:0;
}
int main(){
read(n),read(m);
for(ri i=0;i<=n;i++){
input(i);
}
memset(ok,0,sizeof(ok));
for(ri i=1;i<=p3;i++){
if(pre_calc(1ll*i))ok[i]=1;
}
for(ri i=1;i<=m;i++)
if(ok[i%p3]&&calc(1ll*i)){
ans[++tot]=i;
}
printf("%d\n",tot);
for(ri i=1;i<=tot;i++)printf("%d\n",ans[i]);
return 0;
}

luogu题解P2312解方程--暴力模+秦九韶的更多相关文章

  1. 洛谷P2312 解方程(暴力)

    题意 题目链接 Sol 出这种题会被婊死的吧... 首先不难想到暴力判断,然后发现连读入都是个问题. 对于\(a[i]\)取模之后再判断就行了.注意判断可能会出现误差,可以多找几个模数 #includ ...

  2. 洛谷P2312 解方程题解

    洛谷P2312 解方程题解 题目描述 已知多项式方程: \[a_0+a_1x+a_2x^2+\cdots+a_nx^n=0\] 求这个方程在 \([1,m]\) 内的整数解(\(n\) 和 \(m\) ...

  3. 洛谷 P2312 解方程 题解

    P2312 解方程 题目描述 已知多项式方程: \[a_0+a_1x+a_2x^2+\cdots+a_nx^n=0\] 求这个方程在 [1,m][1,m] 内的整数解(\(n\) 和 \(m\) 均为 ...

  4. bzoj3751 / P2312 解方程

    P2312 解方程 bzoj3751(数据加强) 暴力的一题 数据范围:$\left | a_{i} \right |<=10^{10000}$.连高精都无法解决. 然鹅面对这种题,有一种常规套 ...

  5. [noip2014]P2312 解方程

    P2312 解方程 其实这道题就是求一个1元n次方程在区间[1, m]上的整数解. 我们枚举[1, m]上的所有整数,带进多项式中看看结果是不是0即可. 这里有一个技巧就是秦九韶算法,请读者自行查看学 ...

  6. codevs3732==洛谷 解方程P2312 解方程

    P2312 解方程 195通过 1.6K提交 题目提供者该用户不存在 标签数论(数学相关)高精2014NOIp提高组 难度提高+/省选- 提交该题 讨论 题解 记录   题目描述 已知多项式方程: a ...

  7. 洛谷 P2312 解方程 解题报告

    P2312 解方程 题目描述 已知多项式方程: \(a_0+a_1x+a_2x^2+\cdots+a_nx^n=0\)求这个方程在 \([1,m]\) 内的整数解(\(n\) 和 \(m\) 均为正整 ...

  8. P2312 解方程(随机化)

    P2312 解方程 随机化的通俗解释:当无法得出100%正确的答案时,考虑随机化一波,于是这份代码很大可能会对(几乎不可能出错). 比如这题:把系数都模一个大质数(也可以随机一个质数),然后O(m)跑 ...

  9. 洛谷 P2312 & bzoj 3751 解方程 —— 取模

    题目:https://www.luogu.org/problemnew/show/P2312 https://www.lydsy.com/JudgeOnline/problem.php?id=3751 ...

随机推荐

  1. POJ 2109 -- Power of Cryptography

    Power of Cryptography Time Limit: 1000MS   Memory Limit: 30000K Total Submissions: 26622   Accepted: ...

  2. 【Linux命令】find命令

    [find命令] 说明:find命令用来在指定目录下查找文件.任何位于参数之前的字符串都将被视为欲查找的目录名.如果使用该命令时,不设置任何参数,则find命令将在当前目录下查找子目录与文件.并且将查 ...

  3. 【Makefile】Makefile中的赋值符号=、:=、?=、+=

    Makefile中主要有四个赋值符号: = 是最基本的赋值:= 是覆盖之前的值?= 是如果没有被赋值过就赋予等号后面的值+= 是添加等号后面的值 1.“=” “=”赋值:make会将整个makefil ...

  4. LC 835. Image Overlap

    Two images A and B are given, represented as binary, square matrices of the same size.  (A binary ma ...

  5. Python数据分析中 DataFrame axis=0与axis=1的理解

    python中的axis究竟是如何定义的呢?他们究竟代表是DataFrame的行还是列? 直接上代码people=DataFrame(np.random.randn(5,5), columns=['a ...

  6. flutter webview_flutter 设置cookies

    WebViewController _controller; onWebViewCreated: (WebViewController wvc) { _controller = wvc; } onPa ...

  7. 1、DNS服务基础

    w'indows上名称解析目录: C:\Windows\System32\drivers\etc https://jocent.me/2017/06/18/dns-protocol-principle ...

  8. docker版的zabbix部署

    环境准备:一台server端,两台agent端 server端部署zabbix-server和mariadb服务 agent端部署zabbix-agent服务 一.docker容器里下载zabbix和 ...

  9. vue如何获取并操作DOM元素

    原文地址 方法一: 直接给相应的元素加id,然后再document.getElementById("id");获取,然后设置相应属性或样式 方法二: 使用ref,给相应的元素加re ...

  10. Vue 使用eventBus 实现兄弟组件间的通信

    实现方式:  主要是在相互通信的兄弟组件之中,都引入一个新的vue实例,然后通过分别调用这个实例的事件触发事件广播 和监听来实现通信和参数传递. 需求: a页面tree的增删改后,数据还是之前的老数据 ...