luogu2312 解方程 (数论,hash)

luogu2312 解方程 (数论,hash)

第一次外出学习讲过的题目,然后被讲课人的一番话惊呆了.
这个题,我想着当年全国只有十几个满分.....然后他又说了句我考场A这道题时,用了5个模数
确实不好做想不到.
由于\(a\)非常大.转为以下思路.

\(f(x) = a_0+a_1x+a_2x^2+\cdots+a_nx^n\)
对于\(f(x) = 0\)则\(f(x)\%p = 0\)
\(f(x \% p) = 0\)
然后这里最好是选择素数.由于luogu数据较水,可以直接选择\(1e9 + 7\)水过.
判断会有误,所以这里选择两个数.
bzoj也有这道题,但是极其考验卡常技巧......
卡到吐血.还是没A,算了,不卡了....
记录:

#pragma GCC optimize(2)
#include <iostream>
#include <cstdio>
#define ll long long
const int maxN = 100 + 7;
const int maxM = 1e6 + 7;
const int p = 1e9 + 7;
const int p1 = 20030327;

ll a[maxN],b[maxN];
ll n,m;
bool vis[maxM];

inline bool calc(int x)
{
    long long sum = 0;
    for(int i = n;i >= 1;i --)
    {
        sum = ( (long long) ( a[i] + sum ) * x ) % p;
    }
    sum = ( sum + a[0] ) % p;
    return !sum;
}

inline void read(ll &x1,ll &x2) {
    x1 = 0,x2 = 0;
    int f = 1;
    char c = getchar();
    while(c < '0' || c > '9') {if(c == '-')f = -1;c = getchar();}
    while(c >= '0' && c <= '9') {
        x1 = ( x1 * 10 + c - '0' ) % p;
        x2 = ( x2 * 10 + c - '0') % p1;
        c = getchar();
    }
    x1 *= f;
    x2 *= f;
}

inline bool calc1(int x)
{
    long long sum = 0;
    for(int i = n;i >= 1;i --)
    {
        sum = ( (long long) ( b[i] + sum ) * x ) % p1;
    }
    sum = ( sum + b[0] ) % p1;
    return !sum;
}

void print(int x)
{
    if(x < 0)
    {
        putchar('-');
        x = -x;
    }
    if(x > 9)
    {
        print(x / 10);
    }
    putchar(x % 10 + '0');
}

int main() {
    scanf("%d%d",&n,&m);
    for(int i = 0;i <= n;++ i)
        read(a[i],b[i]);
    int cnt = 0;
    for(int i = 1;i <= m;++ i)
        if(calc(i) && calc1(i)) ++ cnt,vis[i] = true;
    printf("%d\n",cnt);
    for(int i = 1;i <= m;++ i)
        if(vis[i]) print(i),puts("");
    return 0;
}

原文地址:https://www.cnblogs.com/tpgzy/p/9715362.html

时间: 2024-10-08 11:13:14

luogu2312 解方程 (数论,hash)的相关文章

【bzoj3751】[NOIP2014]解方程 数论

题目描述 已知多项式方程: a0+a1*x+a2*x^2+...+an*x^n=0 求这个方程在[1,m]内的整数解(n和m均为正整数). 输入 第一行包含2个整数n.m,每两个整数之间用一个空格隔开. 接下来的n+1行每行包含一个整数,依次为a0,a1,a2,...,an. 输出 第一行输出方程在[1,m]内的整数解的个数. 接下来每行一个整数,按照从小到大的顺序依次输出方程在[1,m]内的一个整数解. 样例输入 2 10 2 -3 1 样例输出 2 1 2 题解 真心不难的数论题 首先高精度

luogu2312 解方程

题目大意 已知多项式方程:a0+a1x+a2x^2+..+anx^n=0求这个方程在[1, m ] 内的整数解(n 和m 均为正整数)ai<=10^10000 题解 枚举多个不太大的质数$p_i$,枚举$x\in[0,p_i-1]$,预处理,用秦九韶算法看看$f(x)$是否为0(f运算时取模).多枚举几个质数进行上述操作,对于$x\in [1,m]$,若对$\forall p_i, f(x\mathrm{mod}p_i)=0$,则$x$为一个解. 注意.如果枚举的$p_i$是个大质数,枚举$x\

BZOJ3751 NOIP2014 解方程(Hash)

题目链接  BZOJ3751 这道题的关键就是选取取模的质数. 我选了4个大概几万的质数,这样刚好不会T 然后统计答案的时候如果对于当前质数,产生了一个解. 那么对于那些对这个质数取模结果为这个数的数也要统计进答案. #include <bits/stdc++.h> using namespace std; #define rep(i, a, b) for (int i(a); i <= (b); ++i) #define dec(i, a, b) for (int i(a); i &g

$Noip2014/Luogu2312$ 解方程

$Luogu$ $Sol$ 枚举解+秦九韶公式计算+取模. $Code$ #include<iostream> #include<cstdio> #include<cstring> #include<queue> #include<algorithm> #define il inline #define Rg register #define go(i,a,b) for(Rg int i=a;i<=b;++i) #define yes(i,

Vijos P1915 解方程 加强版 还是Hash大法好!

题目大意:同解方程 数据范围m<=10^8 O(m)做法见 http://blog.csdn.net/popoqqq/article/details/40984859 O(m)跪了你就当我没辙么? 首先找到一个比较靠谱的第一个质数 将对第一个质数取模为0的值全都存在一个数组里 由于这个是有循环节的 所以我们只需要处理出[0,p-1]中对第一个质数取模为0的数就可以搞出所有了 然后对于这个数组里的所有数用剩余的质数验证一遍就行了 时间复杂度未知 但是第一个质数必须靠谱 如果第一个质数被卡掉那就退化

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

P2312 解方程 195通过 1.6K提交 题目提供者该用户不存在 标签数论(数学相关)高精2014NOIp提高组 难度提高+/省选- 提交该题 讨论 题解 记录 题目描述 已知多项式方程: a0+a1x+a2x^2+..+anx^n=0 求这个方程在[1, m ] 内的整数解(n 和m 均为正整数) 输入输出格式 输入格式: 输入文件名为equation .in. 输入共n + 2 行. 第一行包含2 个整数n .m ,每两个整数之间用一个空格隔开. 接下来的n+1 行每行包含一个整数,依次

[noip2014]P2312 解方程

P2312 解方程 其实这道题就是求一个1元n次方程在区间[1, m]上的整数解. 我们枚举[1, m]上的所有整数,带进多项式中看看结果是不是0即可. 这里有一个技巧就是秦九韶算法,请读者自行查看学习. 时间复杂度O(n*m). 然后你应该可以拿30分. 我们发现这些数都太大了,要开高精度.然后你愉快地拿了50分——复杂度O(n*m*length)会爆炸. 这里我们考虑hash的思想,对结果取模(最好是一个很大的质数P),如果结果是零就说明这是一个解. 应为如果结果是零,那么要么这是一个解,要

用python解方程和微积分

用python解方程: from sympy import * x = Symbol('x')  y = Symbol('y') print solve([2* x - y -3,3* x + y -7],[x, y]) 2. 求极限: 代码中的oo就代表无穷. from sympy import * n = Symbol('n') s = ((n+3)/(n+2))**n print limit(s, n, oo) 3. 求定积分: integrate函数用于积分问题. from sympy 

NOIP201410解方程(C++)

NOIP201410解方程 难度级别:A: 运行时间限制:1000ms: 运行空间限制:51200KB: 代码长度限制:2000000B 试题描述 已知多项式方程: a0+a1*x+a2*x^2+a3*x^3+-+an*x^n=0 求这个方程在[1, m]内的整数解(n 和 m 均为正整数). 输入 输入共 n+2 行.第一行包含 2 个整数 n.m,每两个整数之间用一个空格隔开.接下来的 n+1 行每行包含一个整数,依次为a0,a1,a2-an.  输出 第一行输出方程在[1, m]内的整数解