sgu197:Nice Patterns Strike Back(状压dp+矩阵乘法优化)

题意:

有一个N*M的棋盘,棋盘每个格子要么是黑要么是白,每个2*2的方块颜色不能都相同,求总的颜色种数mod P。

1<=N<=10^100, 1<=M<=5, 1<=P<=10000。

分析:

1<=M<=5很明显的状压dp,但是1<=N<=10^100是我们不能一行行dp,还得用矩阵乘法优化。

对于矩阵G,G[i][j]==1表示j状态可以转移到i状态,我们dfs出矩阵G,然后G^(n-1),接着就是水了。

#include <cstdio>
#define get(s, i) ((s>>(i-1))&1)
using namespace std;
const int MAXL = 109, MAXS = 1<<5;
struct bi
{
	short a[MAXL];
	int l;	
	void read()
	{
		char c = getchar();
		if(c >= '0' && c <= '9') read();
		else return ;
		a[++l] = c-'0';	
	}
	void minus(int k)
	{
		if(a[1] >= k) a[1] -= k;
		else
		{
			a[1] += 10-k;
			a[2]--;
			int i = 2;
			while(a[i] < 0)
			{
				a[i] += 10;
				a[++i]--;	
			}
			for(int i = l; i >= 1; --i)
				if(a[i] != 0) {l = i;break;}
		}
	}
	bool empty()
	{
		if(l == 1 && a[1] == 0) return true;
		else return false;	
	}
	void binary()
	{
		if(a[1]&1) a[1]--;
		for(int i = l; i >= 1; --i)	
		{
			if(a[i]&1) a[i-1] += 10;
			a[i] >>= 1;
		}
		for(int i = l; i >= 1; --i)
			if(a[i] != 0) {l = i;break;}
	}
}n;
int m, mod;
struct mtx
{
	short g[MAXS][MAXS];
	int n, m;
}ini, conse, f;
int ans;

mtx operator * (const mtx a, const mtx b)
{
	mtx tmp = {{0}, 0, 0};
	tmp.n = a.n, tmp.m = b.m;
	for(int i = 0; i < a.n; ++i)
		for(int k = 0; k < a.m; ++k)
			for(int j = 0; j < b.m; ++j)
				tmp.g[i][j] = (tmp.g[i][j]+a.g[i][k]*b.g[k][j]%mod)%mod;	
	return tmp;
}

void make(int ss, int s, int k)
{
	if(k == m+1) 
	{
		ini.g[ss][s] = 1;	
		return ;
	}
	if(k == 1) 
	{
		make(ss, s, k+1);
		make(ss+1, s, k+1);
	}
	else 
	{
		int a = get(ss, k-1), b = get(s, k-1), c = get(s, k);
		if(a == b && b == c)
			make(ss+((a^1)<<(k-1)), s, k+1);
		else
		{
			make(ss+(1<<(k-1)), s, k+1);	
			make(ss, s, k+1);
		}
	}
}

int main()
{
	n.read();scanf("%d%d", &m, &mod);
	ini.n = ini.m = 1<<m;f.n = 1<<m, f.m = 1;
	for(int i = 0; i < (1<<m); ++i) 
	{
		make(0, i, 1);
		f.g[i][0] = 1;	
	}
	n.minus(1);
	
	conse.n = conse.m = 1<<m;
	for(int i = 0; i < conse.n; ++i)
		conse.g[i][i] = 1;
	while(!n.empty())
	{
		if(n.a[1]&1) conse = conse*ini;
		ini = ini*ini;
		n.binary(); 	
	}
	f = conse*f;
	for(int i = 0; i < f.n; ++i)
		ans = (ans+f.g[i][0])%mod;
	printf("%d\n", ans);
	return 0;
}


内容概要:本文详细介绍了如何利用Simulink进行自动代码生成,在STM32平台上实现带57次谐波抑制功能的霍尔场定向控制(FOC)。首先,文章讲解了所需的软件环境准备,包括MATLAB/Simulink及其硬件支持包的安装。接着,阐述了构建永磁同步电机(PMSM)霍尔FOC控制模型的具体步骤,涵盖电机模型、坐标变换模块(如Clark和Park变换)、PI调节器、SVPWM模块以及用于抑制特定谐波的陷波器的设计。随后,描述了硬件目标配置、代码生成过程中的注意事项,以及生成后的C代码结构。此外,还讨论了霍尔传感器的位置估算、谐波补偿器的实现细节、ADC配置技巧、PWM死区时间和换相逻辑的优化。最后,分享了一些实用的工程集成经验,并推荐了几篇有助于深入了解相关技术和优化控制效果的研究论文。 适合人群:从事电机控制系统开发的技术人员,尤其是那些希望掌握基于Simulink的自动代码生成技术,以提高开发效率和控制精度的专业人士。 使用场景及目标:适用于需要精确控制永磁同步电机的应用场合,特别是在面对高次谐波干扰导致的电流波形失真问题时。通过采用文中提供的解决方案,可以显著改善系统的稳定性和性能,降低噪声水平,提升用户体验。 其他说明:文中不仅提供了详细的理论解释和技术指导,还包括了许多实践经验教训,如霍尔传感器处理、谐波抑制策略的选择、代码生成配置等方面的实际案例。这对于初学者来说是非常宝贵的参考资料。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值