APIO2015

还没有写完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 ;
}