解析

比較妙的一道題

首先的一個技巧是把一條線段拆成左和右兩個端點
順便也就有了離散化
利用k的關鍵限制，壓縮每個點覆蓋線段的狀態
需要動態更新每一位所對應的線段是哪條

代碼

#include<bits/stdc++.h> using namespace std; #define ll long long const int N=2e5+100; const int mod=2009; ll read(){ll x=0,f=1;char c=getchar();while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}while(isdigit(c)){x=x*10+(c^48);c=getchar();}return x*f; } inline void Max(int &x,int y){if(x<y) x=y;} int n,m,k; int bit[260],mi[9]; int dp[N][260],tot; struct node{int pl,op,id; }p[N]; bool cmp (node a,node b){if(a.pl!=b.pl) return a.pl<b.pl;else return a.op<b.op; } int nam[9]; int main(){#ifndef ONLINE_JUDGEfreopen("a.in","r",stdin);freopen("a.out","w",stdout);#endifn=read();m=read();k=read();mi[0]=1;for(int i=1;i<=k;i++) mi[i]=mi[i-1]<<1;for(int i=1;i<mi[k];i++) bit[i]=bit[i-(i&-i)]+1;for(int i=1;i<=n;i++){int l=read(),r=read();p[++tot]=(node){l,1,i};p[++tot]=(node){r+1,-1,i};}sort(p+1,p+1+tot,cmp);memset(dp,-0x3f,sizeof(dp));memset(nam,-1,sizeof(nam));dp[0][0]=0;p[tot+1].pl=p[tot].pl;for(int i=1;i<=tot;i++){if(p[i].op==1){int o(-1);for(int j=0;j<k;j++){if(nam[j]==-1){nam[j]=p[i].id;o=j;break;}}//assert(o!=-1);//printf("i=%d id=%d\n",p[i].id,o);for(int s=0;s<mi[k];s++){if(s&mi[o]) continue;Max(dp[i][s],dp[i-1][s]+(bit[s]&1)*(p[i+1].pl-p[i].pl));Max(dp[i][s|mi[o]],dp[i-1][s]+(bit[s|mi[o]]&1)*(p[i+1].pl-p[i].pl));}}else{int o(-1);for(int j=0;j<k;j++){if(nam[j]==p[i].id){nam[j]=-1;o=j;break;}}for(int s=0;s<mi[k];s++){if(s&mi[o]) Max(dp[i][s-mi[o]],dp[i-1][s]+(bit[s^mi[o]]&1)*(p[i+1].pl-p[i].pl));else Max(dp[i][s],dp[i-1][s]+(bit[s]&1)*(p[i+1].pl-p[i].pl));}}//printf("\ni=%d pl=%d id=%d op=%d\n",i,p[i].pl,p[i].id,p[i].op);//for(int s=0;s<mi[k];s++) printf(" s=%d dp=%d\n",s,dp[i][s]);}//fprintf(stderr,"check=%d\n",dp[4][6]);printf("%d\n",dp[tot][0]); } /* 3 50000 013 002 001 */

總結

以上是生活随笔為你收集整理的CF1313D：Happy New Year（状压dp）的全部內容，希望文章能夠幫你解決所遇到的問題。

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