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:你能不能结合后缀自动机的性质来思考如何解决本题? ...

随机推荐

Linux：Centos7升级内核（转）
更新前,内核版本为: uname -r 3.10.0-327.10.1.el7.x86_64 升级的方法: 1.导入key rpm --import https://www.elrepo.org/RP ...
编译安装MariaDB-10.0.21
一.源码编译安装gcc-5.1.0 1.下载gcc源码包 Download (HTTP): http://ftpmirror.gnu.org/gcc/gcc-5.2.0/gcc-5.2.0.tar.b ...
php int 与 datetime 转换
数据库日期类型是int类型的,该查询结果是datetime类型的 select from_unixtime( `dateline` ) from cdb_posts 如果原来类型是datetime类型 ...
java中全面的单例模式多种实现方式总结
单例模式的思想想整理一些 java 并发相关的知识,不知道从哪开始,想起了单例模式中要考虑的线程安全,就从单例模式开始吧. 以前写过单例模式,这里再重新汇总补充整理一下,单例模式的多种实现. 单例模 ...
Mybatis入门配置
MyBatis 本是apache的一个开源项目iBatis, 2010年这个项目由apache software foundation 迁移到了google code,并且改名为MyBatis .20 ...
PL/SQL编程-块编程
(1). 简单分类 |————过程(存储过程) | ...
vue-cli脚手架build目录中的webpack.base.conf.js配置文件
转载自:http://www.cnblogs.com/ye-hcj/p/7082620.html webpack.base.conf.js配置文件// 引入nodejs路径模块 var path = ...
asp.net时间显示
DateTime dt = DateTime.Now;// Label1.Text = dt.ToString();//2005-11-5 13:21:25// Label2.Text = ...
git 常用的撤销操作
git 开发时的使用流程,从同事那里学习的, 第一步,先建立一个新的分支,可以在GITLAB或通过命令git checkout -b newbranch 或通过IDE插件第二步,让项目切换到该分支, ...
CIFAR和SVHN在各CNN论文中的结果
CIFAR和SVHN结果加粗表示原论文中该网络的最优结果. 可以看出DenseNet-BC优于ResNeXt优于DenseNet优于WRN优于FractalNet优于ResNetv2优于ResNet ...

hihocoder1457

hihocoder1457的更多相关文章

随机推荐

热门专题