【题解】[CJOI2019]Cipher

【题解】[CJOI2019]Cipher

题目描述

给定你\(p\)进制数\(s\),\(p \le 9+26\),求对于十进制数\(k\),求\(k^s \equiv ? \mod m\)

数据范围

\(len(s)\le10^5\)

\(Solution\)

我给的题意其实很简略了,代码里一些细枝末节有所不同,比如\(s\)读入和\(p\)的确定。

考虑题目\(k^s\)代数变形:,\(a_i\)表示\(s\)\(i\)位上的数。

\[ k^{\Sigma a_ip^{i-1}}=\prod k^{a_ip^{i-1}}=\prod ({k^{p^{i-1}}})^{a_i} \dots \ (1) \]

其中,可以有恒等式
\[ k^{p^{i+1}}=(k^{p^i})^p \]

记为
\[ f(i)=k^{p^i}, \]

则有

\[ f(i+1)=f(i)^p \]

所以

\[ (1)=\prod f(i-1)^{a_i} \]

预处理\(f(i)​\) ,快速幂处理所有的次方,复杂度\(O(nlogn)​\)

#include<bits/stdc++.h>

using namespace std;
#define RP(t,a,b) for(register int t=(a),edd=(b);t<=edd;++t)
#define DRP(t,a,b) for(register int t=(a),edd=(b);t>=edd;--t)
#define ERP(t,a) for(register int t=head[a];t;t=e[t].nx)
#define Max(a,b) ((a)<(b)?(b):(a))
#define Min(a,b) ((a)<(b)?(a):(b))
#define midd register int mid=(l+r)>>1
#define TMP template < class ccf >
typedef long long ll;
int cnt;
#define int ll
TMP inline ccf qr(ccf b){
    char c=getchar();
    int q=1;
    ccf x=0;
    while(c<48||c>57)
    q=c==45?-1:q,c=getchar();
    while(c>=48&&c<=57)
    x=x*10+c-48,c=getchar();
    return q==-1?-x:x;
}
ll k,mod;
ll p;
string t1;
const int maxn=1e5+15;
int data[maxn];


inline int trans(char x){
    if(x>=48&&x<=57)
    return x-48;
    return x-'a'+10;
}

inline ll ksm(ll base,ll x){
    ll ret=1;
    while(x){
    if(x&1ll)
        ret*=base,ret%=mod;
    base*=base;
    base%=mod;
    x>>=1;
    }
    return ret%mod;
}

inline void wk(){
    ll ans=1;
    ll base=k;
    RP(t,1,cnt){
    ans*=ksm(base,data[t]);
    base=ksm(base,p);
    ans%=mod;
    }
    cout<<ans%mod<<endl;
    return;
}

string qaq;
inline void qrqr(){
    p=0;
    cnt=0;
    memset(data,0,sizeof data);
    DRP(t,qaq.length()-1,0){
    cnt++;
    data[cnt]=trans(qaq[t]);
    p=Max(p,data[cnt]+1);
    }
}

signed main(){

    freopen("cipher.in","r",stdin);
    freopen("cipher.out","w",stdout);
    
    k=qr(1ll);
    mod=qr(1ll);
    
    while(cin>>qaq)
    qrqr(),wk();
    
    return 0;
}

转载于:https://www.cnblogs.com/winlere/p/10359868.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值