hdu2769:枚举+同余方程
2024-10-15 22:49:10
题意:有一个随机数生成器 x[i+1]=(a*x[i]+b)%10001
已知 x1,x3,x5...求 x2,x4,x6......
x的个数为 2n (n<=10000) a,b也在 0到10000之间
分析,有 a,b两个未知数,不好解方程,只能通过枚举
比赛的时候想了一下枚举a,b感觉复杂度略大 没敢写(后来据说姿势优美的暴力也能过)
正解是枚举 a,求b
带入x1,x3即可得到一个很简单的exgcd同余方程,对0,10000的每一个a,解出一个b进行验证即可
代码:
#include <iostream>
#include <stdio.h>
#include<string.h>
#include<algorithm>
#include<string>
#include<ctype.h>
using namespace std;
#define mod 10001
long long in[];
long long out[];
long long exgcd(long long a,long long b,long long &x,long long &y)
{
if(!b)
{
x=;
y=;
return a;
}
long long tt=exgcd(b,a%b,x,y);
long long t;
t=x;
x=y;
y=(t-a/b*y);
return tt; }
int main()
{
int n;
scanf("%d",&n);
for(int i=;i<n;i++)
scanf("%I64d",in+i);
int ok=;
for(long long i=;i<=;i++)
{
long long c=((in[]-i*i*in[]%)%+)%;
long long a=i+;
long long b=-;
long long x,y;
long long d=exgcd(a,b,x,y);
if(c%d)
continue;
long long t=abs(b/d);
x*=c/d;
x=(x%t+t)%t;
for(int j=;j<n;j++)
{
out[j]=(i*in[j]+x)%;
if(j==n-)
{
ok=;
break;
}
if(((i*out[j]+x)%)!=in[j+])
break;
}
if(ok)
break;
}
if(!ok)
{
puts("bug");
}
for(int i=;i<n;i++)
{
printf("%I64d\n",out[i]);
}
return ;
}
最新文章
- mysql数据库引擎 MyISAM和 InnoDB区别
- 最小/大费用最大流模板(codevs1914)
- 20145223《Java程序程序设计》第1周学习总结
- 夺命雷公狗ThinkPHP项目之----企业网站22之网站前台中间层(解决代码冗余)
- win7下环境搭建
- C# 中对WinForm窗体中的控件快速设置TableIndex次序
- 简单的C#线程开发实例(隔一秒改变一下Label的Text)
- Linux tar 解压缩命令
- 用HiveDB横向切分MySQL数据库
- Git使用操作指南和GitHub
- Docker基于已有的镜像制新的镜像
- springBoot系列教程03:redis的集成及使用
- HTML中的Hack条件注释语句
- MATLAB实现二值化函数
- [福大软工] Z班——个人技术博客评分
- 剑指offer(61)序列化二叉树
- 04-接口隔离原则(ISP)
- HashMap中的hash算法总结
- 【叶问】 MySQL常用的sql调优手段或工具有哪些
- debezium mongodb 集成测试