HDU 4990 Ordered Subsequence --数据结构优化DP
生活随笔
收集整理的這篇文章主要介紹了
HDU 4990 Ordered Subsequence --数据结构优化DP
小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.
題意:給一串?dāng)?shù)字,問長度為m的嚴(yán)格上升子序列有多少個
解法:首先可以離散化為10000以內(nèi),再進(jìn)行dp,令dp[i][j]為以第i個元素結(jié)尾的長度為j的上升子序列的個數(shù),
則有dp[i][j] = SUM(dp[k][j-1]) ?(a[k] < a[i] && k < i)
不可能直接遍歷,所以考慮優(yōu)化,可以看出dp方程相當(dāng)于一個區(qū)間求和,所以可以用樹狀數(shù)組來優(yōu)化。
代碼:
#include <iostream> #include <cmath> #include <iostream> #include <cstdio> #include <cstring> #include <cstdlib> #include <cmath> #include <algorithm> #define SMod 123456789 #define lll __int64 using namespace std; #define N 10007lll c[N],a[N],b[N]; int n,m; lll dp[N][104];int lowbit(int x) {return x & (-x); }void modify(int pos,lll val) {while(pos <= n){c[pos] += val;pos += lowbit(pos);} }lll getsum(int pos) {lll res = 0;while(pos > 0){res = (res+c[pos])%SMod;pos -= lowbit(pos);}return res; }int main() {int i,j;while(scanf("%d%d",&n,&m)!=EOF){for(i=1;i<=n;i++){scanf("%I64d",&a[i]);b[i] = a[i];}sort(b+1,b+n+1);memset(dp,0,sizeof(dp));for(i=1;i<=n;i++)dp[i][1] = 1;for(j=2;j<=m;j++){memset(c,0,sizeof(c));for(i=1;i<=n;i++){int ind = lower_bound(b+1,b+n+1,a[i])-b;dp[i][j] = getsum(ind-1);modify(ind,dp[i][j-1]);}}lll ans = 0;for(i=1;i<=n;i++)ans = (ans + dp[i][m])%SMod;printf("%I64d\n",ans);}return 0; } View Code?
轉(zhuǎn)載于:https://www.cnblogs.com/whatbeg/p/3961018.html
總結(jié)
以上是生活随笔為你收集整理的HDU 4990 Ordered Subsequence --数据结构优化DP的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: xcap 发包工具 使用帮助
- 下一篇: 暴力破解-----token验证