组合数问题（复数快速幂， 组合数）

牛客

思路
参考题解， 众所周知组合数是满足二项式定理的, 这里有：
$(x+y{)}^n=C_n^0{x}_0{y}_n+C_n^{1}x{y}^{n-1}+...+C_n^{n-1}{x}_{n-1}y+C_n^n{x}_n{y}_0$
考虑复数有
$i^1=i$
$i^2=-1$
$i^3=-i$
$i^4=1$
因此当$x=1,y=-i$，该二项式为：
$（1+i{)}^n=C_n^0{i}^n+C_n^1{i}^{n-1}+C_n^2{i}^{n-2}+C_n^3{i}^{n-3}+...+C_n^{n-1}i+C_n^n{i}^0$
当$x=1,y=-i$,该二项式为
$(1-i{)}^n=C_n^0(-i{)}^n+C_n^1(-i{)}^{n-1}+C_n^2(-i{)}^{n-2}+C_n^3(-i{)}^{n-3}+...+C_n^{n-1}(-i)+C_n^n(-i{)}^0=（1+i{)}^n=C_n^0{i}^n-C_n^1{i}^{n-1}+C_n^2{i}^{n-2}-C_n^3{i}^{n-3}+...-C_n^{n-1}i+C_n^n{i}^0$

两者相加有
$(1+i{)}^n+(1-i{)}^n=2(C_n^0-C_n^2+C_n^4-C_n^6+...+C_n^n)$
$(1+1{)}^n+(1-1{)}^n=2(C_n^0+C_n^2+C_n^4+C_n^6+...+C_n^n)=2^n$

显然
$2^n+(1+i{)}^n+(1-i{)}^n=4\ast (C_n^0+C_n^4+C_n^8+C_n^{12}+....)$
利用复数快速幂模板可以很快求解

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;

const int mod = 998244353;

typedef long long ll;
struct Complex {
	ll a, b; // a + bi
	Complex(ll _a, ll _b) {
		a = _a, b = _b;
	}
};

Complex mul(Complex A, Complex B) {
	ll a = ((A.a * B.a - A.b * B.b) % mod + mod) % mod;
	ll b = ((A.a * B.b + A.b * B.a) % mod + mod) % mod;
	return Complex(a, b); 
} 

Complex C_qp(Complex A, ll b) {
	Complex ans(1, 0);
	while(b) {
		if(b & 1) ans = mul(ans, A);
		A = mul(A, A);
		b >>= 1;
	}
	return ans;
}

ll qp(ll a, ll b) {
	ll ans = 1;
	while(b) {
		if(b & 1) ans = ans * a % mod;
		a = a * a % mod;
		b >>= 1;
	}
	return ans;
}

ll n;
int main()
{
	scanf("%lld", &n);
	ll res = qp(2, n);
	
	Complex A(1, 1);
	Complex B(1, -1);
	A = C_qp(A, n);
	B = C_qp(B, n);
	
	res = (res + A.a + B.a) % mod;
	res = res * qp(4, mod - 2) % mod;
	res = (res + mod) % mod;
	
	printf("%lld\n", res);
	
	return 0;
}