0

コンテキストを提供するために、私はProgramming Praxis Bingo Challengeに取り組んでおり、このコードをどれだけ速く実行できるかを確認したいと考えていました。

static void fisher_yates(T& source) {   
    const size_t len = source.size();
    for(size_t i = 1; i < len;++i) {
        std::swap(source[i],source[rand() % (i+1)]);
    }
}
std::array<int,25> generate_table() {
    std::array<int,25> bingo_grid;
    for(int i = 0 ; i < 25;++i) {
        switch(i) {
        case 0: case 1: case 2: case 3: case 4:
            bingo_grid[i] = rand() % 15 + 1;
            break;
        case 5: case 6: case 7: case 8: case 9:
            bingo_grid[i] = rand() % 15 + 16;
            break;
        case 10: case 11: case 12: case 13: case 14:
            bingo_grid[i] = rand() % 15 + 31;
            break;
        case 15: case 16: case 17: case 18: case 19:
            bingo_grid[i] = rand() % 15 + 46;
            break;
        case 20: case 21: case 22: case 23: case 24:
            bingo_grid[i] = rand() % 15 + 61;
            break;
        }
    }
    bingo_grid[12] = 0;
    return bingo_grid;
}

bool is_bingoed(const std::array<int,25>& grid) { 
    // Check columns
    if(grid[0] == 0) {
        if(grid[1] == 0 && grid[2] == 0 && grid[3] == 0 && grid[4] == 0)
            return true;
        if(grid[0]  == 0 && grid[6]  == 0 && grid[18]  == 0 && grid[24]  == 0)
            return true;
        if(grid[5] == 0 && grid[10] == 0 && grid[15] == 0 && grid[20] == 0)
            return true;
    }
    if(grid[1] == 0) {
        if(grid[6]  == 0 && grid[11]  == 0 && grid[16]  == 0 && grid[21]  == 0)
            return true;
    }
    if(grid[2] == 0) {  
        if(grid[7] == 0 && grid[17]  == 0 && grid[22]  == 0)
            return true;
    }
    if(grid[3] == 0) {
        if(grid[8]  == 0 && grid[13]  == 0 && grid[18]  == 0 && grid[23]  == 0)
            return true;
    }
    if(grid[4] == 0) {
        if(grid[9]  == 0 && grid[14]  == 0 && grid[19]  == 0 && grid[24]  == 0)
            return true;
        if(grid[8]  == 0 && grid[16]  == 0 && grid[21]  == 0)
            return true;
    }
    if(grid[6] == 0) {
        if(grid[6]  == 0 && grid[7]  == 0 && grid[8]  == 0 && grid[9]  == 0)
            return true;
    }
    if(grid[12] == 0) {
        if(grid[10]  == 0 && grid[11]  == 0 && grid[13]  == 0 && grid[14] == 0)
            return true;
    }
    if(grid[18] == 0) {
        if(grid[15]  == 0 && grid[16]  == 0 && grid[17]  == 0 && grid[19]  == 0)
            return true;
    }   
    return false;
}

static bool mark_card(const int card,std::array<int,25>& bingo_grid) {
    for(auto &i : bingo_grid)
        if(card == i) {
            i = 0;
            return true;
        }
        return false;
}

int play_game() {
    // Bingo is 5 columns, each column(n) is random permutation of 1-15*n
    // Fisher-Yates to generate random permutations

    // Create 500 playing cards
    const int max = 500;
    std::vector<std::array<int,25>> bingo_cards;
    bingo_cards.reserve(max);
    for(int i = 0; i<max;++i) {
        bingo_cards.push_back(generate_table());
        //display_bingo(bingo_cards[i]);
    }
    // Random shuffle 75 cards
    auto iter = boost::counting_range(1,76);
    std::vector<int> cards(std::begin(iter),std::end(iter));
    fisher_yates(cards);
    bool is_finished = false;
    int counter = 0;
    for(auto card : cards) {
        for(auto& playing_card : bingo_cards) {
            if(mark_card(card,playing_card)) {
                //display_bingo(playing_card);
                if(is_bingoed(playing_card))
                    return counter;
            }
        }
        counter++;
    }
    return counter;
}
int bingo() {
    srand(time(NULL));
    int total = 0;
    for(int i = 0 ; i < 10000;i++) {
        total+=play_game();
    }
    boost::singleton_pool<boost::pool_allocator_tag, sizeof(int)>::release_memory();
    return total / 10000;
}

元のバージョンでは、boost::multi_array を使用してグリッドを表現していました。プロファイリングの後、それを std::array に変更したところ、大幅に高速化されました。次に、fisher_yates shuffle を使用してビンゴ カードを生成する方法から、乱数ジェネレーターを使用する方法に移行しました。最後に、is_bingoed テスト関数を変更して、呼び出しごとのチェック回数を減らし、ゲーム オーバー チェックを高速化しました。

これはすべて役に立ちました。現在、このコードをプロファイリングすると、generate_table 関数が 72% の時間を占め、mark_card() が 18%、is_bingoed() が約 6% です。いずれかの速度を向上させるために何ができるかを確認するためのヒントを探しています。

is_bingoed() で最初に考えたのは、SSE 組み込み関数を使用して 0 と比較することです (おそらく XOR を使用しますか?) が、generate_table() または mark_car() についてのアイデアはありません。これは楽しみのための自己挑戦ですが、他の人はどう思うでしょうか?

現在のタイミングは、2Ghz Q6660 で 4.6 秒かかります (元の 35 秒から短縮)。

4

2 に答える 2

1

最も高価な関数 に注目するだけgenerate_tableで、コードのこの部分を単純化し、分岐を少なくすることができます。これが役立つ場合があります。

for(int i = 0 ; i < 25;++i) {
    switch(i) {
    case 0: case 1: case 2: case 3: case 4:
        bingo_grid[i] = rand() % 15 + 1;
        break;
    case 5: case 6: case 7: case 8: case 9:
        bingo_grid[i] = rand() % 15 + 16;
        break;
    case 10: case 11: case 12: case 13: case 14:
        bingo_grid[i] = rand() % 15 + 31;
        break;
    case 15: case 16: case 17: case 18: case 19:
        bingo_grid[i] = rand() % 15 + 46;
        break;
    case 20: case 21: case 22: case 23: case 24:
        bingo_grid[i] = rand() % 15 + 61;
        break;
    }
}

例えば

for(int i = 0 ; i < 25;++i) {
    int r = rand() % 15 + 1;
    bingo_grid[i] = r + (i / 5) * 15;
}

それを超えて、より高速な rand() を見て、除算とモジュロを取り除くことができるかどうかも確認します。

別の注意として、bingo_grid で重複する数字を防ぐものがないという点で、アルゴリズムに欠陥がある可能性があります。

于 2013-06-26T08:02:35.357 に答える
0

is_bingoed() メソッドを SSE 命令 (Agner Fog のライブラリを使用) を使用するように変更し、Paul R の generate_table() を使用すると、時間がわずか 1.05 秒に短縮されました。Intel の fast_rand()関数を使用すると、0.38 秒に短縮されました。そこで、興味のある他の人のためにコードの変更を貼り付けようと思いました。

static unsigned int g_seed;


//Used to seed the generator.
inline void fast_srand( int seed )
{
    g_seed = seed;
}


//fastrand routine returns one integer, similar output value range as C lib.
inline int fastrand()
{
    g_seed = (214013*g_seed+2531011);
    return (g_seed>>16)&0x7FFF;

}
bool is_bingoed(const std::array<int,25>& grid) { 
    // Check columns
    Vec8i vec(grid[0],grid[1],grid[2],grid[3],grid[4],0,0,0);
    Vec4i vec2(grid[6],grid[18],grid[24],0);
    Vec4i vec3(grid[5],grid[10],grid[15],20);
    Vec8i vec4(grid[1],grid[6],grid[11],grid[16],grid[21],0,0,0);
    Vec4i vec5(grid[2],grid[7],grid[17],grid[22]);
    Vec8i vec6(grid[3],grid[8],grid[13],grid[18],grid[23],0,0,0);
    Vec8i vec7(grid[4],grid[9],grid[14],grid[19],grid[24],0,0,0);
    Vec4i vec8(grid[8],grid[16],grid[21],grid[4]);
    Vec4i vec9(grid[6],grid[7],grid[8],grid[9]);
    Vec8i vec10(grid[12],grid[10],grid[11],grid[13],grid[14],0,0,0);
    Vec8i vec11(grid[18],grid[15],grid[16],grid[17],grid[19],0,0,0);
    if(horizontal_and(vec) && horizontal_and(vec2) && horizontal_and(vec3) && horizontal_and(vec4) &&
        horizontal_and(vec5) && horizontal_and(vec6) && horizontal_and(vec7) && horizontal_and(vec8)) {
            return false;
    }
    if(horizontal_and(vec9) && horizontal_and(vec10) && horizontal_and(vec11)) {
            return false;
    }
    return true;
}
于 2013-07-11T22:13:58.943 に答える