题意:
有个机器人,然后给你一个1~n的环:
机器人初始的时候在1,现在给你m(1e6)个操作,每次操作给出一个步数x,机器人可以等概率的选择顺时针或者逆时针走x。
现在问你,m次操作以后,机器人停落在 [ l, r ] 区间的概率是多少。
解析:
状态转移方程,come on滚动数组滚起来,( ̄Q ̄)╯。
dp[cur ^ 1][runClock(i, x)] += dp[cur][i] * 0.5; dp[cur ^ 1][runAntiClock(i, x)] += dp[cur][i] * 0.5;
开始想歪了,又是一看别人的状态转移就懂了- -。
代码:
#pragma comment(linker, "/STACK:1677721600")
#include <map>
#include <set>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <climits>
#include <cassert>
#include <iostream>
#include <algorithm>
#define pb push_back
#define mp make_pair
#define LL long long
#define lson lo,mi,rt<<1
#define rson mi+1,hi,rt<<1|1
#define Min(a,b) ((a)<(b)?(a):(b))
#define Max(a,b) ((a)>(b)?(a):(b))
#define mem(a,b) memset(a,b,sizeof(a))
#define FIN freopen("in.txt", "r", stdin)
#define FOUT freopen("out.txt", "w", stdout)
#define rep(i,a,b) for(int i=(a); i<=(b); i++)
#define dec(i,a,b) for(int i=(a); i>=(b); i--)
using namespace std;
const int mod = 1e9 + 7;
const double eps = 1e-8;
const double ee = exp(1.0);
const int inf = 0x3f3f3f3f;
const int maxn = 200 + 10;
const double pi = acos(-1.0);
const LL iinf = 0x3f3f3f3f3f3f3f3f;
int readT()
{
char c;
int ret = 0,flg = 0;
while(c = getchar(), (c < '0' || c > '9') && c != '-');
if(c == '-') flg = 1;
else ret = c ^ 48;
while( c = getchar(), c >= '0' && c <= '9') ret = ret * 10 + (c ^ 48);
return flg ? - ret : ret;
}
int n, m, l, r;
double dp[2][maxn];
int runClock(int res, int x)
{
return (res + x) % n;
}
int runAntiClock(int res, int x)
{
return ((res - x) % n + n) % n;
}
int main()
{
#ifdef LOCAL
FIN;
#endif // LOCAL
while (~scanf("%d%d%d%d", &n, &m, &l, &r))
{
if (!n && !m && !l && !r)
break;
mem(dp, 0);
dp[0][0] = 1.0;
int cur = 0;
while (m--)
{
int x = readT();
rep(i, 0, n - 1)
{
dp[cur ^ 1][i] = 0;
}
rep(i, 0, n - 1)
{
if (dp[cur][i] == 0)
continue;
dp[cur ^ 1][runClock(i, x)] += dp[cur][i] * 0.5;
dp[cur ^ 1][runAntiClock(i, x)] += dp[cur][i] * 0.5;
}
cur ^= 1;
}
if (l == 1 && r == n)
{
printf("1.0000\n");
continue;
}
double ans = 0;
rep(i, l - 1, r - 1)
{
ans += dp[cur][i];
}
printf("%.4lf\n", ans);
}
return 0;
}