解析
比较妙的一道题
首先的一个技巧是把一条线段拆成左和右两个端点
顺便也就有了离散化
利用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
*/