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

题目大意:同解方程 数据范围m<=10^8

O(m)做法见 http://blog.csdn.net/popoqqq/article/details/40984859

O(m)跪了你就当我没辙么?

首先找到一个比较靠谱的第一个质数 将对第一个质数取模为0的值全都存在一个数组里

由于这个是有循环节的 所以我们只需要处理出[0,p-1]中对第一个质数取模为0的数就可以搞出所有了

然后对于这个数组里的所有数用剩余的质数验证一遍就行了

时间复杂度未知 但是第一个质数必须靠谱 如果第一个质数被卡掉那就退化成O(m)了

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define M 110
using namespace std;
typedef long long ll;
const int prime[]={30011,11261,14843,19997,10007,21893};
int n,m,stack[1001001],top;
int ans[110],tot;
ll a[M][6],f[30011][6];
inline ll F(int x,int j)
{
    int i;
    ll re=0;
    for(i=n;~i;i--)
        re=(re*x+a[i][j])%prime[j];
    return re;
}
inline void Input(int x)
{
    static char s[10100];
    int i,j;
    bool flag=false;
    scanf("%s",s+1);
    for(i=1;s[i];i++)
    {
        if(s[i]=='-')
            flag=true;
        else
            for(j=0;j<6;j++)
                a[x][j]=( (a[x][j]<<1) + (a[x][j]<<3) + s[i]-'0' )%prime[j];
    }
    if(flag)
        for(j=0;j<6;j++)
            a[x][j]=prime[j]-a[x][j];
}
int main()
{

	//freopen("equation.in","r",stdin);

    int i,j;
    cin>>n>>m;
    for(i=0;i<=n;i++)
        Input(i);

    for(j=0;j<6;j++)
        for(i=0;i<prime[j];i++)
            f[i][j]=F(i,j);

    for(i=0;i<prime[0];i++)
    {
    	if(f[i%prime[0]][0]==0)
            stack[++top]=i;
    }

    for(i=1;i<=top;i++)
    {
    	if(stack[i]+prime[0]>m)
    		break;
    	stack[++top]=stack[i]+prime[0];
    }

    for(i=1;i<=top;i++)
    {
    	if(stack[i]>m)
    		break;
    	for(j=1;j<6;j++)
    		if(f[stack[i]%prime[j]][j])
    			break;
    	if(j==6)
    		ans[++tot]=stack[i];
    }

    cout<<tot<<endl;
    for(i=1;i<=tot;i++)
        printf("%d\n",ans[i]);
}
时间: 2024-09-29 20:58:00

Vijos P1915 解方程 加强版 还是Hash大法好!的相关文章

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数据较水,可以直接选择\

[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]内的整数解

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 行每行包含一个整数,依次

HDU 4793 Collision --解方程

题意: 给一个圆盘,圆心为(0,0),半径为Rm, 然后给一个圆形区域,圆心同此圆盘,半径为R(R>Rm),一枚硬币(圆形),圆心为(x,y),半径为r,一定在圆形区域外面,速度向量为(vx,vy),硬币向圆盘撞过去,碰到圆盘后会以相反方向相同速度回来(好像有点违背物理规律啊,但是题目是这样,没办法).问硬币某一部分在圆形区域内的总时间. 解法: 解方程,求 (x+vx*t,y+vy*t) 代入圆形区域方程是否有解,如果没解,说明硬币运动轨迹与圆形区域都不相交,答案为0 如果有解,再看代入圆盘有

【NOIP之旅】NOIP2014 day2 T3 解方程

3.解方程 (equation.cpp/c/pas) [问题描述] 已知多项式方程: 求这个方程在[1, m]内的整数解(n和m均为正整数).   [输入] 输入文件名为equation.in. 输入共n+2行. 第一行包含2个整数n.m,每两个整数之间用一个空格隔开. 接下来的n+1行每行包含一个整数,依次为a0,a1,a2,……,an. [输出] 输出文件名为equation.out. 第一行输出方程在[1, m]内的整数解的个数. 接下来每行一个整数,按照从小到大的顺序依次输出方程在[1,

1743: 解方程

1743: 解方程 Description 一群奥特曼打败了一群小怪兽,已知所有的奥特曼均有x1个头.y1条腿(变异奥特曼),所有的小怪兽均有x2个头.y2条腿.战场上一共有q个头,w条腿,问有多少奥特曼,有多少个小怪兽? Input 输入数据有多组每组包含6个正整数,分别为,x1,y1,x2,y2,q,w :(0<=q,w<=1000000000):输入数据保证有唯一解.读到0 0 0 0 0 0结束. Output 输出占一行,包含两个正整数,分别为:奥特曼和小怪兽的数目: Sample

解方程与高斯消元

解决多元一次方程组: n个方程n个未知数高斯消元法:思路:1.把方程未知数的系数写成一个矩阵2.消去x1,用除第一行以外的其他行减去第一行,写在原位上3.同理第二步,消去第三行的x24.原来的矩阵会变成一个上三角矩阵5.最后直接解方程*如果当前方程需要消元的位置系数是0,就把当前方程和下面第一个非0的方程换一下位置反正也不影响答案 代码实现: void gauss() { juzhen:m for(int i=1;i<=n;++i) { if(m[i][i]==0) { for(int j=i+