Project Euler #356

传送门ber~

Duan2baka在外地找不到脚头的地方

题中$g(x)=x^3-2^n{x}^2+n$可以看作是降幂的形式，即$x^3=2^n{x}^2-n$
$x^k=x^{k-3}\cdot x^3=x^{k-3}\cdot (2^n{x}^2-n)=2^n{x}^{k-1}-n{x}^{k-3}$
令$A(k)=x^k$，那么有$A(k)=2^{n}A(k+1)-nf(x-3)$，可以矩阵乘法求得$A(0),A(1),A(2)$的系数
可是很遗憾，即使知道系数，我们也不能很轻松的求得$A(1)$和$A(2)$的值
我们换个思路，设$g(x)$的三个根分别为$a,b,c$，令$f(x)=a^x+b^x+c^x$，递推式和上面的很明显是一样的

$f(0)=3$显然，接下来我们要算$f(1)$和$f(2)$的值
$(x-a)(x-b)(x-c)=x^3-2^n{x}^2+n$，待定系数，$a+b+c=2^n$,$ab+bc+ca=0$
所以$f(1)=a+b+c=2^n$,$f(2)=f(1{)}^2-2(ab+bc+ca)=4^n$
我们可以很轻松的求出$f(987654321)$
那怎么通过$a^{987654321}+b^{987654321}+c^{987654321}$求出最大根$c$的$c^{987654321}$？

我们观察一下根的分布
几何画板发现三个根分别在$(-1,0)(0,1)(2^n-1,2^n)$
让我们瞎jb证一下。。。。
$g(0)=n>0,g(1)=1-2^n+n<0,g(-1)=-1-2^n+n<0$，所以$(-1,0)(0,1)$存在两根;
化简方程$g(x)=x^2(x-2^n)+n$，$x$取$2^n$时$g(x)=n>0$，$x$取$2^n-1$时$g(x)=-(2^n-1{)}^2+n<0$，所以必存在一根
设$a<b<c$，$a^{987654321}+b^{987654321}$必$<1$

$a^3-2^n{a}^2+n=0$，所以$-a^3-2^n{a}^2+n>0$，所以$|b|>|a|$，$0<a^{987654321}+b^{987654321}<1$

所以所求的$\lfloor c^{987654321}\rfloor =f(987654321)-1$

代码如下：

#include<bits/stdc++.h>
#define int long long
#define N 1000020
using namespace std;
inline int read(){
    int x=0,f=1;char c;
    do c=getchar(),f=f=c=='-'?-1:f; while(!isdigit(c));
    do x=(x<<3)+(x<<1)+c-'0',c=getchar(); while(isdigit(c));
    return x*f;
}
const int MOD=1e8;
const int paower=987654321;
int pw[N];
int ans;
struct Matrix{
    int a[4][4];
    Matrix(){memset(a,0,sizeof a);}
    inline void init(int n){
        memset(a,0,sizeof a);
        a[2][1]=1;a[3][2]=1;
        a[1][3]=MOD-n;a[3][3]=pw[n];
    }
    inline Matrix operator * (const Matrix &b){
        Matrix c;
        memset(c.a,0,sizeof c.a);
        for(int t=1;t<=3;t++)
            for(int i=1;i<=3;i++)
                for(int j=1;j<=3;j++)
                    (c.a[i][j]+=1ll*a[i][t]*b.a[t][j]%MOD)%=MOD;
        return c;
    } 
}a,tmp;
inline Matrix qpow(Matrix a,int k){
    Matrix sum=a;k--;
    while(k){
        if(k&1) sum=sum*a;
        a=a*a;
        k>>=1;
    }
    return sum;
}
 main(){
    for(int i=pw[0]=1;i<=30;i++) pw[i]=pw[i-1]<<1;
    for(int i=1;i<=30;i++){
        a.init(i);
        tmp=qpow(a,paower-2);
        ans=(ans+3ll*tmp.a[1][3]%MOD+1ll*pw[i]*tmp.a[2][3]%MOD+1ll*pw[i]*pw[i]%MOD*tmp.a[3][3]%MOD+MOD-1)%MOD;
    }
    printf("%lld\n",ans);
    return 0;
}