【SDOI2013 R1 Day1】随机数生成器

传送门:NKOJ2637
这是一道基础的数学题,然而本蒟蒻WA了好久(果然还是太菜了吗qaq)。
以下是题解:
这里写图片描述

注意要特判:

放上渣代码:
这里写图片描述
我不太明白改成注释里的为何会WA掉,求解答。=-=

/*By F.F.Chopin*/ 
#include<stdio.h>
#include<iostream>
#include<cstdio>
#include<map>
#include<cmath>
#define LL long long
using namespace std;
map<LL,LL>MMap;
int T;
LL p,A,B,x1,t,c,d;
LL ksm(LL a,LL b){
    LL ans=1;
    while(b){
        if(b&1)ans=ans*a%p;
        b>>=1,a=a*a%p; 
    }
    return ans;
}
void bsgs(LL x,LL y){
    LL m=ceil(sqrt(p-1)),ex=ksm(ksm(x,m),p-2),k=1;
    MMap.clear();
    MMap[1]=0;
    //MMap[1]=m+1;
    for(LL i=1;i<m;i++){
        k=k*x%p;
        if(!MMap.count(k))MMap[k]=i;
        //if(!MMap[k])MMap[k]=i;
    }
    for(LL i=0;i<m;i++){
        if(MMap.count(y)){
            printf("%lld\n",i*m+MMap[y]+1);
            return;
        }
        y=y*ex%p;
    }
    puts("-1");
    return;
}
bool Special_Judge(){
    if(x1==t)return puts("1"),1;
    if(!A){
        if(B==t)return puts("2"),1;
        else return puts("-1"),1;
    }
    if(A==1){
        if(!B)return puts("-1"),1;
        LL tmp=(B+(t-x1))%p*ksm(B,p-2)%p;
        if(tmp<=0)tmp+=p;
        printf("%lld\n",tmp);
        return 1;
    }
    return 0;
}
int main(){
    scanf("%d",&T);
    while(T--){
        scanf("%lld%lld%lld%lld%lld",&p,&A,&B,&x1,&t);
        if(Special_Judge())continue;//Naive!
        c=ksm((x1*(A-1)%p`
B+p)%p,p-2);
        d=(B+(A-1)*t+p)%p;
        bsgs(A,c*d%p);
    }
    return 0;
}
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值