A^X mod P

最新推荐文章于 2022-05-07 22:12:42 发布

i是一只程序猿

最新推荐文章于 2022-05-07 22:12:42 发布

阅读量924

点赞数

分类专栏： ACM

本文链接：https://blog.csdn.net/zh9406/article/details/44925783

版权

ACM 专栏收录该内容

92 篇文章 0 订阅

订阅专栏

这是山东省赛的一道题，题意就和题目一样，一开始看觉得，卧槽，这尼玛数学，记忆化一下，用以前的结果往下计算，可是我错了。。。。

主要是当时卡在弱A上了，哭死。。。。。。。简直对不住队友。。。。。。。还是最后hardbird的几何模板找的及时，我敲上在最后6分钟A掉了。。。。

A^f = A^（ax+b）因为a和b是一个数量级的

这题的思路是，预处理出A的b的次方，以及(A^b)^x，然后求解，时间复杂度就减少到n了

#include<cstdio>
typedef long long LL;
using namespace std;
const int maxn =  100000;
LL n,A,m,a,b,P,K,ans,f,T;
LL dp1[100010],dp2[100010];
void solve()
{
    dp1[0]=dp2[0]=1;
    dp1[1]=A;
    for(int i=2; i<=maxn; i++)
    {
        dp1[i]=(dp1[i-1]*A)%P;
    }
    dp2[1]=dp1[100000];
    for(int i=1; i<=maxn/10; i++)
    {
        dp2[i]=(dp2[i-1]*dp2[1])%P;
    }
}

int main()
{
    scanf("%lld",&T);
    for(LL cas=1; cas<=T; cas++)
    {
        ans=0;
        scanf("%lld%lld%lld%lld%lld%lld%lld",&n,&A,&K,&a,&b,&m,&P);
        solve();
        f=K;
        for(LL i=1; i<=n; i++)
        {
            ans=(ans+dp1[f%maxn]*dp2[f/maxn])%P;
            f=(a*f+b)%m;
        }
        printf("Case #%lld: %lld\n",cas,ans);//注意下格式啊，忘了个#号，wa了n次
    }
    return 0;
}