Side Transmutations(组合数学+快速幂)

原题:

题目不短,但是容易理解。
最后题目要统计的是“不相等”的字符串个数,是一道组合计数的题目,显然需要我们把各个部分的情形数乘起来得到答案。

最容易算的部分是剪掉两端各b[m]后留下的不能被翻转的“中间”部分,每个字符都有 ∣ A ∣ \left | A \right | A种可能,故共有 ∣ A ∣ n − 2 ∗ b m \left | A \right |^{n-2*b_{m}} An2bm种情形

两边可翻转的部分没办法一眼看出计算方法,我们可以从简单的m=1入手思考
设翻转的字符数为r,则一段长度为r的序列共有 S = ∣ A ∣ r S=\left | A \right |^r S=Ar种可能
我们把左边要翻转的串叫做L,右边叫做R
情形一:R翻转后=L,这种情况恰有S种
情形二:R翻转后≠L,那么翻转前后视为同一个串,这种情况有 C S 2 C_{S}^{2} CS2
故一共有 C S 2 + S = S ( S + 1 ) 2 C_{S}^{2}+S=\frac{S(S+1)}{2} CS2+S=2S(S+1)种情形

再对m进行推广,注意到 b i + 1 b_{i+1} bi+1的翻转对于 b i b_i bi内的情形数不产生影响,因此只需要计算 b i + 1 b_i+1 bi+1 b i + 1 b_{i+1} bi+1的情形数就行了(即 r = b i + 1 − b i r=b_{i+1}-b_i r=bi+1bi时的情形数)

最后审视全部计算过程,发现只需要快速幂就足够完成计算了,时间复杂度 O ( m l o g n ) O(m logn) O(mlogn)

代码:

#include <bits/stdc++.h>
typedef long long ll;
using namespace std;
const ll PP=998244353;
const int MAXN=2e5+10;
int n,m,A,i,b[MAXN];
ll ans,t,p,x,f[10000];

ll ex(ll x,int y) //calc x*(A^y) 
{
	ll i=0;
	while (y)
	{
		if (y&1) x=(x*f[i])%PP;
		y=y>>1;
		i++;
	}
	return x;
}	

int main()
{
	cin>>n>>m>>A;
	for (i=1;i<=m;i++) scanf("%d",&b[i]);
	A=A%PP;
	f[0]=A%PP;b[0]=0;
	for (i=1;i<=100;i++) f[i]=(f[i-1]*f[i-1])%PP;
	ans=1;t=1;
	for (i=1;i<=m;i++)
	{
		t=ex(1,b[i]-b[i-1]);
		t=(t*(t+1)/2)%PP;
		ans=(ans*t)%PP;
	}
	ans=ex(ans,n-2*b[m]);
	cout<<ans;
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
#include<graphics.h> #define cube_size 20 //设置单个网格的大小 #define side 80 //设置一共有多少个网格 bool is_run = true; bool is_pause = false; bool field[side + 1][side + 1]; bool ass[side + 1][side + 1]; int d_x[3] = { -1,0,1 }; int d_y[3] = { -1,0,1 }; bool update(int x, int y) { int sum = 0; for (int i = 0; i < 3; i++) for (int j = 0; j < 3; j++) { int xt = x + d_x[i], yt = y + d_y[j]; if (field[xt][yt]) sum += 1; } if (field[x][y]) { if (sum < 3 || sum>4) return false; else return true; } else { if (sum == 3 || sum == 4) return true; else return false; } } int main() { initgraph(cube_size * side, cube_size * side); setbkcolor(BLACK); cleardevice(); setfillcolor(GREEN); setlinecolor(WHITE); for (int i = 0; i < side; i++) { line(i * cube_size, 0, i * cube_size, side * cube_size); line(0, i * cube_size, side * cube_size, i * cube_size);; } while (is_run) { Sleep(500); // 画图 for (int i = 0; i < side; i++) for (int j = 0; j < side; j++) { if (field[i][j]) fillrectangle(i * cube_size, j * cube_size, (i + 1) * cube_size, (j + 1) * cube_size); else clearrectangle(i * cube_size + 1, j * cube_size + 1, (i + 1) * cube_size - 1, (j + 1) * cube_size - 1); } //鼠标操控游戏 MOUSEMSG m; MOUSEMSG* pm = &m; while (PeekMouseMsg(pm, 1)) { if (m.mkRButton) is_pause = true; } while (is_pause)//鼠标右键暂停 { m = GetMouseMsg(); if (m.mkLButton) { int x = m.x / cube_size, y = m.y / cube_size; if (field[x][y]) { field[x][y] = false; clearrectangle(x * cube_size + 1, y * cube_size + 1, (x + 1) * cube_size - 1, (y + 1) * cube_size - 1); } else { field[x][y] = true; fillrectangle(x * cube_size, y * cube_size, (x + 1) * cube_size, (y + 1) * cube_size); } } else if (m.mkRButton)//ctrl键加鼠标右键初始化网格 { if (m.mkCtrl) { for (int i = 0; i < side; i++) for (int j = 0; j < side; j++) { field[i][j] = false; fillrectangle(i * cube_size, j * cube_size, (i + 1) * cube_size, (j + 1) * cube_size); } } else is_pause = false; } } for (int i = 0; i < side; i++) for (int j = 0; j < side; j++) ass[i][j] = update(i, j); for (int i = 0; i < side; i++) for (int j = 0; j < side; j++) field[i][j] = ass[i][j]; } return 0; }请对上述代码添加以下功能,定义类Automaton对仿真建模,并且可以暂停和恢复仿真
06-11

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值