1009-Paperfolding = HDU6822
题目大意
根据从左到右的长度为n的字符串执行一条操作指令,该字符串仅包含L,R,U,D的4个不同字符。
- L指令表示将其从左向右折叠,
- R指令表示从右到左折叠,
- U指令是指从上到下折叠,
- D指令表示从下到上对折。
例如,如果从左到右折叠纸张,则应让左侧重叠在右侧,并且不能旋转。
现在,在折叠后的纸张可见部分的中心水平切割(完全切割),然后垂直切割(完全切割)。切割后所得纸张的张数为num(S)。
请参阅示例和图片以更好地理解。
现在给定一个非负整数n,可以生成
4
n
4^n
4n个不同的字符串S(字符串仅包含L,R,U,D的4个不同字符)。
求出被拆分的纸张张数的期望值,即
E
(
n
u
m
(
S
)
)
m
o
d
998244353
E(num(S))\ mod\ 998244353
E(num(S)) mod 998244353。
可以证明答案可以用
P
/
Q
P/Q
P/Q表示,其中
P
P
P和
Q
Q
Q是互质整数,并打印
P
∗
Q
−
1
m
o
d
998244353
P*Q^{-1} mod\ 998244353
P∗Q−1mod 998244353的值。
Input
第一行包含一个整数 T ( 1 ≤ T ≤ 1 0 5 ) T(1≤T≤10^5) T(1≤T≤105),即测试用例的数量。
接下来的T行中的每行都包含数字 n ( 0 ≤ n ≤ 1 0 18 ) n(0≤n≤10^{18}) n(0≤n≤1018)。
Output
对于每个测试用例,将答案打印在一行中。
Sample Input
2
0
1
Sample Output
4
6
解题思路
- 模拟一下可以看出水平对折和垂直对折的答案相对独立
- 对于 x x x次水平对折和 y y y次垂直对折
- 答案是 ( 2 x + 1 ) ( 2 y + 1 ) (2^x+1)(2^y+1) (2x+1)(2y+1)
- 因为通过反向逐操作还原,可以看到刀的痕迹的数量变化是每次在某一维倍增的
- 因此,相当于一张纸,水平和竖直分别切了 2 x , 2 y 2^x,2^y 2x,2y刀
- 所以数学期望为 E ( x ) = 1 2 n ∑ i = 0 n C n i ( 2 i + 1 ) ( 2 n − i + 1 ) = 1 + 2 n + 2 × 3 n / 2 n E(x)=\frac{1}{2^n}\sum_{i=0}^nC_n^i(2^i+1)(2^{n-i}+1)=1+2^n+2×3^n/2^n E(x)=2n1∑i=0nCni(2i+1)(2n−i+1)=1+2n+2×3n/2n
AC代码
#include<iostream>
using namespace std;
#define mod 998244353
long long quickpow( long long a, long long b, long long c) {
if (b < 0) return 0;
long long ret = 1;
a %= c;
while (b) {
if (b & 1) ret = (ret * a) % c;
b >>= 1;
a = (a * a) % c;
}
return ret;
}
long long inv( long long a, long long b) {
return quickpow(a, b - 2, b);
}
int main(){
int T;
long long n,ans=0;
cin>>T;
while(T--){
cin>>n;
long long _2n=quickpow(2,n,mod);
long long _3n=quickpow(3,n,mod);
long long x=inv(_2n,mod);
ans=(1%mod+_2n%mod+((2*_3n)%mod*x%mod)%mod)%mod;
cout<<ans<<endl;
}
return 0;
}