当前位置:网站首页>神奇快速幂

神奇快速幂

2022-07-07 21:53:00 sophilex

 思路;
纯模拟肯定是送,有很多人都是打表找规律,找循环节然后过去的。

在网上看到一个秀翻我的思路:把一次传染过程看作乘法操作,操作对象为标记数组,那么经历了k次操作,也就是乘了k次,这就是一个变形的快速幂

那么只要按照快速幂敲一下就好了

整体复杂度大概是m^2*log(k),4e7的水平,理论上是能过的

but,时间还是卡的很死的(可能常数有点大),必须得把所有杂七杂八的优化全部加上去:快读,register,inline。。。

优化到极致了,才能过。。。(亲身体会)

虽然但是,这个思路确实是我第一次见,很优秀

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define lowbit(x) x&(-x)
const int N=1500;
ll k;
int m,n;
inline ll read() {
    register char ch;
    while(!isdigit(ch=getchar()));
    register ll x=ch^'0';
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}
bool mas[N],tp[N],ans[N];
inline void mul(bool a[],bool b[])
{
	memset(tp,0,sizeof tp);
	for(register int i=0;i<m;++i)
	{
		for(register int j=0;j<m;++j)
		{
			tp[i*j%m]|=a[i]&&b[j];
		}
	}
	std::copy(&tp[0],&tp[m],a);
}
int main()
{//ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
	k=read();m=read();n=read();
	for(register int i=1;i<=n;++i)
	{
		//a=read();
		mas[read()]=1;
	}
	ans[1]=1;
	while(k)
	{
		if(k&1) mul(ans,mas);
		mul(mas,mas);
		k>>=1;
	}
	for(register int i=0;i<m;++i) if(ans[i]) printf("%d ",i);
	return 0; 
 }

原网站

版权声明
本文为[sophilex]所创,转载请带上原文链接,感谢
https://blog.csdn.net/sophilex/article/details/125646744