题目大意: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);
}
}