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
*/

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章