HDU 4576 Robot 解题报告（循环矩阵）

    题目大意：m条指令，每次走w步。可是可以顺时针走，也可以逆时针走，问最终机器人在l，r范围内的概率。

    解题报告：n挺大，不过矩阵仍然是循环矩阵，可以在n^2复杂度内求出。

    对于走w步时矩阵的构建要注意，w可能大于n。而且顺时针和逆时针可能走到相同的位置。所以是在可能的位置+0.5。

    代码如下：

#include<cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

int mp[110];
double a[210];
double b[210];
double c[210];
int n,m;

void mul(double*x,double*y)
{
	memset(c,0,sizeof(c));
	for(int i=0;i<n;i++)
		for(int k=0;k<n;k++)
			c[i]+=x[k]*y[k>=i?(k-i):(n+k-i)];
	for(int i=0;i<n;i++)
		x[i]=c[i];
}

int main()
{

	int l,r;
	while(~scanf("%d%d%d%d",&n,&m,&l,&r) && (n|m|l|r))
	{
		memset(mp,0,sizeof(mp));
		memset(a,0,sizeof(a));
		a[0]=1;

		for(int i=0;i<m;i++)
		{
			int t;
			scanf("%d",&t);
			mp[t]++;
		}

		for(int i=1;i<=100;i++) if(mp[i])
		{
			int k=mp[i];
			memset(b,0,sizeof(b));
			b[i%n]+=0.5;
			b[(n-i%n)%n]+=0.5;
			while(k)
			{
				if(k&1)
					mul(a,b);
				mul(b,b);
				k>>=1;
			}
		}

		double ans=0;
		for(int i=l-1;i<r;i++)
			ans+=a[i];
		printf("%.4lf\n",ans);
	}
}