南昌icpc2019网络赛

今天的比赛打的好气啊，b题最短路模板题提议理解错了一直wa，赛后一看题解瞬间感觉自己好傻逼，，，还有H提矩阵快速幂最初想的简单了把地推矩阵写出来之后交了一发1998ms tle后来乱改了一下还是tle，赛后看了一下大佬的解题思路，就说找循环节，如果有就直接按循环节输出，没有就直接输出最后一个值，抑或值保存在数组里面。b题代码就不贴了，大家估计看了提议就都会了
在此直接附上H’题代码：

#include
using namespace std;
const int N = 1e6 + 50;
typedef long long ll;
struct Mat{
    ll m[2][2];
};
const int mod = 998244353;
Mat Mul(Mat a,Mat b)
{
    Mat c;
    for(int i = 0;i < 2;i++)
    {
        for(int j = 0;j < 2;j++)
        {
            c.m[i][j] = 0;
        }
    }
    for(int i = 0;i < 2;i++)
    {
        for(int j = 0;j < 2;j++)
        {
            for(int k = 0;k < 2;k++)
            {
                c.m[i][j] = (c.m[i][j]%mod + a.m[i][k]*b.m[k][j])%mod;
            }
        }
    }
    return c;
}
Mat pow(Mat x,ll y)
{
    Mat ans;
    ans.m[0][0] = ans.m[1][1] = 1;
    ans.m[0][1] = ans.m[1][0] = 0;
    while(y)
    {
        if(y&1)
        {
            ans = Mul(ans,x);
        }
        x = Mul(x,x);
        y/=2;
    }
    return ans;
}
mapmp;
ll ans[N];
int main()
{
    Mat a;
    a.m[0][0] = 3;a.m[0][1] = 2;
    a.m[1][0] = 1;a.m[1][1] = 0;
    int q;ll n;
    int now = 0;
    int yu = 0;
    int len = 0;
    scanf("%d%lld",&q,&n);
    for(int i = 1;i <= q;i++)
    {
        Mat tmp;
        ll c;
        tmp = pow(a,n-1);
        c = tmp.m[0][0];
        ans[i] = ans[i-1]^c;
        n = n^(c*c);
        if(mp[n])
        {
            now = mp[n];
            yu = q-i;//还有多少项对len取余加上now就是最后的结果
            len = i - now;//循环节长度
            break;
        }
        mp[n] = i;
    }
    if(len)
    {
        yu = yu%len;
        printf("%lld\n",ans[yu+now]);
    }
    else
    {
        printf("%lld\n",ans[q]);
    }
}