传送门

怎么说呢,一道不可多得的反演题吧,具体解释之后再补

 #include <bits/stdc++.h>
 using namespace std;
 #define rep(i,a,b) for(int i=a;i<=b;++i)
 typedef long long ll;
 ;
 ll mul(ll x,ll y,ll p) {
     x%=p; y%=p;
     return (x*y-(ll)((long double)x/p*y+0.5)*p+p)%p;
 }
 ll _pow(ll x,ll n,ll p) {
     ll ret=;
     ,x=mul(x,x,p)) ) ret=mul(ret,x,p);
     return ret;
 }
 ll tp[]={2LL,3LL,5LL,7LL,13LL,61LL};
 bool MR(ll n) {
     ) return false;
     rep(i,,) if(n==tp[i]) return true;
     rep(i,,) ) return false;
     rep(i,,) {
         ll tmp=n-;)) tmp>>=;
         ll s=_pow(tp[i],tmp,n);
         &&s!=&&tmp!=n-) tmp<<=,s=mul(s,s,n);
         &&!(tmp&)) return false;
     }
     return true;
 }
 ll PR(ll n,ll c) {
     ll i=,k=2LL,x,y; x=y=1LL+rand()%(n-);
     ) {
         x=(mul(x,x,n)+c)%n;
         ll d=__gcd((y-x+n)%n,n);
         &&d!=n) return d;
         if(x==y) return n;
         ;
     }
 }
 int op[maxn],len,cnt,T;
 ll n,P,K,ans,gt[maxn];
 inline void fct(ll n) {
     ) return;
     if(MR(n)){gt[++len]=n;return;}
     ll p=n;
     ;p==n;--c) p=PR(p,c);
     fct(p); fct(n/p);
 }
 ll fpow(ll x,ll n,ll p) {
     ll ret=;
     ,x=x*x%p)
         ) ret=ret*x%p;
     return ret;
 }
 ll g(ll n) {,P);}
 ll f(ll n) {?n%P:(n>>)%P;}
 inline void dfs(int dp,ll d,ll pro) {
     ) {
         )&&(d&)==) return;
         (ans+=1LL*g(n/d)*f(n/d)%P*pro%P)%=P;
         return;
     }
     dfs(dp+,d,pro); pro=1LL*pro*(+P-gt[dp]%P)%P;
     rep(i,,op[dp]) d*=gt[dp],dfs(dp+,d,pro);
 }
 int main() {
 #ifndef ONLINE_JUDGE
     freopen("25.in","r",stdin);
 #endif
     scanf();
     while(T--) {
         scanf("%lld%lld%lld",&n,&K,&P);K%=P;
         len=cnt=;++cnt;
         memset(gt,,,sizeof(op));
         fct(n);
         sort(gt+,gt++len);
         rep(i,,len) {
             ;
             ++op[cnt];
         }
         ans=;dfs(,1LL,1LL);printf("%lld\n",ans);
     }
     ;
 }

[BZOJ 5330][SDOI2018] 反回文串的更多相关文章

  1. BZOJ 5330 Luogu P4607 [SDOI2018]反回文串 (莫比乌斯反演、Pollard Rho算法)

    题目链接 (BZOJ) https://www.lydsy.com/JudgeOnline/problem.php?id=5330 (Luogu) https://www.luogu.org/prob ...

  2. [BZOJ5330][SDOI2018]反回文串

    luogu bzoj sol 枚举一个长度为\(n\)为回文串,它的所有循环位移都可以产生贡献. 但是这样算重了.重复的地方在于可能多个回文串循环同构,或者可能有的回文串经过小于\(n\)次循环位移后 ...

  3. [SDOI2018]反回文串

    题意 问有多少个长度为\(N\)且字符集大小为\(K\)的字符串可以通过回文串旋转 (把第一个字符移到最后)若干次得到.\(K\le N≤10^{18}\) 做法 ARC64F的加强版 设\(h(d) ...

  4. 【SDOI2018】反回文串(【ARC064 F】Rotated Palindromes 加强版)

    题意 给你一个正整数 \(n\),求有多少字符集为 \(1\) 到 \(k\) 之间整数的字符串,使得该字符串可以由一个长度为 \(n\) 的回文串循环移位得到. ARC原题 \(100\%\) 的数 ...

  5. BZOJ 3676: [Apio2014]回文串

    3676: [Apio2014]回文串 Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 2013  Solved: 863[Submit][Status ...

  6. bzoj 3676: [Apio2014]回文串 回文自动机

    3676: [Apio2014]回文串 Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 844  Solved: 331[Submit][Status] ...

  7. 字符串(马拉车算法,后缀数组,稀疏表):BZOJ 3676 [Apio2014]回文串

    Description 考虑一个只包含小写拉丁字母的字符串s.我们定义s的一个子串t的“出 现值”为t在s中的出现次数乘以t的长度.请你求出s的所有回文子串中的最 大出现值. Input 输入只有一行 ...

  8. ●BZOJ 3676 [Apio2014]回文串

    题链: http://www.lydsy.com/JudgeOnline/problem.php?id=3676 题解: 后缀数组,Manacher,二分 首先有一个结论:一个串的本质不同的回文串的个 ...

  9. 「SDOI 2018」反回文串

    题目大意: 求字符集大小为$k$长度为$n$的经循环移位后为回文串的数量. 题解: 这题是D1里最神的吧 考虑一个长度为$n$回文串,将其循环移位后所有的串都是满足要求的串. 但是显然这样计算会算重. ...

随机推荐

  1. 关于PHP的一个坑爹问题(页面刷新)

    最近在用PHP做一个服务端和一个客户端,在快要完工的时候,出现了一个重大问题---- 当在客户端手动输入IP和端口的时候,一按连接,OK,连接成功,嘻嘻,就在我自以为大功告成的时候,来了个晴天霹雳,一 ...

  2. 136. Single Number唯一的一个只出现了一次的数字

    [抄题]: Given a non-empty array of integers, every element appears twice except for one. Find that sin ...

  3. SSL认证

    SSL认证 单向认证 1.发一串消息个对方 2.对方用私钥加密后返回 3.本方用对方的公钥解密,验证消息是否正确, 如果消息相同,则本方认可对方 双向认证 本方认证对方 对方认证本方

  4. Centos里没有lsb_release

    查看Centos操作系统版本,输入指令 lsb_release -a 报无此命令 解决办法,安装lsb_release 1.执行指令:yum install -y redhat-lsb 2.安装完毕后 ...

  5. 修改QPushButton北京颜色和字体背景

    项目需要修改按钮背景的颜色 QPalette pal = startBtn.palette(); //startBtn是我已经定义好的QPushButton对象 pal.setColor(QPalet ...

  6. java日期和时间转换字符

    日期和时间转换字符 字符 描述 例子 c 完整的日期和时间 Mon May 04 09:51:52 CDT 2009 F ISO 8601 格式日期 2004-02-09 D U.S. 格式日期 (月 ...

  7. npm link和react native的问题

    问题说明: 需要自己开发一个ReactNative插件,这个插件在独立git仓库,那么怎么把这个插件安装到主项目的依赖里,并且方便对插件的修改调试 方案一: 把插件发布到npm仓库,每次主项目通过np ...

  8. BOOL运算符号(从C#入门经典第五版中摘录)

    只总结自己觉得难的哈: (1) var1=!var2;    //(非) (2) var1=var2&var3;    //(与) (3)var1=var2|var3;    //(或) (4 ...

  9. 跨域问题hbuilder

    1.借助jquery-jsonp插件 $.jsonp({ url: url, data: { 'name': usd, 'passwd': pass }, callbackParameter: &qu ...

  10. CentOS7下源码包方式安装rabbitmq

    1.先安装erlang http://www.cnblogs.com/justphp/p/6093880.html 2.下载rabbitmq rpm包: wget http://www.rabbitm ...