2019南昌网络赛H.The Nth Item

最新推荐文章于 2021-03-26 19:58:44 发布

*ACoder*

最新推荐文章于 2021-03-26 19:58:44 发布

阅读量141

点赞数

分类专栏：特殊技巧

本文链接：https://blog.csdn.net/FSAHFGSADHSAKNDAS/article/details/100746929

版权

特殊技巧专栏收录该内容

17 篇文章 0 订阅

订阅专栏

链接

点击跳转

题解

特征根法解除斐波那契额通项：

根据递推式 $F (n) = 3 F (n - 1) + 2 F (n - 2)$
设 $x^2 = 3x + 2$

解得 $x_1=\frac{3-\sqrt {17} }{2},x_2=\frac{3+\sqrt {17} }{2}$

设 $F(n) = c_1x_1^n + c_2x_2^n$

将 $F (0) = 0, F (1) = 1$ 带入，解得 $c1=-\frac{\sqrt {17} }{17},c2=\frac{\sqrt {17} }{17}$

所以 $-\frac{\sqrt {17} }{17} (\frac{3-\sqrt {17} }{2})^n + \frac{\sqrt {17} }{17} (\frac{3+\sqrt {17} }{2})^n$

首先可以对 $n$ 欧拉降幂，然后用 $b i g - s m a l l$ 的方法预处理 $x_1,x_2$ 的幂，每次就可以 $O (1)$ 回答了

代码

#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
#define iinf 0x3f3f3f3f
#define linf (1ll<<60)
#define eps 1e-8
#define maxn 100010
#define cl(x) memset(x,0,sizeof(x))
#define rep(_,__) for(_=1;_<=(__);_++)
#define em(x) emplace(x)
#define emb(x) emplace_back(x)
#define emf(x) emplace_front(x)
#define fi first
#define se second
#define mod 998244353ll
using namespace std;
using namespace __gnu_pbds;
typedef long long ll;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
ll read(ll x=0)
{
    ll c, f(1);
    for(c=getchar();!isdigit(c);c=getchar())if(c=='-')f=-f;
    for(;isdigit(c);c=getchar())x=x*10+c-0x30;
    return f*x;
}
ll fastpow(ll a, ll b, ll c)
{
    ll t(a%c), ans(1ll);
    for(;b;b>>=1,t=t*t%c)if(b&1)ans=ans*t%c;
    return ans;
}
ll Q, N, r=473844410, a, b, big[2][maxn], small[2][maxn], _r=fastpow(r,mod-2,mod);
int main()
{
    ll i, S=100000, lastans=0, ans=0;
    a=(3-r+mod)*fastpow(2,mod-2,mod)%mod;
    b=(3+r)*fastpow(2,mod-2,mod)%mod;
    big[0][0]=big[1][0]=1;
    small[0][0]=small[1][0]=1;
    rep(i,S)
    {
        big[0][i]=fastpow(a,S*i,mod)%mod;
        big[1][i]=fastpow(b,S*i,mod)%mod;
        small[0][i]=small[0][i-1]*a%mod;
        small[1][i]=small[1][i-1]*b%mod;
    }
    Q=read(), N=read();
    while(Q--)
    {
        N ^= lastans*lastans;
        ll n = N % (mod-1);
        lastans = _r *( -big[0][n/S]*small[0][n%S]%mod + big[1][n/S]*small[1][n%S]%mod + mod ) %mod;
        ans ^= lastans;
    }
    printf("%lld\n",ans);
    return 0;
}