洛谷 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


不想做题,看到这题,总想暴力。
还要加高精度,想想就觉得麻烦,苟且地看了题解。
首先给hzwer一个赞。
就是去几个质数,也不用太大,类似于等价表达式那道题,不过更简单。

一下代码转载自hzwer
#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int N=105;
int n,m,a[11][N],pre[11][N],ans[1000005],res[11][30005],mod[]={0,11261,19997,22877,21893,14843};
char ch[10005];
int calc(int t,int x)
{
	int sum=0;
	for(int i=0;i<=n;i++)
		sum=(sum+a[t][i]*pre[t][i])%mod[t];
	if(sum<0)
		sum+=mod[t];
	return sum;
}
bool jud(int x)
{
	for(int i=1;i<=5;i++)
		if(res[i][x%mod[i]]!=0)
			return 0;
	return 1;
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=0;i<=n;i++)
	{
		scanf("%s",ch+1);
		int l=strlen(ch+1);
		bool flg=0;
		for(int j=1;j<=5;j++)
			if(ch[1]!='-')
				a[j][i]=ch[1]-'0';
			else
				a[j][i]=0,flg=1;
		for(int j=1;j<=5;j++)
		{
			for(int k=2;k<=l;k++)
				a[j][i]=(a[j][i]*10+ch[k]-'0')%mod[j];
			if(flg)
				a[j][i]*=-1;
		}
	}
	for(int i=1;i<=5;i++)
		for(int j=1;j<=mod[i];j++)
		{
			pre[i][0]=1;
			for(int k=1;k<=n;k++)
				pre[i][k]=(pre[i][k-1]*j)%mod[i];
			res[i][j]=calc(i,j);
		}
	for(int i=1;i<=m;i++)
		if(jud(i))
			ans[++ans[0]]=i;
	printf("%d\n",ans[0]);
	for(int i=1;i<=ans[0];i++)
		printf("%d\n",ans[i]);
	return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值