補寫一下

poj3171 設f[i]表示覆蓋L~i的最小花費，把區間按左端點排序，枚舉區間，f[a[i].r]=min{f[a[i].l~(a[top].r-1)]}+a[i].c (當然還要和原值比較的)

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;struct node{int l,r,d;}a[31000];
bool cmp(node n1,node n2){return n1.r<n2.r;}struct trnode
{int l,r,lc,rc,c;
}tr[210000];int trlen;
void bt(int l,int r)
{int now=++trlen;tr[now].l=l;tr[now].r=r;tr[now].c=(1<<30);tr[now].lc=tr[now].rc=-1;if(l<r){int mid=(l+r)/2;tr[now].lc=trlen+1;bt(l,mid);tr[now].rc=trlen+1;bt(mid+1,r);}
}
void change(int now,int p,int k)
{if(tr[now].l==tr[now].r){tr[now].c=k;return ;}int lc=tr[now].lc,rc=tr[now].rc;int mid=(tr[now].l+tr[now].r)/2;if(p<=mid)change(lc,p,k);else      change(rc,p,k);tr[now].c=min(tr[lc].c,tr[rc].c);
}
int getmin(int now,int l,int r)
{if(tr[now].l==l&&tr[now].r==r)return tr[now].c;int lc=tr[now].lc,rc=tr[now].rc;int mid=(tr[now].l+tr[now].r)/2;if(r<=mid)  return getmin(lc,l,r);else if(mid+1<=l)return getmin(rc,l,r);else return min(getmin(lc,l,mid),getmin(rc,mid+1,r));
}int f[110000];
int main()
{int n,L,R;scanf("%d%d%d",&n,&L,&R);L++,R++;for(int i=1;i<=n;i++){scanf("%d%d%d",&a[i].l,&a[i].r,&a[i].d);a[i].l++,a[i].r++;if(a[i].r<L||R<a[i].l)i--,n--;if(a[i].l<L)a[i].l=L;if(R<a[i].r)a[i].r=R;}sort(a+1,a+n+1,cmp);int top=1; trlen=0,bt(L,R);memset(f,63,sizeof(f));for(int i=L;i<=R;i++){while(top<=n&&a[top].r==i){if(a[top].l==L){if(a[top].d<f[i])f[i]=a[top].d, change(1,i,f[i]);}else {int k=getmin(1,a[top].l-1,a[top].r-1)+a[top].d;if(k<f[i])f[i]=k, change(1,i,f[i]);}top++;}}if(f[R]==f[0])printf("-1\n");else printf("%d\n",f[R]);return 0;
}

poj3171

hdu5542 f[i][j]表示枚舉到第幾個位置，長度為j的嚴格上升序列有多少（注意一定包括第i個位置）

f[i][j]=∑(k<j&&a[k]<a[j])f[k][j-1] 這里先離散化，然后開m個樹狀數組記錄，想想cdq分治，時間維有序第一個限制不管，第二個用樹狀數組解決。實際操作中，是不用定義數組的，直接插入樹狀數組中相應位置即可。

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
const int mod=1e9+7;int lslen;int s[1100][1100];
int lowbit(int x){return x&-x;}
void change(int x,int w,int k)
{while(x<=lslen){s[w][x]=(s[w][x]+k)%mod;x+=lowbit(x);}
}
int getsum(int x,int w)
{int ret=0;while(x>0){ret=(ret+s[w][x])%mod;x-=lowbit(x);}return ret;
}//----------------bit---------------------int a[1100],ls[1100];
int main()
{int T,T_T=0;scanf("%d",&T);while(T--){int n,m;scanf("%d%d",&n,&m);lslen=0;for(int i=1;i<=n;i++)scanf("%d",&a[i]), ls[++lslen]=a[i];sort(ls+1,ls+lslen+1);lslen=unique(ls+1,ls+lslen+1)-ls-1;for(int i=1;i<=n;i++)a[i]=lower_bound(ls+1,ls+lslen+1,a[i])-ls;memset(s,0,sizeof(s));for(int i=1;i<=n;i++){change(a[i],1,1);int li=min(i,m);for(int j=2;j<=li;j++)change(a[i],j,getsum(a[i]-1,j-1));}printf("Case #%d: %d\n",++T_T,getsum(lslen,m));}return 0;
}

hdu5542

?

轉載于:https://www.cnblogs.com/AKCqhzdy/p/9468658.html

總結

以上是生活随笔為你收集整理的0x58 数据结构优化DP的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。