The 2019 Asia Nanchang First Round Online Programming Contest H. The Nth Item 矩阵快速幂+打表优化

题意：q次查询，输出所有的查询异或的结果

每次查询是的第n项,查询的n是上一次的结果的平方异或上一次的n

题解:首先要发现查询fn的循环节是mod-1，然后这样还不够因为查询次数太多了时间还是不行

这时可以把取模数取小一点（例如100，或者1000），会发现q次查询ans在过程中会来回震荡，即两个数来回变

这时候就要跳出循环了

然后别忘了，因为发现循环节是mod-1，在算快速幂的时候别忘了模mod-1

最后，数组不要开太大，开内存会占很多时间

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
#define mem(s) memset(s, 0, sizeof(s))
const int INF = 0x3f3f3f3f;
const double eps = 1e-8;
const int maxn = 1e6+5;
const int mod = 998244353;
struct Mat
{
    ll m[3][3];
};//存储结构体
Mat a,e; //a是输入的矩阵，e是输出的矩阵
Mat Mul(Mat x,Mat y)
{
    Mat c;
    memset(c.m,0,sizeof(c.m));
    for(int i=1;i<=2;++i){
        for(int j=1;j<=2;++j){
            for(int k=1;k<=2;++k){
                c.m[i][j] = (c.m[i][j] + x.m[i][k]*y.m[k][j])%mod;
            }
        }
    }
    return c;
}
inline Mat pow(Mat x,ll y)//矩阵快速幂
{   
    y%=998244352;
    Mat ans = e;
    if(y==0)return ans.m[2][1]=0,ans;
    if(y==1)return ans.m[2][1]=1,ans;
    if(y==2)return ans.m[2][1]=3,ans;
    while(y){
        if(y&1) ans = Mul(ans,x);
        x = Mul(x,x);
        y>>=1;
    }
    return ans;
}
void init(){
    e.m[1][1]=1;e.m[1][2]=0;
    e.m[2][1]=0;e.m[2][2]=1;
    a.m[1][1]=3;a.m[1][2]=2;
    a.m[2][1]=1;a.m[2][2]=0;
}
int main() 
{
    // #ifndef ONLINE_JUDGE    //if not define 如果没有定义这个的话就执行下面
	// freopen("input.in", "r", stdin);   //只改变输入流的文件指针，读入这个文件的内容（必须要有input这个文件）stdin是标准输入流的文件指针
	// freopen("output1.out", "w", stdout);  //只改变输出流的文件指针，写入output内（如果没有output这个文件就会自动生成）stdout是标准输出流的文件指针
	// #endif
    ll q,n,ans=0;
    cin>>q>>n;
    ll ln=-1,lans=-1,lln=-1;
    for(int i=1;i<=q;i++){
        init();
        Mat f=pow(a,n);
        ll tmp=f.m[2][1];
        lln=ln;ln=n;lans=ans;
        n=n^(1ll*tmp*tmp);
        ans^=tmp;
        if(n==lln){
            if((q-i+1)&1){printf("%lld\n",ans);return 0;}
            else {printf("%lld\n",lans);return 0;}
        }
    }
    cout<<ans<<endl;
    return 0;
}