数列

题目描述
数列{xn}的递推公式如下:
给定 x 0 , a, b, c, n, m,求 x n mod m 的值。
输入格式
一行,六个整数,分别表示 x 0 , a, b, c, n, m。
输出格式
一行,一个整数,表示 x n mod m 的值。
输入样例
1 1 1 1 5 1000000000
输出样例
133904603

数据规模与约定

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long lol;
lol a,b,c,n,Mod,x,now,last,f[1000001],st,len,ans;
lol qpow(lol x,lol y,lol p)
{
  lol res=1;
  while (y)
    {
      if (y&1) res=res*x%p;
      x=x*x%p;
      y>>=1;
    }
  return res;
}
int vis[1000001];
int main()
{
    int i;
    lol p,inv;
    cin>>x>>a>>b>>c>>n>>Mod;
    memset(vis,-1,sizeof(vis));
    if (n<=1000000)
    {
        x%=Mod;
        last=x;
        for (i=1; i<=n; i++)
        {
            now=(a*last%Mod*last%Mod+b*last%Mod)%Mod+c;
            now%=Mod;
            last=now;
        }
        cout<<last;
    }
    else if (Mod<=10000000)
    {
        x%=Mod;
        last=x;
        vis[x]=0;
        f[0]=x;
        for (i=1; i; i++)
        {
            now=(a*last%Mod*last%Mod+b*last%Mod)%Mod+c;
            now%=Mod;
            last=now;
            f[i]=now;
            if (vis[now]!=-1) break;
            vis[now]=i;
        }
        st=vis[f[i]];
        len=i-st;
        while (n>=i)
        {
            n-=len;
        }
        cout<<f[n];
    }
    else
      {
    p=qpow(2,n,Mod-1);
    inv=qpow(2*a,Mod-2,Mod);
    now=a*((x+b*inv%Mod)%Mod)%Mod;
    ans=qpow(now,p,Mod);
    ans=(ans*qpow(a,Mod-2,Mod)%Mod-b*inv%Mod+Mod)%Mod;
    cout<<ans;
      }
}