hdu 4712 2种解发（状态dp| 随机化） 我又相信爱情了~~(第一发)

  这个题目我们省赛的时候就已经做看过一次，但是那时候知道你n^2 一定是超时的，然后就gg了，这次在warm up上是发现了状态dp 然后估计了下复杂度大概是（2^20*20吧，然后就开始写了后面发现只向后面推冒失得不到答案，然后果断机智选择spfa的能向前个状态，然后发现这个时间是硬伤，第一组数据差不多2s才出来。这能交嘛？然后自己yy出向后退也可以，但是要记录每一个点最短和次短长度，对于已经有点最短明显是0 次短是inf 后面重每个地方推就可以了判断下是不是来自同一个父亲，然后就是tel了，然后就是果断的各种改小参数了，然后过了好久才ac 哎。。。。。

#include<cstdio>
#include<cstring>
#include<queue>
#include<iostream>
#define inf 0x3f3f3f3f
#define M 1100000
using namespace std;
int is_hash[M];
int tmp;
int fa[M][2];
int dp[M][2];
int big[M];
inline int change(char p1){
    if(p1>='0'&&p1<='9') return p1-'0';
    return p1-'A'+10;
}
int main(){
    int cas;
    big[0]=0;
    big[1]=1;
    for(int i=2;i<(1<<20);i++){
       big[i]=i&(-i);
    }
    while(~scanf("%d",&cas)){
        while(cas--){
            int flag=0;
            int n;scanf("%d",&n);
            char str[10];
            memset(dp,inf,sizeof(dp));
            memset(fa,-1,sizeof(fa));
            for(int i=0;i<n;i++){
                scanf("%s",str);
                int t1=0;
                t1=(change(str[0])<<16)+(change(str[1])<<12)+(change(str[2])<<8)+(change(str[3])<<4)+(change(str[4]));
                dp[t1][0]=0;
                fa[t1][0]=t1;
            }
            if(flag){
                puts("0");
                continue;
            }
            else{
                int tmin=inf;
                int k=0;
                int next=20;
                for(int i=(1<<20)-1;i>=0;i--){
                    if(dp[i][0]>=tmin) continue;
                    int o2=i;
                    int o1;
                    int o;
                    while(o2){
                        int o1=big[o2];
                            o2=o2^o1;
                        o=i^o1;
                        if(fa[o][0]!=fa[i][0]){
                            if(dp[o][1]>dp[i][0]+1){
                                dp[o][1]=dp[i][0]+1;
                                fa[o][1]=fa[i][0];
                            }
                            if(dp[o][1]<dp[o][0]){
                                    tmp=dp[o][1];
                                    dp[o][1]=dp[o][0];
                                    dp[o][0]=tmp;
                                    tmp=fa[o][1];
                                    fa[o][1]=fa[o][0];
                                    fa[o][0]=tmp;
                            }
                        }
                    }
                    if(dp[i][0]+dp[i][1]<tmin){
                        tmin=dp[i][0]+dp[i][1];
                    }
                }
                printf("%d\n",tmin);
            }
        }
    }
    return 0;
}
/*
5
1
0000B
*/

后面听被人说这个随机算法，然后发现真的能过啊，太坑了，但是听了大牛的讲解后发现答案真的分布的很大，发现比dp块不知道多少倍啊，看来还真有随机算法能过的题目啊，哎，长见识了，因为rand要从0-99999 自己改写了下，然后果断过了

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<ctime>
#define M 1700000
using namespace std;
int first[M];
int get_long_rand(int n){
    int t1=rand()%10000;
    int t2=rand()%10000;
    return (t1*10000+t2)%n;
}
void init(){
    for(int i=0;i<M;i++){
        int t1=i;
        int tsum=0;
        while(t1){
            if(t1&1) tsum++;
            t1>>=1;
        }
        first[i]=tsum;
    }
}
int f(char p){
    if(p>='0'&&p<='9') return p-'0';
    return p-'A'+10;
}
int a[M];
int main(){
    srand(time(0));
    init();
    int cas;while(~scanf("%d",&cas)){
        while(cas--){
            int n;scanf("%d",&n);
            for(int i=0;i<n;i++){
                char str[10];
                scanf("%s",str);
                int t1=0;
                t1+=(f(str[0])<<16);
                t1+=(f(str[1])<<12);
                t1+=(f(str[2])<<8);
                t1+=(f(str[3])<<4);
                t1+=(f(str[4]));
                a[i]=t1;
            }
            int tsum=0;
            int mark=60000;
            int tmin=20;
            if(n>50000) tmin=4;
            while(mark--){
                int t1=get_long_rand(n);
                int t2=get_long_rand(n);
                if(t1==t2) continue;
                t1=a[t1];t2=a[t2];
                if(first[t1^t2]<tmin) tmin=first[t1^t2];
            }
            printf("%d\n",tmin);
        }
    }
    return 0;
}
/*
2
2
00001
00002
*/