Description

B 君有兩個好朋友，他們叫寧寧和冉冉。有一天，冉冉遇到了一個有趣的題目：輸入 b;d;n，求

Input

一行三個整數 b;d;n

Output

一行一個數表示模 7528443412579576937 之後的結果。

Sample Input

Sample Output

HINT

其中 0<b^2< = d<(b+1)2< = 10^18,n< = 10^18，並且 b mod 2=1,d mod 4=1

#include<bits/stdc++.h>
#define ll unsigned long long
using namespace std;
const ll ha=7528443412579576937ll;

inline ll add(ll x,ll y){
	x+=y;
	return x>=ha?x-ha:x;
}

inline ll ksc(ll x,ll y){
	ll an=0;
	for(;y;y>>=1,x=add(x,x)) if(y&1) an=add(an,x);
	return an;
}

inline ll ksm(ll x,ll y){
	ll an=1;
	for(;y;y>>=1,x=ksc(x,x)) if(y&1) an=ksc(an,x);
	return an;
}

const ll inv=ksm(2,ha-2);
const ll INV=ksc(inv,inv);
ll B,D,N;
struct node{
	ll a[2][2];
	
	node operator *(const node &u)const{
		node r;
		for(int i=0;i<=1;i++)
		    for(int j=0;j<=1;j++){
		    	r.a[i][j]=add(ksc(a[i][0],u.a[0][j]),ksc(a[i][1],u.a[1][j]));
			}
		return r;
	}
}ans,x;

inline void solve(){
	ans.a[0][0]=ans.a[1][1]=1;
	ans.a[0][1]=ans.a[1][0]=0;
	
	ll O=N;
	N--;
	for(;N;N>>=1,x=x*x) if(N&1) ans=ans*x;
	
	ll an=0;
	an=add(ksc(2,ans.a[0][1]),ksc(B,ans.a[1][1]));
	
	if(ksc(B,B)!=D&&!(O&1)) an=add(an,ha-1);
	printf("%lld\n",an);
}

int main(){
	scanf("%lld%lld%lld",&B,&D,&N);
	if(!N){
		puts("1");
		return 0;
	}
    
	x.a[0][0]=0;
	x.a[1][0]=1;
	x.a[1][1]=B;
	x.a[0][1]=(D-ksc(B,B))>>2;
	
	solve();
	
	return 0;
}

[JLOI2015]有意義的字符串