zoukankan      html  css  js  c++  java
  • [NOIp提高组2014]解方程

    题目大意

    求方程

    [sum_{i=0}^{n}a_ix^i=0 ]

    ([1,m])内的整数解

    (1 leq |a_i| leq 10^{10000},a_i eq 0,1 leq n leq 100,1 leq m<10^6)

    解题思路

    最朴素的做法就是尝试所有解,判断左边多项式值是否为零

    但还有一个高精问题(如果你真要用高精我也不拦你)

    可以考虑一种类哈希做法,如果所有运算都对一个大质数取模,答案的正确性基本可以保证,而且可以省掉高精

    实测计算左侧多项式秦九昭算法和直接算的时间相差不大(只要你不是用快速幂算的)

    (O(nm))做法,貌似(10^8)再加一堆*%运算,卡常?

    没错,我在Calc函数内写ret*=x;ret+=X[i];ret%=P就被卡掉QAQ

    #include<iostream>
    #include<cstdio>
    #include<vector>
    
    const long long P=998244353;
    
    int n,m;
    
    long long X[200];
    
    int all;
    std::vector<int> V;
    
    void read(long long &x){
    	char ch;
    	while (isspace(ch=getchar()));
    	bool flag=false;
    	if (ch=='-'){flag=true;ch=getchar();}
    	x=ch&15;
    	while (isdigit(ch=getchar())){x=(x<<1)+(x<<3)+(ch&15);x%=P;}
    	if (flag) x=-x+P;
    }
    
    long long Calc(int x){
    	long long ret=X[n];
    	for (int i=n-1;i>=0;i--) ret=(ret*x+X[i])%P;
    	return ret;
    }
    
    int main(){
    	scanf("%d%d",&n,&m);
    	for (int i=0;i<=n;i++) read(X[i]);
    	for (int x=1;x<=m;x++){
    		if (!Calc(x)){
    			all++;
    			V.push_back(x);
    		}
    	}
    	printf("%d
    ",all);
    	for (std::vector<int>::iterator it=V.begin();it!=V.end();it++) printf("%d
    ",*it);
    }
    
  • 相关阅读:
    Count the Buildings HDU
    Airport UVA
    17南宁网络赛
    git Please move or remove them before you can merge.
    php ajax bootstrap多文件上传图片预览,ajax上传文件
    thinkphp5多语言
    ueeditor 百度编译器使用onchange效果
    mysql sql_mode=only_full_group_by问题?
    mac navicate 2013
    mac Nginx+CI出现404错误
  • 原文地址:https://www.cnblogs.com/ytxytx/p/9499662.html
Copyright © 2011-2022 走看看