函数

最新推荐文章于 2024-07-30 22:55:11 发布

magical_qting

最新推荐文章于 2024-07-30 22:55:11 发布

阅读量507

点赞数 1

CC 4.0 BY-SA版权

分类专栏：数论文章标签：数论

本文链接：https://blog.csdn.net/magical_qting/article/details/48132719

数论专栏收录该内容

3 篇文章

订阅专栏

本文介绍了一个数学挑战问题，即根据给定的n、q和p，计算Fn%p的值，其中Fn由q和n定义。文章给出了计算公式及算法实现思路，包括等比数列求和的方法。

Description

东东上幼儿园学习了数学，刚刚学会1+1=2，就让教信息学竞赛的爸爸给他出题，并且扬言越难越好，好吧，为了让东东知道自己对数学理解的肤浅，东东爸给出下面的题目:

$F_n=\frac {q^n+1}{q+1}$ （n为奇数）

$F_n=\frac{q^n -1}{q+1}$ （n为偶数）

给定n和q，求 $F_n$ 的值。
东东不畏艰难困苦，终于没有算出来，同时也懂得了数学的博大精深。你能算出来么？

Input

第一行三个整数x、y、z来描述q，q= $x^y$ +z；
$第二行一个正整数n；$
$第三行一个正整数p。$

Output

输出 $F_n$ % p 的值。

Hint

对于60%的数据， 0<=x,y<=50000,1<=z<=50000, 0<=n<=105, 1<=p<108
对于100%的数据，0<=x,y<=50000,1<=z<=50000,0<=n<=109,1<=p<108

Sample Input

3 0 5
3
70546

Sample Output

Key To Problem

设s= $q^2$
当n为奇数时
$F_n=q^{n-1}-q^{n-2}+q^{n-3}-q^{n-4}+……+q^2-q^1+q^0$
则 $F_n=(s^0+s^1+s^2+……+s^{n/2-1})*(s-q)$
当n为偶数时
$F_n=q^{n-1}-q^{n-2}+q^{n-3}-q^{n-4}+……+q^1-q^0$
则 $F_n=(s^0+s^1+s^2+……+s^{n/2-1})*(q-1)$
之后就是简单的等比数列求和啦，可以用exgcd来搞，然而我用的是二分，不过这些都是细节就不要在意啦

CODE

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
ll a1,a2,a3;
int n,MOD;

ll pow(ll x,int y)
{
    ll ans=1;
    while(y)
    {
        if(y&1)
            ans=(x*ans)%MOD;
        x=(x*x)%MOD;
        y>>=1;
    }
    return ans;
}

ll ef(ll x,int sum)
{
    if(sum<0)
        return 0;
    if(sum==0)
        return 1;
    ll ans=0;
    if(sum%2==0)
        ans=((pow(x,sum/2+1)+1)*ef(x,sum/2-1)%MOD+pow(x,sum/2))%MOD;
    else
        ans=(pow(x,sum/2+1)+1)*ef(x,sum/2)%MOD;
    return ans%MOD;
}

int main()
{
//  freopen("function.in","r",stdin);
//  freopen("function.out","w",stdout);
    cin>>a1>>a2>>a3>>n>>MOD;
    ll q=(pow(a1,(int)a2)+a3)%MOD;
    ll s=(q*q)%MOD;
    if(n%2==1)
    {
        ll a=ef(s,n/2-1)*(s-q+MOD)%MOD+1;
        cout<<a%MOD<<endl;
    }else
    {
        ll a=ef(s,n/2-1)*(q-1);
        cout<<a%MOD<<endl;
    }
    return 0;
}