今天的比赛打的好气啊,b题最短路模板题提议理解错了一直wa,赛后一看题解瞬间感觉自己好傻逼,,,还有H提矩阵快速幂最初想的简单了把地推矩阵写出来之后交了一发1998ms tle后来乱改了一下还是tle,赛后看了一下大佬的解题思路,就说找循环节,如果有就直接按循环节输出,没有就直接输出最后一个值,抑或值保存在数组里面。b题代码就不贴了,大家估计看了提议就都会了
在此直接附上H’题代码:
#include<bits/stdc++.h>
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;
}
map<ll,int>mp;
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]);
}
}