Password
Description
Rivest是密码学专家。近日他正在研究一种数列E = {E[1],E[2],……,E[n]},
且E[1] = E[2] = p(p为一个质数),E[i] = E[i-2]*E[i-1] (若2<i<=n)。
例如{2,2,4,8,32,256,8192,……}就是p = 2的数列。在此基础上他又设计了一种加密算法,该算法可以通过一个密钥q (q < p)将一个正整数n加密成另外一个正整数d,计算公式为:d = E[n] mod q。现在Rivest想对一组数据进行加密,但他对程序设计不太感兴趣,请你帮助他设计一个数据加密程序。
Input
第一行读入m,p。其中m表示数据个数,p用来生成数列E。 以下有m行,每行有2个整数n,q。n为待加密数据,q为密钥。 数据范围: 0 < p n< 2^31 0 < q < p 0 < m <= 5000。
Output
将加密后的数据按顺序输出到文件 第i行输出第i个加密后的数据。 输入样例1 2 7 4 5 4 6 输入样例2 4 7 2 4 7 1 6 5 9 3
Sample Input
输入样例1
2 7
4 5
4 6
输入样例2
4 7
2 4
7 1
6 5
9 3
Sample Output
输出样例1
3
1
输出样例2
3
0
1
1
时空限制
1s,64MB
【题解】
我心中的题目难度排名:数学>玄学>信息学……千古难题数学题,考试的时候刚开始连式子都没推,直接暴力递推骗来49分。后来一直在打大模拟,最后临交卷回来看一眼好像是斐波那契,起码应该拿矩阵快速幂和普通快速幂搞一搞,但是没有时间了就没有打。
欧拉定理,我的理解就是用来降幂,不过这道题好像模了phi[c]之后并没有加phi[c],此等天机本蒟蒻不明觉厉。这个式子推一推就发现p的指数是斐波那契数列,但是斐波那契数列增长很快,所以只能用欧拉定理降幂。具体做法是预处理或直接求q的欧拉函数,把它作为求斐波那契数列的矩阵快速幂里的模数,然后得到降幂后的指数再普通快速幂求p的幂。
结合快速幂和欧拉定理,这道题大概并不太难。但是在考场上一是没有仔细想,二来也根本没有降幂的意识。数论就像英语音标一样,学了一遍又一遍,还是忘的多会的少。改题的过程更痛苦,几乎是重学欧拉函数(欧拉定理我学过吗?),原来学的数论只为了做那么两三道题,过后连点定义都说不明白。提高数学水平,是迫在眉睫的问题了。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#define ll long long
using namespace std;
const int sj=;
int m;
ll n,p,q,temp,k,phi[sj],s[sj],ge;
bool v[sj];
void prime(ll x)
{
for(ll i=;i<x;i++)
{
if(!v[i])
{
s[ge++]=i;
phi[i]=i-;
}
for(ll j=;j<ge&&i*s[j]<x;j++)
{
ll mb=i*s[j];
v[mb]=;
if(i%s[j]==)
{
phi[mb]=phi[i]*s[j];
break;
}
else phi[mb]=phi[i]*(s[j]-);
}
}
phi[]=;
}
ll fb(ll x,ll y)
{
ll a[][]={},ans[][]={},f[][]={};
a[][]=a[][]=a[][]=ans[][]=ans[][]=;
while(x)
{
if(x&)
{
memset(f,,sizeof(f));
for(int i=;i<;i++)
for(int j=;j<;j++)
for(int l=;l<;l++)
f[i][j]+=ans[i][l]*a[l][j]%y;
memcpy(ans,f,sizeof(f));
}
x>>=;
memset(f,,sizeof(f));
for(int i=;i<;i++)
for(int j=;j<;j++)
for(int l=;l<;l++)
f[i][j]+=a[i][l]*a[l][j]%y;
memcpy(f,a,sizeof(a));
}
memset(f,,sizeof(f));
memset(a,,sizeof(a));
f[][]=;
f[][]=;
for(int i=;i<;i++)
for(int j=;j<;j++)
for(int l=;l<;l++)
a[i][j]+=f[i][l]*ans[l][j]%y;
return a[][]%y;
}
ll ph(ll x)
{
if(x<sj) return phi[x];
ll temp=x;
for(ll i=;i*i<=x;i++)
if(x%i==)
{
temp=temp-temp/i;
while(x%i==)
x/=i;
}
if(x>)
temp=temp-temp/x;
return temp;
}
ll ksm(ll x,ll y,ll z)
{
x%=z;
ll jg=;
while(y)
{
if(y&) jg=jg*x%z;
x=x*x%z;
y>>=;
}
return jg%z;
}
int main()
{
prime(sj-);
scanf("%d%lld",&m,&p);
for(int i=;i<=m;i++)
{
scanf("%lld%lld",&n,&q);
printf("%lld\n",ksm(p,fb((n-),ph(q)),q));
}
return ;
}
password
Password的更多相关文章
- 打开程序总是会提示“Enter password to unlock your login keyring” ,如何成功关掉?
p { margin-bottom: 0.1in; line-height: 120% } 一.一开始我是按照网友所说的 : rm -f ~/.gnome2/keyrings/login.keyrin ...
- your password has expired.to log in you must change it
今天应用挂了,log提示密码过期.客户端连接不上. 打开mysql,执行sql语句提示密码过期 执行set password=new password('123456'); 提示成功,但客户端仍然连接 ...
- MySql Access denied for user 'root'@'localhost' (using password:YES) 解决方案
关于昨天下午说的MySQL服务无法启动的问题,解决之后没有进入数据库,就直接关闭了电脑. 今早打开电脑,开始-运行 输入"mysql -uroot -pmyadmin"后出现以下错 ...
- [上架] iOS "app-specific password" 上架问题
当你的 Apple ID 改用双重认证密码时,上架 iOS App 需要去建立一个专用密码来登入 Apple ID 才能上架. 如果使用 Application Loader 上传时,得到这个讯息: ...
- [LeetCode] Strong Password Checker 密码强度检查器
A password is considered strong if below conditions are all met: It has at least 6 characters and at ...
- mysql 错误 ERROR 1372 (HY000): Password hash should be a 41-digit hexadecimal number 解决办法
MySQL创建用户(包括密码)时,会提示ERROR 1372 (HY000): Password hash should be a 41-digit hexadecimal number: 问题原因: ...
- phpmyadmin #1045 - Access denied for user 'root'@'localhost' (using password: NO)
phpmyadmin访问遇到1045问题 #1045 - Access denied for user 'root'@'localhost' (using password: NO) 解决办法 找到p ...
- 保留password模式文本框textbox内的数据不丢失。
在asp.net 2.0环境下,使用textbox,提交到服务器再传回,如果textbox是password模式的,那么textbox内的密码(星号),就没有了! protected override ...
- Windows mysql提示:1045 access denied for user 'root'@'localhost' using password yes
Windows mysql提示:1045 access denied for user 'root'@'localhost' using password yes http://blog.csdn.n ...
- ERROR 1819 (HY000): Your password does not satisfy the current policy requirements
为了加强安全性,MySQL5.7为root用户随机生成了一个密码,在error log中,关于error log的位置,如果安装的是RPM包,则默认是/var/log/mysqld.log. 一般可通 ...
随机推荐
- 关于MATLAB处理大数据坐标文件2017529
今天我们离成功又近了一步,因为又失败了两次 第一次使用了所有特征,理由:前天的特征使用的是取单个特征测试超过85分的特征,结果出现过拟合现象. 本次使用所有特征是为了和昨天的结果作比较. 结果稍好:比 ...
- php取出数组中的最大值
<?php /** * @param $arr * @return mixed * php取出数组中的最大值(方法一) */ function getMax($arr){ $max=$arr[0 ...
- USACO Dynamic Programming (1)
首先看一下题目: Introduction Dynamic programming is a confusing name for a programming technique that drama ...
- ==,=和equals()区别
equals和=,==的区别 一. ==和equals的区别 1. ==是运算符 2. equals是String对象的方法 一般有两种类型的比较 1. 基本数据类型的比较 2. 引用对象的比较 ...
- H3CNE实验:配置交换机接口
第1步:配置交换机端口 <H3C>system-view System View: return to User View with Ctrl+Z. [H3C]interface Giga ...
- voa 2015 / 4 / 26
Now, Words and Their Stories, a VOA Special English program about American expressions. I'm Rich Kle ...
- DELPHI XE5安装
1.安装XE5 2. HNFJ-DPADCW-BDWCFU-FPNN QDF4-CTSDHV-RDFCFE-FEAN HNFK-BCN8NN-78N53D-H4RS 破解补丁使用方法: (1).复制压 ...
- Socket异步通信及心跳包同时响应逻辑分析。
有段时间没有更博了,刚好最近在做Socket通信的项目,原理大致内容:[二维码-(加logo)]-->提供主机地址和端口号信息(直接使用[ThoughtWorks.QRCode.dll]比较简单 ...
- WinForm中AssemblyInfo.cs文件参数具体讲解
在.NET中有一个配置文件AssemblyInfo.cs主要用来设定生成的有关程序集的常规信息dll文件的一些参数,下面是默认的AssemblyInfo.cs文件的内容具体介绍 //是否符合公共语言规 ...
- 初学angular
1.angular 表达式 2.ng-app ng-init ng-model ng-repeat ng-model是用于表单元素的,支持双向绑定.对普通元素无效: ng-bind用于普通 ...