Codeforces Round #131 (Div. 1) (又是各种dp)
2013-11-11 23:56
120 查看
http://codeforces.com/contest/213
A. Game
水题,拓扑排序贪心恶搞
B.
Numbers
给你一个n,要你用一些数字组成一些数,数的长度不大于n,给你0~9的cnt(i)表示组成的数种的i至少要cnt(i)个,组成的数没有前导0。
思路:首先考虑组成的数中没有0的情况,那么相当于是问1~9能组成多少个长度不大于n的数,只不过加了cnt的限制条件。这种做法是很经典的做法,相当于有x种颜色的球,每个球给出数量,问有多少种不同的排列方式。
设组成的数中i的个数为sum(i)个,总数为tot,那么很显然不同排列方式有tot!
/ (sigma(sum(i)!))。
对于这题有个不同的地方就是不能有前导0,如果照普通方法计算的话很多情况是不合法的,刚开始我是枚举总共放多少0的,不过这样子复杂度过不去,然后想了个有点巧妙的方法,先算出所有情况包括有前导0的,然后假设我放了一个0在最前面了,也就是说计算的长度都减1!比如我本来要算长度为n的,现在就是算长度为n-1的了,具体有些小细节的,这题坑的地方是居然没有把单独的0当做合法情况!
C. Relay Race
很不错的dp!
给你一个矩阵,问你一只乌龟从左上角走两遍到右下角,(x,
y)只可以走到(x+1,y),(x,y+1),中间经过的所有点的值的总和最大为多少,重复走的只计算一次。
思路:现在考虑成两只不同乌龟一起同步走的状态,设dp[i][j][k]表示左边的乌龟在第i行第j列并且右边的乌龟在j+k列,由于是同步走的,所以可以计算出右边乌龟具体的位置,然后状态转移只有四种,很容易写了接下来~具体见代码~
A. Game
水题,拓扑排序贪心恶搞
B.
Numbers
给你一个n,要你用一些数字组成一些数,数的长度不大于n,给你0~9的cnt(i)表示组成的数种的i至少要cnt(i)个,组成的数没有前导0。
思路:首先考虑组成的数中没有0的情况,那么相当于是问1~9能组成多少个长度不大于n的数,只不过加了cnt的限制条件。这种做法是很经典的做法,相当于有x种颜色的球,每个球给出数量,问有多少种不同的排列方式。
设组成的数中i的个数为sum(i)个,总数为tot,那么很显然不同排列方式有tot!
/ (sigma(sum(i)!))。
对于这题有个不同的地方就是不能有前导0,如果照普通方法计算的话很多情况是不合法的,刚开始我是枚举总共放多少0的,不过这样子复杂度过不去,然后想了个有点巧妙的方法,先算出所有情况包括有前导0的,然后假设我放了一个0在最前面了,也就是说计算的长度都减1!比如我本来要算长度为n的,现在就是算长度为n-1的了,具体有些小细节的,这题坑的地方是居然没有把单独的0当做合法情况!
// Author : JayYe Created Time: 2013-11-9 20:29:54 #include <stdio.h> #include <string.h> #include <algorithm> using namespace std; typedef __int64 ll; const int maxn = 100 + 10; const int mod = 1000000007; ll fac[maxn], dp[12][maxn], ff[maxn]; int cnt[maxn]; ll exgcd(ll a, ll b, ll &x, ll &y) { if(!b) { x = 1; y = 0; return a; } ll ret = exgcd(b, a%b, y, x); y -= a/b*x; return ret; } ll inv(ll a, ll mod) { ll x, y, d = exgcd(a, mod, x, y); if(x < 0) x += mod; return x; } void init(int n) { fac[0] = ff[0] = 1; for(int i = 1;i <= n; i++) fac[i] = fac[i-1]*i%mod; for(int i = 1;i <= n; i++) ff[i] = inv(fac[i], mod); } ll cal(int n, int tot) { for(int i = 0;i < 10; i++) for(int j = 0;j <= n; j++) dp[i][j] = 0; for(int i = cnt[0];i <= n; i++) dp[0][i] = ff[i]; for(int i = 1;i < 10; i++) { for(int j = cnt[i];j <= n; j++) { for(int k = 0;k <= n-j; k++) if(dp[i-1][k]) { dp[i][j+k] = (dp[i][j+k] + dp[i-1][k]*ff[j]%mod)%mod; } } } ll ret = 0; for(int i = 1;i <= n; i++) { ret = ret + dp[9][i]*fac[i]%mod; ret %= mod; } return ret%mod; } int main() { init(100); int n; scanf("%d", &n); int tot = 0; for(int i = 0;i < 10; i++) { scanf("%d", &cnt[i]); tot += cnt[i]; } if(tot > n) return puts("0"), 0; ll ans = cal(n, tot); if(!tot) ans--; if(cnt[0]) { cnt[0]--; ans -= cal(n-1, tot-1); } else { ans -= cal(n-1, tot); } printf("%I64d\n", (ans%mod + mod)%mod); return 0; }
C. Relay Race
很不错的dp!
给你一个矩阵,问你一只乌龟从左上角走两遍到右下角,(x,
y)只可以走到(x+1,y),(x,y+1),中间经过的所有点的值的总和最大为多少,重复走的只计算一次。
思路:现在考虑成两只不同乌龟一起同步走的状态,设dp[i][j][k]表示左边的乌龟在第i行第j列并且右边的乌龟在j+k列,由于是同步走的,所以可以计算出右边乌龟具体的位置,然后状态转移只有四种,很容易写了接下来~具体见代码~
// Author : JayYe Created Time: 2013-11-10 11:00:46 #include <stdio.h> #include <string.h> #include <algorithm> using namespace std; const int maxn = 300 + 5; const int INF = 1<<30; int a[maxn][maxn], dp[maxn][maxn][maxn]; int fix(int x1, int y1, int x2, int y2) { if(x1 == x2 && y1 == y2) return a[x1][y1]; return a[x1][y1] + a[x2][y2]; } int main() { int n; scanf("%d", &n); for(int i = 1;i <= n; i++) for(int j = 1;j <= n; j++) scanf("%d", &a[i][j]); for(int i = 1;i <= n; i++) for(int j = 1;j <= n; j++) for(int k = 0;k <= n; k++) dp[i][j][k] = -INF; dp[1][1][0] = a[1][1]; for(int i = 1;i <= n; i++) { for(int j = 1;j <= n; j++) { for(int k = 0;k <= n; k++) if(dp[i][j][k] > -INF){ int x1 = i, y1 = j, y2 = y1 + k, x2 = i + j - y2; dp[i][j+1][k] = max(dp[i][j+1][k], dp[i][j][k] + fix(i, y1+1, x2, y2+1)); if(k > 0) dp[i][j+1][k-1] = max(dp[i][j+1][k-1], dp[i][j][k] + fix(i, y1+1, x2+1, y2)); dp[i+1][j][k] = max(dp[i+1][j][k], dp[i][j][k] + fix(i+1, j, x2+1, y2)); dp[i+1][j][k+1] = max(dp[i+1][j][k+1], dp[i][j][k] + a[i+1][j] + a[x2][y2+1]); } } } printf("%d\n", dp [0]); return 0; }
相关文章推荐
- Codeforces Round #129 (Div. 1) (各种好dp)
- SRM 511 DIV1 500pt(DP)
- Codeforces Round #426 (Div. 2) D. The Bakery(线段树维护dp)
- Codeforces Round #418 (Div. 2) -- C. An impassioned circulation of affection(DP预处理)
- codeforces #260 DIV 2 C题Boredom(DP)
- Codeforces Round #396(Div. 2)C. Mahmoud and a Message【dp】
- Codeforces Round #396 (Div. 2) C题Mahmoud and a Message(dp)解题报告
- Codeforces Round #420 (Div. 2) E. Okabe and El Psy Kongroo dp+矩阵快速幂
- Codeforces Round #272 (Div. 1)D(字符串DP)
- TCSRM 593 div2(1000)(dp)
- Codeforces Round #260 (Div. 1) A. Boredom (简单dp)
- div 中如何加各种边框(转)
- srm 302 div2 1000(简单题,bfs,dp)
- 汇编入门学习笔记 (七)—— dp,div,dup
- Codeforces Round #408 (Div. 2) E. Exam Cheating 滚动dp
- CSS/html各种div 等元素居中的几种方法
- Codeforces Round #374 (Div. 2) C bfs+dp
- TopCoder SRM 666 Div2 Problem 999 - WalkOverATreeDiv2 (树形DP)
- Codeforces Round #367 (Div. 2) A(暴力) B(二分查找) C(DP) D(01字典树)
- Codeforces Round #367 (Div. 2) C. Hard problem (dp)