[BZOJ 5330][SDOI2018] 反回文串
怎么说呢,一道不可多得的反演题吧,具体解释之后再补
#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] 反回文串的更多相关文章
- BZOJ 5330 Luogu P4607 [SDOI2018]反回文串 (莫比乌斯反演、Pollard Rho算法)
题目链接 (BZOJ) https://www.lydsy.com/JudgeOnline/problem.php?id=5330 (Luogu) https://www.luogu.org/prob ...
- [BZOJ5330][SDOI2018]反回文串
luogu bzoj sol 枚举一个长度为\(n\)为回文串,它的所有循环位移都可以产生贡献. 但是这样算重了.重复的地方在于可能多个回文串循环同构,或者可能有的回文串经过小于\(n\)次循环位移后 ...
- [SDOI2018]反回文串
题意 问有多少个长度为\(N\)且字符集大小为\(K\)的字符串可以通过回文串旋转 (把第一个字符移到最后)若干次得到.\(K\le N≤10^{18}\) 做法 ARC64F的加强版 设\(h(d) ...
- 【SDOI2018】反回文串(【ARC064 F】Rotated Palindromes 加强版)
题意 给你一个正整数 \(n\),求有多少字符集为 \(1\) 到 \(k\) 之间整数的字符串,使得该字符串可以由一个长度为 \(n\) 的回文串循环移位得到. ARC原题 \(100\%\) 的数 ...
- BZOJ 3676: [Apio2014]回文串
3676: [Apio2014]回文串 Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 2013 Solved: 863[Submit][Status ...
- bzoj 3676: [Apio2014]回文串 回文自动机
3676: [Apio2014]回文串 Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 844 Solved: 331[Submit][Status] ...
- 字符串(马拉车算法,后缀数组,稀疏表):BZOJ 3676 [Apio2014]回文串
Description 考虑一个只包含小写拉丁字母的字符串s.我们定义s的一个子串t的“出 现值”为t在s中的出现次数乘以t的长度.请你求出s的所有回文子串中的最 大出现值. Input 输入只有一行 ...
- ●BZOJ 3676 [Apio2014]回文串
题链: http://www.lydsy.com/JudgeOnline/problem.php?id=3676 题解: 后缀数组,Manacher,二分 首先有一个结论:一个串的本质不同的回文串的个 ...
- 「SDOI 2018」反回文串
题目大意: 求字符集大小为$k$长度为$n$的经循环移位后为回文串的数量. 题解: 这题是D1里最神的吧 考虑一个长度为$n$回文串,将其循环移位后所有的串都是满足要求的串. 但是显然这样计算会算重. ...
随机推荐
- 关于PHP的一个坑爹问题(页面刷新)
最近在用PHP做一个服务端和一个客户端,在快要完工的时候,出现了一个重大问题---- 当在客户端手动输入IP和端口的时候,一按连接,OK,连接成功,嘻嘻,就在我自以为大功告成的时候,来了个晴天霹雳,一 ...
- 136. Single Number唯一的一个只出现了一次的数字
[抄题]: Given a non-empty array of integers, every element appears twice except for one. Find that sin ...
- SSL认证
SSL认证 单向认证 1.发一串消息个对方 2.对方用私钥加密后返回 3.本方用对方的公钥解密,验证消息是否正确, 如果消息相同,则本方认可对方 双向认证 本方认证对方 对方认证本方
- Centos里没有lsb_release
查看Centos操作系统版本,输入指令 lsb_release -a 报无此命令 解决办法,安装lsb_release 1.执行指令:yum install -y redhat-lsb 2.安装完毕后 ...
- 修改QPushButton北京颜色和字体背景
项目需要修改按钮背景的颜色 QPalette pal = startBtn.palette(); //startBtn是我已经定义好的QPushButton对象 pal.setColor(QPalet ...
- java日期和时间转换字符
日期和时间转换字符 字符 描述 例子 c 完整的日期和时间 Mon May 04 09:51:52 CDT 2009 F ISO 8601 格式日期 2004-02-09 D U.S. 格式日期 (月 ...
- npm link和react native的问题
问题说明: 需要自己开发一个ReactNative插件,这个插件在独立git仓库,那么怎么把这个插件安装到主项目的依赖里,并且方便对插件的修改调试 方案一: 把插件发布到npm仓库,每次主项目通过np ...
- BOOL运算符号(从C#入门经典第五版中摘录)
只总结自己觉得难的哈: (1) var1=!var2; //(非) (2) var1=var2&var3; //(与) (3)var1=var2|var3; //(或) (4 ...
- 跨域问题hbuilder
1.借助jquery-jsonp插件 $.jsonp({ url: url, data: { 'name': usd, 'passwd': pass }, callbackParameter: &qu ...
- CentOS7下源码包方式安装rabbitmq
1.先安装erlang http://www.cnblogs.com/justphp/p/6093880.html 2.下载rabbitmq rpm包: wget http://www.rabbitm ...