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的更多相关文章

  1. 【后缀自动机】【拓扑排序】【动态规划】hihocoder1457 后缀自动机四·重复旋律7

    解题方法提示 小Hi:我们已经学习了后缀自动机,今天我们再来看这道有意思的题. 小Ho:好!这道题目让我们求的是若干的数字串所有不同子串的和. 小Hi:你能不能结合后缀自动机的性质来思考如何解决本题? ...

随机推荐

  1. [golang note] 变量常量

    变量 • 变量声明 √ golang变量声明的关键字为var. √ golang变量声明时类型信息放置在变量名之后. ▶ 单个变量声明 ▪ 语法如下 var name type ▪ 示例如下 var ...

  2. docker swarm overlay stack 服务部署记录

    项目xxx(后端),xxx-ui前端(前后端分离的项目) 依赖mysql,elasticsearch.分别制作了四个镜像来做这件事.希望可以制作跨主机的部署,使用了swarm,以下是学习记录. 参考 ...

  3. iOS 绘制一个表盘时钟，秒针效果可以“扫秒/游走”

    最近自己 也尝试写了一个表盘时钟,初衷源于等车时候一个老奶奶问时间,我打开手机,时间数字对我来说相对敏感,但是老奶奶是看不清的,我想识别 还是看表盘 老远 看时针分针角度就可以识别当前时间. 于是我想 ...

  4. 对 META标签 的一点点了解

    ①META标签是啥 META标签,是HTML语言head区的一个辅助性标签.在几乎所有的page里,我们都可以看 到类似下面这段html代码: ---------------------------- ...

  5. Hive环境安装

    说明: (Hbase依赖于Hadoop,同时需要把元数据存放在mysql中),mysql自行安装 Hadoop2.0安装参考我的博客: https://www.cnblogs.com/654wangz ...

  6. linux 常用命令总结(三)

    1. setup       // 进入相应配置界面,按空格键选择相关功能 2. ll       // 列出当前目录下详细内容 :等价与ls -all 3. clear        // 清理当前 ...

  7. Qt 编码问题QTextCodec

    在学习计算机语言的时候, 关于字体编码问题, 一直是大家开始学习新语言比较头痛的问题, 在这边总结一下关于Qt图形框架开发的编码问题. 一般在Window开发环境里,是GBK编码,在Linux开发环境 ...

  8. pyinstaller 打包生成的exe文件，在其他电脑上报错

    解决方法: 1.第一种情况,在打包的时候不要加参数-w,看一下执行exe文件后出现的报错再看下一步的行动 2.应该是需要装一个VC 2015 x64(下载地址:https://www.microsof ...

  9. elasticsearch 5.x集群安装

    1. 下载 wget https://artifacts.elastic.co/downloads/elasticsearch/elasticsearch-5.5.1.tar.gz 2. 解压 为便于 ...

  10. https nginx配置

    cd /saas/conf/nginx/ mkdir key cd key 创建key: openssl req -nodes -newkey rsa:2048 -keyout server.key ...