首页 > 代码库 > 明天补完
明天补完
/* cogs 943. [東方S3] 铃仙?优昙华院?稻叶 概率dp 貌似做麻烦了 邻接矩阵和链式前向星都用上了。。。 f[pos][i][j]表示 第i秒 在i点 上一个经过的点是j 方程: f[pos][i][j]=Σf[pos-1][j][k]*(__out[j]+1-(map[j][k]==1))+f[pos-1][i][j]*(__out[i]+1-(map[i][j]==1) 前面的求和考虑的是上一秒从其他的一个节点走过来 后面的考虑的是上一秒选择停留在原地 复杂度O(T * N^3) */ #include <cstdio> void read (int &now) { register char word = getchar (); for (now = 0; word < ‘0‘ || word > ‘9‘; word = getchar ()); for (; word >= ‘0‘ && word <= ‘9‘; now = now * 10 + word - ‘0‘, word = getchar ()); } #define Max 55 bool map[Max][Max]; int __out[Max * 200]; double dp[Max * 10][Max][Max]; int __next[Max * 200]; int __to[Max * 200]; int edge_list[Max * 200]; int Edge_Count ; inline void AddEdge (int from, int to) { Edge_Count ++; __next[Edge_Count] = edge_list[from]; __to[Edge_Count] = to; edge_list[from] = Edge_Count; } int main (int argc, char *argv[]) { freopen ("reisen.in", "r", stdin); freopen ("reisen.out", "w", stdout); int N, M, T; read (N); read (M); read (T); int x, y; for (register int i = 1; i <= M; i ++) { read (x); read (y); map[x][y] = true; __out[x] ++; AddEdge (x, y); } dp[0][1][1] = 1.00; register bool flag; for (register int pos = 0, i, j, Flan; pos <= T; pos ++) for (i = 1; i <= N; i ++) for (j = 1; j <= N; j ++) if (dp[pos][i][j] > 0.00) { flag = false; if (map[i][j]) flag = true; for (Flan = edge_list[i]; Flan; Flan = __next[Flan]) { if (__to[Flan] == j) continue; dp[pos + 1][__to[Flan]][i] += dp[pos][i][j] * 1.00 / (double)(__out[i] + !flag); } dp[pos + 1][i][j] += dp[pos][i][j] * 1.00 / (double)(__out[i] + !flag); } register double Answer; for (register int i = 1, j; i <= N; i ++) { Answer = 0.00; for (j = 1; j <= N; j ++) Answer += dp[T][i][j]; printf ("%.3lf\n", Answer * 100); } return 0; }
#include <cstdio> void read (int &now) { register char word = getchar (); for (now = 0; word < ‘0‘ || word > ‘9‘; word = getchar ()); for (; word >= ‘0‘ && word <= ‘9‘; now = now * 10 + word - ‘0‘, word = getchar ()); } int N, M, R, C; #define Max 1020 long long sum[Max][Max]; inline long long max (long long a, long long b) { return a > b ? a : b; } int main (int argc, char *argv[]) { freopen("aya.in","r",stdin); freopen("aya.out","w",stdout); read (N); read (M); read (R); read (C); int x; for (register int i = 1, j; i <= N; i ++) for (j = 1; j <= M; j ++) { read (x); sum[i][j] = sum[i][j - 1] + sum[i - 1][j] + x - sum[i - 1][j - 1]; } long long Answer = -1; for (register int i = R, j; i <= N; i ++) for (j = C; j <= M; j ++) Answer = max (Answer, sum[i][j] - sum[i - R][j] - sum[i][j - C]+ sum[i - R][j - C]); printf ("%d", Answer); return 0; }
明天补完
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。