还没有写完APIO2015的题目，打算今天写一写。
T1：
按位DP，DP时要保证已确定的位为0。
前4组设f[n][k]表示把前n个分成k组是否合法。
最后一组设g[n]表示把前n个最少分为多少组才符合题意。
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cctype>
#define rep(i,s,t) for(int i=s;i<=t;i++)
#define ren for(int i=first[x];i!=-1;i=next[i])
using namespace std;
inline int read() {
int x=,f=;char c=getchar();
for(;!isdigit(c);c=getchar()) if(c=='-') f=-;
for(;isdigit(c);c=getchar()) x=x*+c-'';
return x*f;
}
const int maxn=;
typedef long long ll;
int n,A,B,f[][],g[maxn];
ll t,S[maxn];
int check(int x) {
if(A==) {
g[]=;
rep(i,,n) {
g[i]=<<;
rep(j,,i-) if(!((S[i]-S[j])&t)) g[i]=min(g[i],g[j]+);
}
return g[n]<=B;
}
memset(f,,sizeof(f));
f[][]=;
rep(i,,n) rep(k,,B) rep(j,,i-) if(!((S[i]-S[j])&t)&&f[j][k-]) f[i][k]=;
rep(i,A,B) if(f[n][i]) return ;
return ;
}
int main() {
n=read();A=read();B=read();
rep(i,,n) S[i]=S[i-]+read();
ll ans=;
for(int i=;i>=;i--) {
t^=(1ll<<i);
if(!check(i)) ans^=(1ll<<i),t^=(1ll<<i);
}
printf("%lld\n",ans);
return ;
}
T3：
k=1时排序求中位数。
k=2时设垮桥的所有两点坐标分别为x1i,x2i。设桥修在了p1,p2的位置，对于每对点(x1,x2)，当(x1+x2)/2接近p1时肯定走p1，(x1+x2)/2接近p2时肯定走p2。故将所有点对按x1+x2升序排序，这样肯定存在一条分割线，使得分割线左侧的走左边一座桥，分割线右侧的走右边一座桥。这样我们实现一个可以支持插入和查询中位数的数据结构即可。Treap什么的就行了，但我DB地用了动态开节点的线段树。
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cctype>
#define rep(s,t) for(int i=s;i<=t;i++)
#define ren for(int i=first[x];i!=-1;i=next[i])
using namespace std;
inline int read() {
int x=,f=;char c=getchar();
for(;!isdigit(c);c=getchar()) if(c=='-') f=-;
for(;isdigit(c);c=getchar()) x=x*+c-'';
return x*f;
}
inline char getc() {
char c=getchar();
while(!isalpha(c)) c=getchar();
return c;
}
const int maxn=;
typedef long long ll;
struct Point {
int x1,x2;
bool operator < (const Point& ths) const {return x1+x2<ths.x1+ths.x2;}
}A[maxn];
int k,n,tmp[maxn];
ll ans,mn=(1ll<<),f[maxn],g[maxn],sumv[];
int ls[],rs[],s[],ToT,root;
void insert(int& o,int l,int r,int pos) {
if(!o) o=++ToT,s[o]=ls[o]=rs[o]=sumv[o]=;
s[o]++;sumv[o]+=pos;if(l==r) return;
int mid=l+r>>;
if(pos<=mid) insert(ls[o],l,mid,pos);
else insert(rs[o],mid+,r,pos);
}
int kth(int o,int l,int r,int k) {
if(l==r) return l;
int k2=s[ls[o]],mid=l+r>>;
if(k2>=k) return kth(ls[o],l,mid,k);
return kth(rs[o],mid+,r,k-k2);
}
ll tot,sum;
void query(int o,int l,int r,int ql,int qr) {
if(!o) return;
if(ql<=l&&r<=qr) {
tot+=s[o];sum+=sumv[o];
return;
}
int mid=l+r>>;
if(ql<=mid) query(ls[o],l,mid,ql,qr);
if(qr>mid) query(rs[o],mid+,r,ql,qr);
}
ll query(int o) {
if(!o) return ;
int v=kth(o,,,s[o]/);
ll ret=;
tot=sum=;query(o,,,,v);ret+=tot*v-sum;
tot=sum=;query(o,,,v,);ret-=tot*v-sum;
return ret;
}
int main() {
k=read();n=read();int t=;
rep(,n) {
char t1,t2;int x1,x2;
t1=getc();x1=read();t2=getc();x2=read();
if(t1==t2) ans+=abs(x1-x2);
else A[++t]=(Point){x1,x2};
}
n=t;sort(A+,A+n+);ans+=n;
if(k==) {
t=;
rep(,n) tmp[t++]=A[i].x1,tmp[t++]=A[i].x2;
sort(tmp,tmp+t);
rep(,t-) ans+=abs(tmp[t/]-tmp[i]);
printf("%lld\n",ans);
return ;
}
ToT=root=;
rep(,n) {
insert(root,,,A[i].x1);insert(root,,,A[i].x2);
f[i]=query(root);
}
ToT=root=;
for(int i=n;i;i--) {
g[i]=query(root);
insert(root,,,A[i].x1);insert(root,,,A[i].x2);
}
rep(,n) mn=min(mn,f[i]+g[i]);if(!n) mn=;
printf("%lld\n",ans+mn);
return ;
}