uva 1386

最新推荐文章于 2020-09-09 22:34:38 发布

chenhongwei6447

最新推荐文章于 2020-09-09 22:34:38 发布

阅读量348

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/chenhongwei6447/article/details/51291557

题目很恶心，读不懂，百度的题意。

题目大意：一个细胞自动机，有n个格子，每个格子的取值为0~m-1。给定距离d，

每次操作后每个格子的值将变为到它距离不超过d的所有格子在操作之前的值之和除以m的余数。

普通的矩阵乘法会超时，然后就看到一个黑科技，循环矩阵

循环矩阵乘法那里细节方面感觉还是有点不大清楚

//  Created by Chenhongwei on 2016-05-01 Sunday 19:03
//  Copyright (c) 2016 Chenhongwei. All rights reserved.

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <climits>
#include <queue>
#include <cmath>
#include <map>
#include <set>
#include <stack>
#include <vector>
#include <sstream>
#include <algorithm>
#define root 1,n,1
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
using namespace std;
const int inf=1e9;
const int mod=1e9+7;
const int maxn=1e5+100;
typedef long long ll;
typedef unsigned long long ull;
ll n,m,d,k;
ll a[550];
ll tmp[550],s[550];
void mul(ll a[],ll b[])
{
	ll c[550];
	for(int i=0;i<n;i++)
	{
		c[i]=0;
		for(int j=0;j<n;j++)
			c[i]=(c[i]+a[j]*b[(i-j+n)%n])%m;
	}
	memcpy(a,c,sizeof c);
}
void fast_pow()
{
	memset(s,0,sizeof s);
	memset(tmp,0,sizeof tmp);
	s[0]=1;
	for(int i=0;i<=d;i++)
		tmp[i]=tmp[n-i]=1;
	while(k)
	{
		if(k&1)
			mul(s,tmp);
		mul(tmp,tmp);
		k>>=1;
	}
}
int main()
{
	//ios::sync_with_stdio(false);
    // freopen("in.txt","r",stdin);
	//freopen("out.txt","w",stdout);
	while(scanf("%lld%lld%lld%lld",&n,&m,&d,&k)!=EOF)
	{
		ll ans[550];
		for(int i=0;i<n;i++)
			scanf("%lld",&a[i]);
		fast_pow();
		for(int i=0;i<n;i++)
		{
			ans[i]=0;
			for(int j=0;j<n;j++)
				ans[i]=(ans[i]+a[j]*s[(i-j+n)%n])%m;
		}
		for(int i=0;i<n;i++)
			printf("%lld%c",ans[i],i==n-1?'\n':' ');
	}
	return 0;
}