ACM模版

描述

題解

AC 自動機 + 高斯消元。

程式碼

#include <queue>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;

const int MAX_LEN = 1010;
const int MAXN = 120;
const int MAXM = 6;
const int MAXF = 15
 
;
const double EPS = 1e-6;
const double MAGIC = 1.0 / 6;

int n, L, sz;
int fp[MAXF];
int tag[MAXN];
int fail[MAXN];
int trie[MAXN][MAXM];
queue<int> qi;

/* *********** Aho *********** */
struct Aho
{
    int root;

    int newnode()                           //  靜態建立新節點
    {
        memset(trie[sz], -1
 
, sizeof(trie[sz]));
        tag[sz] = 0;
        sz++;
        return sz - 1;
    }

    void init()
    {
        sz = 0;
        newnode();
    }

    void insert(int s[], int id)            //  插入字串構建ac自動機，構建trie樹
    {
        int len = L, p = 0;
        for (int i = 0; i < len; i++)
        {
            int
 
 _id = s[i] - 1;
            if (trie[p][_id] == -1)
            {
                trie[p][_id] = newnode();
            }
            p = trie[p][_id];
        }
        tag[p] = id;                        //  結束標記
        fp[id] = p;
    }

    void getfail()                          //  構建自動機fail指標
    {
        while (!qi.empty())
        {
            qi.pop();
        }

        fail[root] = root;                  //  root指向root
        for (int i = 0; i < MAXM; i++)
        {
            if (trie[root][i] == -1)        //  第一個字元不存在，指向root
            {
                trie[root][i] = root;
            }
            else                            //  第一個字元的fail指標指向root
            {
                fail[trie[root][i]] = root;
                qi.push(trie[root][i]);     //  並放入佇列，待bfs擴充套件
            }
        }

        while (!qi.empty())
        {
            int u = qi.front();             //  取擴充套件節點
            qi.pop();

            if (tag[fail[u]])
            {
                tag[u] = tag[fail[u]];      //  如果之前是tag，直接標記
            }
            for (int i = 0; i < MAXM; i++)  //  遍歷所有子節點
            {
                if (trie[u][i] == -1)       //  如果不存在，則子節點直接指向fail[u]節點的對應子節點
                {
                    trie[u][i] = trie[fail[u]][i];
                }
                else                        //  如果存在，則該節點的fail指標指向fail[u]節點對應的子節點
                {
                    fail[trie[u][i]] = trie[fail[u]][i];
                    qi.push(trie[u][i]);    //  繼續擴充套件
                }
            }
        }
    }
} aho;

/* *********** Guass *********** */
int equ, var;
double x[MAXN];
double a[MAXN][MAXN];

int Gauss()
{
    int max_r;
    for (int k = 0, col = 0; k < equ && col < var; k++, col++)
    {
        max_r = k;
        for (int i = k + 1; i < equ; i++)
        {
            if (fabs(a[i][col]) > fabs(a[max_r][col]))
            {
                max_r = i;
            }
        }

        if (fabs(a[max_r][col]) < EPS)
        {
            return 0;
        }

        if (k != max_r)
        {
            for (int j = col; j < var; j++)
            {
                swap(a[k][j], a[max_r][j]);
            }
            swap(x[k], x[max_r]);
        }

        x[k] /= a[k][col];
        for (int j = col + 1; j < var; j++)
        {
            a[k][j] /= a[k][col];
        }
        a[k][col] = 1;
        for (int i = 0; i < equ; i++)
        {
            if (i != k)
            {
                x[i] -= x[k] * a[i][k];
                for (int j = col + 1; j < var; j++)
                {
                    a[i][j] -= a[k][j] * a[i][col];
                }
                a[i][col] = 0;
            }
        }
    }
    return 1;
}

void init_guass(int n, int m)
{
    memset(x, 0, sizeof(x));
    memset(a, 0, sizeof(a));
    var = n;
    equ = m;
}

int son[MAXN];

void build()        // 建立方程組
{
    init_guass(sz, sz);

    for (int i = 0; i < sz; i++)
    {
        if (!tag[i])
        {
            for (int j = 0; j < MAXM; j++)
            {
                a[trie[i][j]][i] += MAGIC;
            }
        }
        a[i][i] += -1;
    }

    x[0] += -1;     //  虛擬節點以1的概率轉移到節點0
}

int name[MAXF][MAXF];

int main()
{
    int T;
    cin >> T;

    while (T--)
    {
        cin >> n >> L;

        aho.init();

        for (int i = 0; i < n; i++)
        {
            for (int j = 0; j < L; j++)
            {
                scanf("%d", name[i] + j);
            }
            aho.insert(name[i], i + 1);
        }

        aho.getfail();

        build();
        Gauss();

        for (int i = 1; i <= n; i++)
        {
            printf("%.6f%c", x[fp[i]], i == n ? '\n' : ' ');
        }
    }

    return 0;
}