题目链接

题意 : 给你一个串、要你将其划分成两个串、使得左边的串的本质不同回文子串的个数是右边串的两倍、对于每一个这样子的划分、其对答案的贡献就是左边串的长度、现在要你找出所有这样子的划分、并将贡献乘起来、答案 mod 1e9+7

分析 :

从左到右跑一边回文自动机、对于每个前缀

能够得出其有多少个本质不同的回文子串

本质不同的回文子串的个数实际上就是自动机节点数 - 2

那么跑一遍前缀之后我们能得到所有可作为左边部分串的本质不同回文子串的个数

因为是回文串、所以我们倒着跑一遍、就同样能得到作为右边部分串的本质不同回文子串的个数

最后暴力检查一遍对于每一个位置是否有符合题意的合理划分、如果有就将左边部分长度累乘起来

#include<bits/stdc++.h>
#define LL long long
#define ULL unsigned long long

#define scl(i) scanf("%lld", &i)
#define scll(i, j) scanf("%lld %lld", &i, &j)
#define sclll(i, j, k) scanf("%lld %lld %lld", &i, &j, &k)
#define scllll(i, j, k, l) scanf("%lld %lld %lld %lld", &i, &j, &k, &l)

#define scs(i) scanf("%s", i)
#define sci(i) scanf("%d", &i)
#define scd(i) scanf("%lf", &i)
#define scIl(i) scanf("%I64d", &i)
#define scii(i, j) scanf("%d %d", &i, &j)
#define scdd(i, j) scanf("%lf %lf", &i, &j)
#define scIll(i, j) scanf("%I64d %I64d", &i, &j)
#define sciii(i, j, k) scanf("%d %d %d", &i, &j, &k)
#define scddd(i, j, k) scanf("%lf %lf %lf", &i, &j, &k)
#define scIlll(i, j, k) scanf("%I64d %I64d %I64d", &i, &j, &k)
#define sciiii(i, j, k, l) scanf("%d %d %d %d", &i, &j, &k, &l)
#define scdddd(i, j, k, l) scanf("%lf %lf %lf %lf", &i, &j, &k, &l)
#define scIllll(i, j, k, l) scanf("%I64d %I64d %I64d %I64d", &i, &j, &k, &l)

#define lson l, m, rt<<1
#define rson m+1, r, rt<<1|1
#define lowbit(i) (i & (-i))
#define mem(i, j) memset(i, j, sizeof(i))

#define fir first
#define sec second
#define VI vector<int>
#define ins(i) insert(i)
#define pb(i) push_back(i)
#define pii pair<int, int>
#define VL vector<long long>
#define mk(i, j) make_pair(i, j)
#define all(i) i.begin(), i.end()
#define pll pair<long long, long long>

#define _TIME 0
#define _INPUT 0
#define _OUTPUT 0
clock_t START, END;
void __stTIME();
void __enTIME();
void __IOPUT();
using namespace std;

struct Palindromic_Tree {
    ;
     ;

    deque<int> text;
    int len[maxn];///len[i]表示节点i表示的回文串的长度
    int cnt[maxn];///第 i 节点表示的回文串出现的次数、注意最后调用 getCnt 函数完成计算
    int num[maxn];///表示以节点 i 代表的回文串端点为端点的不同回文串个数
    int nxt[maxn][Letter];///nxt静态指针和字典树类似,指向的串为当前串两端加上同一个字符构成
    int fail[maxn];///失配指针
    int odd, even;///奇偶原始节点
    int pre, suf;
    int totNode;///自动机所有节点的个数
    LL totNum;///串中所有回文串数量

    int newNode(int _len){
        int ret = totNode++;
        len[ret] = _len;
        cnt[ret] = num[ret] = ;
        fail[ret] = ;
        memset(nxt[ret], , sizeof(nxt[ret]));
        return ret;
    }

    inline void init(){
        text.clear();
        totNode = ;
        totNum = ;
        even = newNode();
        odd = newNode(-);
        fail[even] = fail[odd] = odd;
        pre = suf = even;
    }

    int encode(char c){
        return c - 'a';
    }

    template<typename FunT>
    int getFail(int ch, int cur, FunT getNext){
        for(int i = getNext(len[cur]);
            i <  || i >= (int)text.size() || text[i] != ch;
            cur = fail[cur], i = getNext(len[cur]));
        return cur;
    }

    template<typename FunT>
    int Insert(int ch, int last, FunT getNext){
        last = getFail(ch, last, getNext);
        if(!nxt[last][ch]){
            );
            fail[cur] = nxt[getFail(ch, fail[last], getNext)][ch];
            nxt[last][ch] = cur;
            num[cur] = num[fail[cur]] + ;
        }

        last = nxt[last][ch];
        totNum += (LL)num[last];
        cnt[last]++;
        return last;
    }

    int push_front(char c){
        text.push_front(encode(c));
        pre = Insert(encode(c), pre, [](; });
        if(len[pre] == (int)text.size()) suf = pre;
        return totNode;
    }

    int push_back(char c){
        text.push_back(encode(c));
        suf = Insert(encode(c), suf, [; });
        if(len[suf] == (int)text.size()) pre = suf;
        return totNode;
    }

    int push_front(char *s){
        ;
        ; s[i]; i++)
            ret = push_front(s[i]);
        return ret;
    }

    int push_back(char *s){
        ;
        ; s[i]; i++)
            ret = push_back(s[i]);
        return ret;
    }

    inline void getCnt(){
        ; i>=; i--)
            cnt[fail[i]] += cnt[i];
    }

}PAM1, PAM2;

;
;

char str[maxn];
LL L[maxn], R[maxn];

int main(void){__stTIME();__IOPUT();

    int T;
    sci(T);
    while(T--){
        mem(L, );
        mem(R, );
        PAM1.init();
        PAM2.init();
        scs(str);
        int len = strlen(str);
        ; i<len; i++){
            PAM1.push_back(str[i]);
            L[i] = PAM1.totNode - ;
        }
        ; i>=; i--){
            PAM2.push_back(str[i]);
            R[i] = PAM2.totNode - ;
        }

        LL ans = ;
        ; i<len; i++){
             * R[i+]){
                ) ans = 1LL;
                ans = (1LL * (i+) * ans) % mod;
            }
        }

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

__enTIME();;}

void __stTIME()
{
    #if _TIME
        START = clock();
    #endif
}

void __enTIME()
{
    #if _TIME
        END = clock();
        cerr<<"execute time = "<<(double)(END-START)/CLOCKS_PER_SEC<<endl;
    #endif
}

void __IOPUT()
{
    #if _INPUT
        freopen("in.txt", "r", stdin);
    #endif
    #if _OUTPUT
        freopen("out.txt", "w", stdout);
    #endif
}

最新文章

  1. 算法与数据结构(七) AOV网的拓扑排序
  2. 04. Web大前端时代之:HTML5+CSS3入门系列~HTML5 表单
  3. fluent批量处理&mdash;&mdash;模型参数的设置
  4. Arcgis10安装说明
  5. 作业八—Alpha阶段项目总结
  6. hdu 5146 Sequence
  7. hdu 1526(最大匹配)
  8. D&amp;F学数据结构系列——插入排序
  9. 给定一个函数rand()能产生0到n-1之间的等概率随机数,问如何产生0到m-1之间等概率的随机数?
  10. [翻译]AOP编程
  11. 模拟美萍加密狗--Rockey2虚拟狗(四)
  12. Oracle单表的复杂查询
  13. [置顶] ios 时间定时器 NSTimer应用demo
  14. TortoiseGit HTTPS方式保存密码最简单的方法
  15. hdu_5680_zxa and set(想法题)
  16. 关于foo的一个面试题
  17. Java程序设计第一次作业
  18. 【剑指offer】【python】面试题2~5
  19. 使用kbmmw中的随机数替换delphi 10.3 自带的随机数
  20. luogu 1850 换教室 概率+dp

热门文章

  1. selenium获取标签中的文本
  2. 【转】Entity Framework 6 Code First 实践系列(1):实体类配置-根据依赖配置关系和关联
  3. Date及DateFormat用法
  4. 怎样禁用浏览器的Cookie功能
  5. javaIO——PipedReader 和 PipedWriter 实现模拟即时聊天
  6. Linux常用命令(自用)
  7. SpringBoot-JPA入门
  8. 关于SpringMVC映射模型视图的几点小事
  9. Linux下mysql不区分大小写设置
  10. go语言入门(9)文本文件处理