hdu 4963 Dividing a String 中途相遇法 2014 Multi-University Training Contest 9-1004

题意：

给定一个n*2长的由a和b组成的序列T，现在要将序列分成两个长为n的子序列T1和T2，使得T1=T2。序列中每个位置由一个价值ci，T1中的各个元素的位置价值和为s1，T2的为s2，现在要求最小的abs(s1-s2)。

题解：

一个中途相遇法应用的题目

看“爱but的苍蝇”的博客可以看到这次比赛的题解。

将序列分成两半，各n各字符。对于前一半，枚举每个字符属于T1还是T2，共2^n种情况。在一个方案中，T1的子序列为A，T2的子序列为B，我们假定A的长度不大于B的长度。显然，A是B的前缀才符合要求。令C=B-A，即B去掉A前缀后剩余的部分。同理，处理后一半得到一个C‘，就是去掉相同后缀后剩余的部分。当C=C'就是一个方案使得T1=T2。我们记录每个C和C’可以取得的权值，最后排序一下就能求得最小的权值了。复杂度O（n*2^n）

我根据上述处理，然后将符合的C和C‘存到结构体数组e中，结构体中flag标志C和C’，w表示C（C‘）部分的T2权值-T1权值，id表示C的值和长度（用len*1e7+C组成，这里特别注意，长度不同的C是不同的。）。最后根据id和w排序，计算所有相邻的C和C’的权值差，其中最小的就是答案了。

剪枝：

1）上述假设A长度不大于B的长度，那么枚举2^n的时候，可以先判断1的个数和0的个数，减少接近一半的枚举操作。

2）a的个数是奇数的时候就没办法分成两部分，可以直接特判掉。

注意：这个代码还是不够好，在G++上用9671MS过了，C++还是TLE了，哎。。

代码：

#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <queue>
#include <map>
#include <set>
#include <vector>
#include <cctype>
using namespace std;

const int maxn=44;
const int maxm=(1<<20)+10;
const int INF=1e9;
char a[maxn];
int c[maxn],b[maxn],x,y,xx,yy,sx,sy,t,f[maxn],g[maxm],len[maxn];
struct node{
    int id,w,flag;
    bool operator <(const node &a)const{
        if(id==a.id)return w<a.w;
        return id<a.id;
    }
}e[maxm*2];
int bitcount(int i)
{
    int x=0;
    while(i)
    {
        x+=(i&1);
        i=i>>1;
    }
    return x;
}
void init()
{
    for(int i=0;i<=22;i++)
    {
        f[i]=(1<<i)-1;
        len[i]=i*10000000;
    }
    for(int i=0;i<(1<<20);i++)
    {
        g[i]=bitcount(i);
    }
}
void judge()
{
    int i,j,k;

    if((y>>(yy-xx))^x)return ;
    //printf("***%d %d %d %d %d\n",x,y,xx,yy,y>>(yy-xx));
    e[t].flag=0;
    e[t].id=len[yy-xx]+(y&f[yy-xx]);
    e[t++].w=sy-sx;
}
void judge2()
{
    int i,j,k,p;
    if(x^(y&f[xx]))return ;
    e[t].flag=1;
    e[t].id=len[yy-xx]+(y>>xx);
    e[t++].w=sy-sx;
}
int main()
{
    //freopen("D:\\1004.in","r",stdin);
    //freopen("D:\\t1004.out","w",stdout);
    init();
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        if(n==0)break;
        int i,j,k=0,ans=INF;
        scanf("%s",a);
        for(i=0;i<2*n;i++)
        {
            b[i]=a[i]-'a';
            k+=b[i];
        }
        for(i=0;i<n*2;i++)
            scanf("%d",&c[i]);
        if(k%2){printf("-1\n");continue;}
        t=0;
        for(i=0;i<(1<<n);i++)
        {
            if(g[i]>n-g[i])continue;
            x=y=xx=yy=sx=sy=0;
            for(j=n-1;j>=0;j--)
            {
                if((1<<j)&i)
                {
                    x=(x<<1)|b[n-1-j];
                    sx+=c[n-1-j];
                    xx++;
                }
                else
                {
                    y=(y<<1)|b[n-1-j];
                    sy+=c[n-1-j];
                    yy++;
                }
            }
            //cout<<"**"<<xx<<" "<<yy<<endl;
            if(xx<=yy)judge();
        }
        for(i=0;i<(1<<n);i++)
        {
            if(g[i]>n-g[i])continue;
            x=y=xx=yy=sx=sy=0;
            for(j=n-1;j>=0;j--)
            {
                if((1<<j)&i)
                {
                    x=(x<<1)|b[2*n-1-j];
                    sx+=c[2*n-1-j];
                    xx++;
                }
                else
                {
                    y=(y<<1)|b[2*n-1-j];
                    sy+=c[2*n-1-j];
                    yy++;
                }
            }
            //cout<<"**"<<x<<" "<<y<<endl;
            if(xx<=yy)judge2();
        }
        sort(e,e+t);
        int p=-1,q=-1;
        for(i=0;i<t;i++)
        {
            //printf("*%d %d %d\n",e[i].flag,e[i].id,e[i].w);
            if(e[i].flag==0)p=i;
            else q=i;
            if(p==-1||q==-1)continue;
            if(e[p].id==e[q].id)ans=min(ans,abs(e[p].w-e[q].w));
        }
        if(ans==INF)printf("-1\n");
        else printf("%d\n",ans);
    }
    return 0;
}