bzoj 3122 [Sdoi2013]随机数生成器

题目链接: http://www.lydsy.com/JudgeOnline/problem.php?id=3122

感谢my的题解:http://blog.csdn.net/xaphoenix/article/details/50719539


大致题意是:给你个方程,Xn+1=(aXn+b)%P,求最小的n满足Xn=t ,P是质数。

我们将递推式展开得到: Xn+1=a^n*x1+b+a^2*b+···+a^(n-1)*b。后半部分我们用等比数列进行化简。我们设d=(b/(a-1))%P。那么Xn+1=a^n*x1+d*a^n-db。令(x1+d)为p,x1+bd为q。那么q=p*a^n。令x=q/p。则a^n=x(mod P) 到这里后我们利用BSGS进行计算即可。注意本题的细节部分:x1==t,a==0,a==1。

注意要用LL ,一开始我自己的BSGS模板过样例的时候发现少1,,,纠结了一会 那么加一再交吧!就ok。



/**************************************************************
    Problem: 3122
    User: yukun133
    Language: C++
    Result: Accepted
    Time:1832 ms
    Memory:10684 kb
****************************************************************/
 
#include<cstdio>
#include<queue>
#include<cstring>
#include<iostream>
#include<cmath>
#include<map>
using namespace std;
#define LL long long
const int N = 100100;
const int M = 1000100;
 
LL sz,num[N],pre[N],last[M];
LL quickpow(LL a,LL b,LL p){
    LL res=1;a%=p;
    while(b>0){
        if(b & 1) res=(long long)res*a%p;
        b>>=1;
        a=(long long)a*a%p;
    }
    return res;
}

void Hash(LL x,LL y){
    sz++;num[sz]=y;pre[sz]=last[x];last[x]=sz;
}
LL BSGS(LL y,LL z,LL p){ //return -1表示无解   
    memset(last,0,sizeof(last));
    sz=0;z%=p;
    if(z==1) return 0;
    if(y%p==0){
        if(z==0) return 1;
        else return -1;
    }
    LL i,j,mi,limt=sqrt(p),now=z;
    for(i=1;i<=limt+1;i++){
        Hash(now%M,now);
        now=(long long)now*y%p;
    }
    now=1;mi=quickpow(y,limt,p);
    for(i=1;i<=limt+1;i++){
        now=(long long)now*mi%p;
        for(j=last[now%M];j;j=pre[j]){
            if(num[j]==now){
                return i*limt-j+1;
            }
        }
    }
    return -1;
}

void extendgcd(LL a,LL b,LL &d,LL &x,LL &y)  
{  
    if(b==0){d=a;x=1;y=0;return;}  
    extendgcd(b,a%b,d,y,x);  
    y-=x*(a/b);  
}  
LL GetNi(LL A,LL M)  
{  
    LL rex=0,rey=0;  
    LL td=0;  
    extendgcd(A,M,td,rex,rey);  
    return (rex%M+M)%M;  
}  
 
int main()
{
    int T;  
    cin>>T;  
    LL p,a,b,x1,t;
    while (T--)  
    {  
        cin>>p>>a>>b>>x1>>t;  
        if (t==x1)  
        {  
            printf("1\n");  
            continue;  
        }  
        if (a==0)  
        {  
            if (b==t) cout<<2<<endl;  
                else  cout<<-1<<endl;  
            continue;  
        }   
        if (a==1)  
        {  
            t=(t+p-x1+b)%p;  
            t=(t*GetNi(b,p))%p;  
            if (t>0) printf("%d\n",t);  
            else if(b==0) printf("-1\n");  
            else printf("%d\n",p);  
            continue;  
        }  
        LL d=(b*GetNi(a-1,p))%p;  
        t=(t+d)%p;  
        x1=(x1+d)%p;  
        t=(t*GetNi(x1,p))%p; 
        LL ans=BSGS(a,t,p);
        if(~ans) 
            cout<<ans+1<<endl;
        else cout<<-1<<endl;
        
    }  
     
    return 0;
}
 
 
/*
 
3 
7  1 1 3 3
7  2 2 2 0
7  2 2 2 1
 
*/


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
SQLAlchemy 是一个 SQL 工具包和对象关系映射(ORM)库,用于 Python 编程语言。它提供了一个高级的 SQL 工具和对象关系映射工具,允许开发者以 Python 类和对象的形式操作数据库,而无需编写大量的 SQL 语句。SQLAlchemy 建立在 DBAPI 之上,支持多种数据库后端,如 SQLite, MySQL, PostgreSQL 等。 SQLAlchemy 的核心功能: 对象关系映射(ORM): SQLAlchemy 允许开发者使用 Python 类来表示数据库表,使用类的实例表示表中的行。 开发者可以定义类之间的关系(如一对多、多对多),SQLAlchemy 会自动处理这些关系在数据库中的映射。 通过 ORM,开发者可以像操作 Python 对象一样操作数据库,这大大简化了数据库操作的复杂性。 表达式语言: SQLAlchemy 提供了一个丰富的 SQL 表达式语言,允许开发者以 Python 表达式的方式编写复杂的 SQL 查询。 表达式语言提供了对 SQL 语句的灵活控制,同时保持了代码的可读性和可维护性。 数据库引擎和连接池: SQLAlchemy 支持多种数据库后端,并且为每种后端提供了对应的数据库引擎。 它还提供了连接池管理功能,以优化数据库连接的创建、使用和释放。 会话管理: SQLAlchemy 使用会话(Session)来管理对象的持久化状态。 会话提供了一个工作单元(unit of work)和身份映射(identity map)的概念,使得对象的状态管理和查询更加高效。 事件系统: SQLAlchemy 提供了一个事件系统,允许开发者在 ORM 的各个生命周期阶段插入自定义的钩子函数。 这使得开发者可以在对象加载、修改、删除等操作时执行额外的逻辑。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值