• Bzoj3926: [Zjoi2015]诸神眷顾的幻想乡


    题面

    传送门

    Sol

    求多个串的不同子串的个数

    广义后缀自动机
    也就是可以表示所有的串的所有的后缀的自动机

    那么每次建一个串后,另外一个串接在初始节点下面建就好了

    叶子节点最多(20)
    那么对于每个叶子结点遍历一遍树建立(sam)
    注意回溯时要把(last)指回来

    # include <bits/stdc++.h>
    # define RG register
    # define IL inline
    # define Fill(a, b) memset(a, b, sizeof(a))
    using namespace std;
    typedef long long ll;
     
    IL int Input(){
        RG int x = 0, z = 1; RG char c = getchar();
        for(; c < '0' || c > '9'; c = getchar()) z = c == '-' ? -1 : 1;
        for(; c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) + (c ^ 48);
        return x * z;
    }
     
    const int maxn(1e5 + 5);
    const int maxm(4e6 + 5);
    
    int trans[10][maxm], fa[maxm], tot = 1, last, len[maxm];
    int n, first[maxn], cnt, col[maxn], degree[maxn];
    ll ans;
    
    struct Edge{
    	int to, next;
    } edge[maxn << 1];
    
    IL void Add(RG int u, RG int v){
    	edge[cnt] = (Edge){v, first[u]}, first[u] = cnt++;
    }
    
    IL void Extend(RG int c, RG int p){
    	RG int np = ++tot;
    	len[last = np] = len[p] + 1;
    	while(p && !trans[c][p]) trans[c][p] = np, p = fa[p];
    	if(!p) fa[np] = 1;
    	else{
    		RG int q = trans[c][p];
    		if(len[q] == len[p] + 1) fa[np] = q;
    		else{
    			RG int nq = ++tot;
    			fa[nq] = fa[q], len[nq] = len[p] + 1;
    			for(RG int i = 0; i < 10; ++i) trans[i][nq] = trans[i][q];
    			fa[q] = fa[np] = nq;
    			while(p && trans[c][p] == q) trans[c][p] = nq, p = fa[p];
    		}
    	}
    }
    
    IL void Dfs(RG int u, RG int ff, RG int p){
    	Extend(col[u], p), p = last;
    	for(RG int e = first[u]; e != -1; e = edge[e].next)
    		if(edge[e].to != ff) Dfs(edge[e].to, u, p);
    }
    
    int main(RG int argc, RG char *argv[]){
    	n = Input(), Input();
    	for(RG int i = 1; i <= n; ++i) first[i] = -1, col[i] = Input();
    	for(RG int i = 1; i < n; ++i){
    		RG int u = Input(), v = Input();
    		Add(u, v), Add(v, u), ++degree[u], ++degree[v];
    	}
    	for(RG int i = 1; i <= n; ++i)
    		if(degree[i] == 1) last = 1, Dfs(i, 0, 1);
    	for(RG int i = 1; i <= tot; ++i) ans += len[i] - len[fa[i]];
    	printf("%lld
    ", ans);
        return 0;
    }
    
  • 相关阅读:
    python 第三方库大全
    windows CMD实现的信息收集工具
    ip协议是哪一层的协议
    MetaWeblog访问地址
    通过卸载程序触发恶意进程
    IP地址查询接口
    mysql 密码忘记解决办法
    查询IP地址的免费API
    showdan
    【译】android的审计和hacking工具
  • 原文地址:https://www.cnblogs.com/cjoieryl/p/8929095.html
Copyright © 2020-2023  润新知