题意

\(n\) 局石头剪刀布,设每局的贡献为赢的次数与输的次数之 \(\gcd\) ,求期望贡献乘以 \(3^{2n}\) ,定义若 \(xy=0\) 则,\(\gcd(x,y)=x+y\)

思路

不难得出

\[ans=3^n\sum_{i=0}^n\sum_{j=0}^{n-i}{n\choose i}{n-i\choose j}\gcd(i,j)
\]

对于正整数 \(n\) ,有如下表达式

\[n=\sum_{d|n}\varphi(d)
\]

那么

\[ans=3^n\sum_{i=1}^n\sum_{j=1}^{n-i}{n\choose i}{n-i\choose j}\sum_{d|\gcd(i,j)}\varphi(d)+3^n\cdot2\sum_{i=1}^n{n\choose i}i
\]

注意欧拉函数是不取到零的,若要严谨的写,零应该特判。

在外层枚举 \(d\) ,并只枚举为 \(d\) 倍数的 \(i,j\) ,

\[ans=3^n\sum_{d=1}^n\varphi(d)\sum_{i=1}^{\lfloor{n\over d}\rfloor}\sum_{j=1}^{\lfloor{n\over d}\rfloor-i}{n\choose id}{n-id\choose jd}+3^n\cdot2\sum_{i=1}^n{n\choose i}i
\]

另外,利用表达式

\[{n\choose m}={n\over m}{n-1\choose m-1}
\]

还可以化简后式(虽然对复杂度没有影响。。)

\[ans=3^nn!\sum_{d=1}^n\varphi(d)\sum_{i=1}^{\lfloor{n\over d}\rfloor}\sum_{j=1}^{{\lfloor{n\over d}\rfloor}-i}{1\over (id)!}{1\over (jd)!}\cdot{1\over(n-id-jd)!}+6^nn
\]

变成了卷积的形式,后面于 \(i+j\) 有关的项就当作最后乘出多项式的常数。

复杂度的证明如下

\[\begin{array}{}
T(n)&=\displaystyle\sum_{i=1}^n {n\over i}\log{n\over i}\\
&\geq \displaystyle\sum_{i=1}^n {n\over i}\log{n}\\
&= \log n\displaystyle\sum_{i=1}^n {n\over i}\\
&\geq n\log n\ln n
\end{array}
\]

那么此题可解。

代码

#include<bits/stdc++.h>
#define FOR(i,x,y) for(int i=(x),i##END=(y);i<=i##END;++i)
#define DOR(i,x,y) for(int i=(x),i##END=(y);i>=i##END;--i)
using namespace std;
template<typename T,typename _T>inline bool chk_min(T &x,const _T y){return y<x?x=y,1:0;}
template<typename T,typename _T>inline bool chk_max(T &x,const _T y){return x<y?x=y,1:0;}
typedef long long ll;
const double PI=acos(-1.0);
const int N=1<<17|5;
namespace _Maths
{
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
void exgcd(ll a,ll b,ll &x,ll &y)
{
if(!b){x=1,y=0;return;}
exgcd(b,a%b,y,x),y-=a/b*x;
}
ll Pow(ll a,ll p,ll P)
{
ll res=1;
for(;p>0;p>>=1,(a*=a)%=P)if(p&1)(res*=a)%=P;
return res;
}
ll inv(ll a,ll P){ll x,y;exgcd(a,P,x,y);return (x%P+P)%P;}
};
using namespace _Maths;
int P;
struct Complex
{
double x,y;
Complex operator +(const Complex &_)const{return (Complex){x+_.x,y+_.y};}
Complex operator -(const Complex &_)const{return (Complex){x-_.x,y-_.y};}
Complex operator *(const Complex &_)const{return (Complex){x*_.x-y*_.y,x*_.y+y*_.x};}
Complex operator /(const int &_)const{return (Complex){x/_,y/_};}
};
namespace _Polynomial
{
const int K=(1<<15)-1,L=15;
Complex A[N<<1],B[N<<1],C[N<<1],D[N<<1];
Complex w[N<<1];int r[N<<1];
void DFT(Complex *a,int op,int n)
{
FOR(i,0,n-1)if(i<r[i])swap(a[i],a[r[i]]);
for(int i=2;i<=n;i<<=1)
for(int j=0;j<n;j+=i)
for(int k=0;k<i/2;k++)
{
Complex u=a[j+k],t=w[op==1?n/i*k:(n-n/i*k)&(n-1)]*a[j+k+i/2];
a[j+k]=u+t,a[j+k+i/2]=u-t;
}
if(op==-1)FOR(i,0,n-1)a[i]=a[i]/n;
}
void multiply(const int *a,const int *b,int *c,int n1,int n2)
{
int n=1;
while(n<n1+n2-1)n<<=1;
FOR(i,0,n1-1)A[i]=(Complex){a[i]&K,a[i]>>L};
FOR(i,0,n2-1)B[i]=(Complex){b[i]&K,b[i]>>L};
FOR(i,n1,n-1)A[i]=(Complex){0,0};
FOR(i,n2,n-1)B[i]=(Complex){0,0};
FOR(i,0,n-1)r[i]=(r[i>>1]>>1)|((i&1)*(n>>1));
FOR(i,0,n-1)w[i]=(Complex){cos(2*PI*i/n),sin(2*PI*i/n)};
DFT(A,1,n),DFT(B,1,n);
FOR(i,0,n-1)
{
int j=(n-i)&(n-1);
C[i]=(Complex){0.5*(A[i].x+A[j].x),0.5*(A[i].y-A[j].y)}*B[i];
D[i]=(Complex){0.5*(A[i].y+A[j].y),0.5*(A[j].x-A[i].x)}*B[i];
}
DFT(C,-1,n),DFT(D,-1,n);
FOR(i,0,n1+n2-2)
{
ll s=C[i].x+0.5,t=C[i].y+0.5,u=D[i].x+0.5,v=D[i].y+0.5;
c[i]=(s+((t+u)%P<<L)%P+(v%P<<L<<L)%P)%P;
}
}
};
int A[N],B[N],C[N<<1];
int phi[N],fac[N],ifac[N];
int n; void init()
{
fac[0]=fac[1]=1;FOR(i,2,n)fac[i]=(ll)fac[i-1]*i%P;
ifac[0]=ifac[1]=1;FOR(i,2,n)ifac[i]=(ll)(P-P/i)*ifac[P%i]%P;
FOR(i,2,n)ifac[i]=(ll)ifac[i-1]*ifac[i]%P;
FOR(i,1,n)phi[i]=i;
FOR(i,2,n)if(phi[i]==i)
for(int j=i;j<=n;j+=i)
phi[j]=phi[j]/i*(i-1);
} int main()
{
int T;
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&n,&P);
init();
ll ans=0;
FOR(d,1,n)
{
ll sum=0;
FOR(i,1,n/d)A[(i)-1]=ifac[i*d];
FOR(i,1,n/d)B[(i)-1]=ifac[i*d];
_Polynomial::multiply(A,B,C,n/d,n/d);
FOR(i,1,n/d)(sum+=(ll)C[(i)-2]*ifac[n-i*d]%P)%=P;
(ans+=sum*phi[d]%P)%=P;
}
ans=(ans*Pow(3,n,P)%P*fac[n]%P+Pow(6,n,P)*n%P)%P;
printf("%lld\n",(ans+P)%P);
}
return 0;
}

HDU 6088 Rikka with Rock-paper-scissors(NTT+欧拉函数)的更多相关文章

  1. hdu 1286:找新朋友(数论,欧拉函数)

    找新朋友 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Submis ...

  2. HDU 2824.The Euler function-筛选法求欧拉函数

    欧拉函数: φ(n)=n*(1-1/p1)(1-1/p2)....(1-1/pk),其中p1.p2…pk为n的所有素因子.比如:φ(12)=12*(1-1/2)(1-1/3)=4.可以用类似求素数的筛 ...

  3. (hdu step 7.2.1)The Euler function(欧拉函数模板题——求phi[a]到phi[b]的和)

    题目: The Euler function Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Othe ...

  4. HDU 6088 - Rikka with Rock-paper-scissors | 2017 Multi-University Training Contest 5

    思路和任意模数FFT模板都来自 这里 看了一晚上那篇<再探快速傅里叶变换>还是懵得不行,可能水平还没到- - 只能先存个模板了,这题单模数NTT跑了5.9s,没敢写三模数NTT,可能姿势太 ...

  5. 2018 ACM-ICPC 中国大学生程序设计竞赛线上赛 H题 Rock Paper Scissors Lizard Spock.(FFT字符串匹配)

    2018 ACM-ICPC 中国大学生程序设计竞赛线上赛:https://www.jisuanke.com/contest/1227 题目链接:https://nanti.jisuanke.com/t ...

  6. HDU 2824 简单欧拉函数

    1.HDU 2824   The Euler function 2.链接:http://acm.hdu.edu.cn/showproblem.php?pid=2824 3.总结:欧拉函数 题意:求(a ...

  7. HDU 1695 GCD (欧拉函数+容斥原理)

    GCD Time Limit: 6000/3000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submiss ...

  8. HDU 5430 Reflect(欧拉函数)

    题目: http://acm.hdu.edu.cn/showproblem.php?pid=5430 从镜面材质的圆上一点发出一道光线反射NNN次后首次回到起点. 问本质不同的发射的方案数. 输入描述 ...

  9. HDU 4483 Lattice triangle(欧拉函数)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4483 题意:给出一个(n+1)*(n+1)的格子.在这个格子中存在多少个三角形? 思路:反着想,所有情 ...

随机推荐

  1. thinkphp5中使用phpmailer实现发送邮件功能(转载)

    一.开启SMTP服务(使用php发送邮件需要用到SMTP服务,这里以163邮箱的SMTP服务为例). 1.登录163邮箱,在首页上找到“设置”. 2.选择开启的服务,一般都全选,POP3/SMTP/I ...

  2. openssl内核升级

    由于工作需要,防止安全漏洞需要对openssl升级现在整理出centos6.8和ubuntu14.4升级 centos升级openssl 1.首先去OpenSSL的网站 https://www.ope ...

  3. Linux手动添加swap分区

    转自:https://blog.csdn.net/whatday/article/details/51024571 为什么需要swap 根据Redhat公司的建议,Linux系统swap分区最适合的大 ...

  4. Broadcast

    静态注册广播接收器 1. 活动中创建内部类继承BroadcastReceiver实现 onReceive函数 2. new 一个内部类的对象 3. registerReceiver注册内部类 4. 在 ...

  5. 洛谷P3376 【模板】网络最大流

    题目描述 如题,给出一个网络图,以及其源点和汇点,求出其网络最大流. 输入输出格式 输入格式: 第一行包含四个正整数N.M.S.T,分别表示点的个数.有向边的个数.源点序号.汇点序号. 接下来M行每行 ...

  6. 推举算法 AdaBoost 哥德尔奖 Godel Prize

    推举算法 AdaBoost  2003年理论计算机科学界最高奖 哥德尔奖 Godel Prize

  7. Fiddler (进阶)内置命令与断点

    Fiddler 内置命令与断点 命令 对应请求项 介绍 示例 ? All 问号后边跟一个字符串,可以匹配出包含这个字符串的请求 ?google > Body 大于号后面跟一个数字,可以匹配出请求 ...

  8. libvirt虚拟库

    转载自:https://www.ibm.com/developerworks/cn/linux/l-libvirt/index.html Libvirt 虚拟化库剖析   讲到向外扩展计算(比如云计算 ...

  9. .net ML机器学习中遇见错误记录

    避免入坑: 1.错误提示 numClasses must be at least   2 大概是训练模型的数据分类必须是两种,如下错误: 正确数据集如下:

  10. 012-mac下shell,zsh,oh-my-zsh,以及插件

    1.查看当前shell echo $SHELL 2.查看安装的shell cat /etc/shells 查看可知 /bin/bash /bin/csh /bin/ksh /bin/sh /bin/t ...