CF1313D：Happy New Year（状压dp）

解析

比较妙的一道题

首先的一个技巧是把一条线段拆成左和右两个端点
顺便也就有了离散化
利用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_JUDGE
  freopen("a.in","r",stdin);
  freopen("a.out","w",stdout);
  #endif
  n=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
*/