hihocoder1457

http://hihocoder.com/problemset/problem/1457

找不重复子串的和

topo序搞一搞，用父亲更新儿子节点的val，记得乘上节点数

//#pragma comment(linker, "/stack:200000000")

//#pragma GCC optimize("Ofast,no-stack-protector")

//#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")

//#pragma GCC optimize("unroll-loops")

#include<bits/stdc++.h>

#define fi first

#define se second

#define mp make_pair

#define pb push_back

#define pi acos(-1.0)

#define ll long long

#define mod (1000000007)

#define C 0.5772156649

#define ls l,m,rt<<1

#define rs m+1,r,rt<<1|1

#define pil pair<int,ll>

#define pii pair<int,int>

#define ull unsigned long long

#define base 1000000000000000000

#define fio ios::sync_with_stdio(false);cin.tie(0)

using namespace std;

const double g=10.0,eps=1e-;

const int N=+,maxn=+,inf=0x3f3f3f3f,INF=0x3f3f3f3f3f3f3f3f;

inline void sub(ll &x){if(x>=mod)x-=mod;}

char s[N];

int c[N],a[N];

ll val[N],sum[N];

struct SAM{

    int last,cnt;

    int ch[N<<][],fa[N<<],l[N<<];

    int mx[N<<],tmp[N<<];

    void ins(int c)

    {

        int p=last,np=++cnt;last=np;l[np]=l[p]+;

        for(;p&&!ch[p][c];p=fa[p])ch[p][c]=np;

        if(!p)fa[np]=;

        else

        {

            int q=ch[p][c];

            if(l[p]+==l[q])fa[np]=q;

            else

            {

                int nq=++cnt;l[nq]=l[p]+;

                memcpy(ch[nq],ch[q],sizeof(ch[q]));

                fa[nq]=fa[q];fa[q]=fa[np]=nq;

                sum[nq]=1ll*c;

                for(;ch[p][c]==q;p=fa[p])ch[p][c]=nq;

            }

        }

        sum[np]=1ll*c;

    }

    SAM()

    {

        cnt=;

    }

    void build(){

        int len=strlen(s+);

        last=;

        for(int i=;i<=len;i++)ins(s[i]-'');

//        topo();

//        for(int i=1;i<=cnt;i++)mx[i]=l[i],tmp[i]=0;

    }

    void topo(){

        for(int i=;i<=cnt;i++)c[l[i]]++;

        for(int i=;i<=cnt;i++)c[i]+=c[i-];

        for(int i=;i<=cnt;i++)a[c[l[i]]--]=i;

    }

    void debug()

    {

        puts("");

        for(int i=;i<=cnt;i++)

        {

            printf("%d ",i);

            for(int j=;j<;j++)

                printf("%d ",ch[i][j]);

            printf("%d %d %lld %lld\n",fa[i],l[i],val[i],sum[i]);

        }

    }

    void faupdate()

    {

        for(int i=cnt;i;i--)

            tmp[fa[a[i]]]=max(tmp[fa[a[i]]],tmp[a[i]]);

    }

    void maupdate()

    {

        for(int i=;i<=cnt;i++)

            mx[i]=min(mx[i],tmp[i]),tmp[i]=;

    }

    void match()

    {

        int len=strlen(s+),now=,res=;

        for(int i=;i<=len;i++)

        {

            if(ch[now][s[i]-'a'])

            {

                now=ch[now][s[i]-'a'];

                res++;

            }

            else

            {

                while(now&&!ch[now][s[i]-'a'])now=fa[now];

                if(!now)now=,res=;

                else

                {

                    res=l[now]+;

                    now=ch[now][s[i]-'a'];

                }

            }

            tmp[now]=max(tmp[now],res);

        }

        faupdate();

        maupdate();

    }

    int maxlen()

    {

        int ans=;

        for(int i=;i<=cnt;i++)ans=max(ans,mx[i]);

        return ans;

    }

    void cal()

    {

        topo();

        for(int i=;i<=cnt;i++)

        {

            int p=a[i];

            val[p] =  val[p] + 1ll*sum[p] * (l[p] - l[fa[p]])%mod;

            sub(val[p]);

            for(int j=;j<;j++)

            {

                if(ch[p][j])

                {

                    val[ch[p][j]]=val[ch[p][j]]+val[p]*10ll%mod;

                    sub(val[ch[p][j]]);

                }

            }

         //   cout << val[p] << endl;

        }

       // debug();

        ll ans=;

        for(int i=cnt;i;i--)

        {

            ans=ans+val[i];

            sub(ans);

        }

        printf("%lld\n",ans);

    }

}sam;

int main()

{

    int n;

    scanf("%d",&n);

    for(int i=;i<=n;i++)

    {

        scanf("%s",s+);

        sam.build();

    }

    sam.cal();

    return ;

}

/********************

1

101

********************/

版本1

//#pragma comment(linker, "/stack:200000000")

//#pragma GCC optimize("Ofast,no-stack-protector")

//#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")

//#pragma GCC optimize("unroll-loops")

#include<bits/stdc++.h>

#define fi first

#define se second

#define mp make_pair

#define pb push_back

#define pi acos(-1.0)

#define ll long long

#define mod (1000000007)

#define C 0.5772156649

#define ls l,m,rt<<1

#define rs m+1,r,rt<<1|1

#define pil pair<int,ll>

#define pii pair<int,int>

#define ull unsigned long long

#define base 1000000000000000000

#define fio ios::sync_with_stdio(false);cin.tie(0)

using namespace std;

const double g=10.0,eps=1e-;

const int N=+,maxn=+,inf=0x3f3f3f3f,INF=0x3f3f3f3f3f3f3f3f;

inline void sub(ll &x){if(x>=mod)x-=mod;}

char s[N];

int c[N],a[N];

ll val[N],sum[N],sz[N];

struct SAM{

    int last,cnt;

    int ch[N<<][],fa[N<<],l[N<<];

    int mx[N<<],tmp[N<<];

    void ins(int c)

    {

        int p=last,np=++cnt;last=np;l[np]=l[p]+;

        for(;p&&!ch[p][c];p=fa[p])ch[p][c]=np;

        if(!p)fa[np]=;

        else

        {

            int q=ch[p][c];

            if(l[p]+==l[q])fa[np]=q;

            else

            {

                int nq=++cnt;l[nq]=l[p]+;

                memcpy(ch[nq],ch[q],sizeof(ch[q]));

                fa[nq]=fa[q];fa[q]=fa[np]=nq;

                for(;ch[p][c]==q;p=fa[p])ch[p][c]=nq;

            }

        }

        sum[np]=1ll*c;

    }

    SAM()

    {

        cnt=;

    }

    void build(){

        int len=strlen(s+);

        last=;

        for(int i=;i<=len;i++)ins(s[i]-'');

//        topo();

//        for(int i=1;i<=cnt;i++)mx[i]=l[i],tmp[i]=0;

    }

    void topo(){

        for(int i=;i<=cnt;i++)c[l[i]]++;

        for(int i=;i<=cnt;i++)c[i]+=c[i-];

        for(int i=;i<=cnt;i++)a[c[l[i]]--]=i;

    }

    void debug()

    {

        puts("");

        for(int i=;i<=cnt;i++)

        {

            printf("%d ",i);

            for(int j=;j<;j++)

                printf("%d ",ch[i][j]);

            printf("%d %d %lld %lld\n",fa[i],l[i],val[i],sum[i]);

        }

    }

    void faupdate()

    {

        for(int i=cnt;i;i--)

            tmp[fa[a[i]]]=max(tmp[fa[a[i]]],tmp[a[i]]);

    }

    void maupdate()

    {

        for(int i=;i<=cnt;i++)

            mx[i]=min(mx[i],tmp[i]),tmp[i]=;

    }

    void match()

    {

        int len=strlen(s+),now=,res=;

        for(int i=;i<=len;i++)

        {

            if(ch[now][s[i]-'a'])

            {

                now=ch[now][s[i]-'a'];

                res++;

            }

            else

            {

                while(now&&!ch[now][s[i]-'a'])now=fa[now];

                if(!now)now=,res=;

                else

                {

                    res=l[now]+;

                    now=ch[now][s[i]-'a'];

                }

            }

            tmp[now]=max(tmp[now],res);

        }

        faupdate();

        maupdate();

    }

    int maxlen()

    {

        int ans=;

        for(int i=;i<=cnt;i++)ans=max(ans,mx[i]);

        return ans;

    }

    void cal()

    {

        topo();

        sz[]=;

        for(int i=;i<=cnt;i++)

        {

            int p=a[i];

//            val[p] =  val[p] + 1ll*sum[p] * (l[p]-l[fa[p]])%mod;

//            sub(val[p]);

            for(int j=;j<;j++)

            {

                if(ch[p][j])

                {

                    val[ch[p][j]]=val[ch[p][j]]+val[p]*10ll%mod+1ll*j*sz[p]%mod;

                    sub(val[ch[p][j]]);

                    sz[ch[p][j]]+=sz[p];

                }

            }

         //   cout << val[p] << endl;

        }

       // debug();

        ll ans=;

        for(int i=cnt;i;i--)

        {

            ans=ans+val[i];

            sub(ans);

        }

        printf("%lld\n",ans);

    }

}sam;

int main()

{

    int n;

    scanf("%d",&n);

    for(int i=;i<=n;i++)

    {

        scanf("%s",s+);

        sam.build();

//        if(i!=n)sam.ins(10);

    }

    sam.cal();

    return ;

}

/********************

1

987654321123456789

********************/

版本2

hihocoder1457的更多相关文章

【后缀自动机】【拓扑排序】【动态规划】hihocoder1457 后缀自动机四·重复旋律7
解题方法提示小Hi:我们已经学习了后缀自动机,今天我们再来看这道有意思的题. 小Ho:好!这道题目让我们求的是若干的数字串所有不同子串的和. 小Hi:你能不能结合后缀自动机的性质来思考如何解决本题? ...

随机推荐

Oracle数据库查找持有锁的SQL语句，而不是请求锁的SQL语句（原创）
Oracle数据库查找持有锁的SQL语句,而不是请求锁的SQL语句查找活动的事务以及活动事务关联的会话信息 select s.sid 会话ID, s.serial# 会话序列号, s.usernam ...
C++学习笔记--名称空间
名称空间是为了更好的控制名称的作用域,以管理不同的类库,避免发生冲突. 1.创建名称空间如下,使用namespace关键字创建了pers和debts两个名称空间. #ifndef NAMESP_H_ ...
08 网络配置、shh服务、bash命令和元字符
作业一:完成作业未做完的集群架构作业二:临时配置网络(ip,网关,dns)+永久配置作业三:为集群内的机器设定主机名,利用/etc/hosts文件来解析自己的集群中所有的主机名,相应的,集群的配置应 ...
SQL Server WITH ROLLUP、WITH CUBE、GROUPING语句的应用
CUBE:CUBE 生成的结果集显示了所选列中值的所有组合的聚合. ROLLUP:ROLLUP 生成的结果集显示了所选列中值的某一层次结构的聚合. GROUPING:当行由 CUBE 或 ROLLUP ...
转：centos彻底删除文件夹、文件命令
转自:http://www.cnblogs.com/kluan/p/4458296.html centos彻底删除文件夹.文件命令(centos 新建.删除.移动.复制等命令: 1.新建文件夹 mkd ...
Spring入门学习推荐
该作者的 spring 博客很好,作为学习用简介 https://blog.csdn.net/column/details/15088.html?&page=1 https://blog.c ...
maven 项目配置到tomcat不能正常启动
最近使用IntelliJ IDEA搭建公司项目,该项目是maven项目,加载jar和编译的时候没有任何异常,但是部署到tomcat上之后,就会出现如下异常: org.apache.catalina.L ...
ubuntu中mysql版本升级到5.7
0 前言前几天图书馆说服务器(Ubuntu14.04)有安全漏洞,不按时修复会关停. 看了一下漏洞清单,主要是ssh和mysql的版本问题. 把mysql升级了一下,升到了5.7,升级之前还备份了数 ...
NIO 03
1. 客户端要主动去连接:channel.connect(new InetSocketAddress("localhost",8888)); //用channel.finishCo ...
一只代码小白git托管路上的二三事
[经验]一只代码小白git托管路上的二三事写在前面的话寒假的时候,娄老师给我们布置了代码托管的作业,并要求把托管地址发给学委.因假期的时候没有带电脑回家,所以只是在手机上草草注册了,也稀里糊涂就将 ...

hihocoder1457

hihocoder1457的更多相关文章

随机推荐

热门专题