【THUSC 2017】如果奇迹有颜色【polya引理】【矩阵】【计数dp】【BM打表+线性递推】

题意:长度为 nnn 的环染 mmm 种颜色,要求任意相邻 mmm 个元素不能包含全部的颜色。求方案数 模 109+710^9+7109+7,循环同构。

n≤109,m≤7n\leq 10^9,m\leq7n109,m7

为啥我现在天天都在打表啊

先上 polya,对于移动 iii 位的置换不动点个数为 gcd⁡(i,n)\gcd(i,n)gcd(i,n)。而因为这个也是循环的,所以可以把每相邻的 n/gcd⁡(i,n)n/\gcd(i,n)n/gcd(i,n) 个看成一个环。设 f(n)f(n)f(n) 表示循环不同构的方案,显然答案就是

∑d∣nf(d)φ(nd)\sum_{d\mid n}f(d)\varphi(\frac nd)dnf(d)φ(dn)

这样就把循环同构搞掉了。

然后有一个暴力的dp:f(i,S)f(i,S)f(i,S) 表示填了 iii 个数,最后 mmm 个的颜色状态是 SSS,枚举开头 mmm 个转移到最后再看合不合法。

发现这个 dp 是个矩阵的形式。有个很强的结论:存在 nnn 阶矩阵递推等价于存在 nnn 阶线性递推。这个矩阵递推可以有很多含义,包括任意元素、所有元素的和,甚至矩阵本身。

也就是说,尽管这个限制很多,我们也有(xia)理(j)由(b)相(cai)信(xiang)这个 fff 一定存在不超过 mmm^mmm 次的线性递推式。然后你就可以打表BM线性递推一条龙服务了。

但是这个暴力 dp 实在太慢了,你光枚举开始状态都有 mmm^mmm ,后面再怎么都至少有个 mmm^mmm,写得不好可能还要多个 mmm^mmm,可能不能在可预见的时间内跑出来。考虑做点优化。

上面已经看出来了,复杂度瓶颈其实在枚举开头,考虑这部分怎么优化。注意到如果开头枚举的 mmm 个数中有相同的,那么结尾的时候检查就不用到这里了,也就是开头剩下的怎么填都不影响合法性。而之前的还可以重标号搞掉。

具体而言,我们枚举一个 www 表示开头有多少个连续的不同的数,即对于 i∈[1,w],ai=ii\in[1,w],a_i=ii[1,w],ai=i。然后后面一个数钦定与前面的相同,即 ai+1∈[1,w]a_{i+1}\in [1,w]ai+1[1,w]。之后就可以正常转移了。最后再乘上个 mw‾m^{\underline w}mw 加起来就是答案。

这样复杂度到了O(N⋅mm⋅poly⁡(m))O(N\cdot m^m\cdot \operatorname{poly}(m))O(Nmmpoly(m)),其中 NNN 为打前多少个数的表,考场上根据你的梦想和人品来决定,这里打了 100010001000 个。因为懒得进一步优化,就直接跑 55min⁡55\min55min 跑出来了,但我看其他人只跑了几十秒……

然后就是套板子的事情了。

暴力求 fff

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#define MAXN 1000005
using namespace std;
const int MOD=1e9+7,N=1000;
inline int add(const int& x,const int& y){return x+y>=MOD? x+y-MOD:x+y;}
typedef long long ll;
inline int qpow(int a,int p)
{int ans=1;while (p){if (p&1) ans=(ll)ans*a%MOD;a=(ll)a*a%MOD,p>>=1;}return ans;
}
int m,a[10],MAX,vis[MAXN];
inline bool check(int x)
{if (~vis[x]) return vis[x];for (int i=0;i<m;i++) a[i]=0;int t=x;for (int i=0;i<m;i++){a[x%m]=1;x/=m;}for (int i=0;i<m;i++) if (!a[i]) return (vis[t]=1);return (vis[t]=0);
}
int cur[MAXN],t[MAXN],ans[MAXN];
inline void solve()
{for (int w=1;w<m;w++){int st=0;for (int i=1;i<=w;i++) st=st*m+i;for (int k=1;k<=w;k++){for (int i=0;i<MAX;i++) cur[i]=0;cur[st*m+k]=1;int tmp=0;for (int tt=st*m+k,j=1;j<=w+1;tt=(tt*m+(j>w? k:j))%MAX,j++)if (!check(tt))goto end1;tmp=add(tmp,1);end1:;for (int i=m;i>=m-w+1;i--) tmp=(ll)tmp*i%MOD;ans[w+1]=add(ans[w+1],tmp);for (int T=1;T<=N-w-1;T++){for (int i=0;i<MAX;i++) t[i]=0;for (int i=0;i<MAX;i++)if (check(i))for (int v=0;v<m;v++){int j=(i*m+v)%MAX;if (check(j)) t[j]=add(t[j],cur[i]);}for (int i=0;i<MAX;i++) cur[i]=t[i];int tmp=0;for (int i=0;i<MAX;i++){for (int t=i,j=1;j<=w+1;t=(t*m+(j>w? k:j))%MAX,j++)if (!check(t))goto end;tmp=add(tmp,cur[i]);end:;}for (int i=m;i>=m-w+1;i--) tmp=(ll)tmp*i%MOD;ans[T+w+1]=add(ans[T+w+1],tmp);
//				if (T+w+1==m+1) ans[m]=add(ans[m],tmp);}}}
}
int main()
{memset(vis,-1,sizeof(vis));cin>>m;MAX=qpow(m,m);printf("%d\n",N);solve();for (int i=1;i<m;i++) ans[i]=qpow(m,i);for (int i=1;i<=N;i++) printf("%d ",ans[i]); return 0;
}

找递推式:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#define MAXN 2005
using namespace std;
inline int read()
{int ans=0;char c=getchar();while (!isdigit(c)) c=getchar();while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();return ans;
}
const int MOD=1e9+7;
typedef long long ll;
inline int add(const int& x,const int& y){return x+y>=MOD? x+y-MOD:x+y;}
inline int dec(const int& x,const int& y){return x<y? x-y+MOD:x-y;}
inline int qpow(int a,int p)
{int ans=1;while (p){if (p&1) ans=(ll)ans*a%MOD;a=(ll)a*a%MOD,p>>=1;}return ans;
}
int n,m,F[MAXN],R[MAXN],las[MAXN],p,delta,tmp[MAXN];
int main()
{n=read();for (int i=1;i<=n;i++) F[i]=read();for (int k=1;k<=n;k++){int res=0;for (int i=1;i<=R[0];i++) res=(res+(ll)F[k-i]*R[i])%MOD;if (res==F[k]) continue;if (!R[0]){R[0]=p=k,delta=F[k];continue;}memcpy(tmp,R,sizeof(tmp));int x=(ll)dec(F[k],res)*qpow(delta,MOD-2)%MOD;R[k-p]=add(R[k-p],x);for (int i=1;i<=las[0];i++) R[k-p+i]=dec(R[k-p+i],(ll)las[i]*x%MOD);R[0]=max(R[0],k+las[0]-p);memcpy(las,tmp,sizeof(las)),delta=dec(F[k],res),p=k;}printf("%d\n",R[0]);for (int i=1;i<=R[0];i++) printf("%d%c",F[i],",\n"[i==R[0]]);for (int i=1;i<=R[0];i++) printf("%d%c",R[i],",\n"[i==R[0]]);return 0;
}

交上去的程序:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
using namespace std;
const int len[]={0,0,1,7,18,47,134,413};
const int S[][500]=;
const int R[][500]=;
const int MOD=1e9+7;
typedef long long ll;
inline int add(const int& x,const int& y){return x+y>=MOD? x+y-MOD:x+y;}
inline int dec(const int& x,const int& y){return x<y? x-y+MOD:x-y;}
inline int qpow(int a,int p)
{int ans=1;while (p){if (p&1) ans=(ll)ans*a%MOD;a=(ll)a*a%MOD;p>>=1;}return ans;
}
int m,F[1000],ans[1000],cur[1000],t[1000];
inline void mul(int* F,int* G)
{for (int i=0;i<=(len[m]<<1);i++) t[i]=0;for (int i=0;i<=len[m];i++)for (int j=0;j<=len[m];j++)t[i+j]=(t[i+j]+(ll)F[i]*G[j])%MOD;for (int i=0;i<=(len[m]<<1);i++) F[i]=t[i];for (int i=(len[m]<<1);i>len[m];i--)if (F[i]){for (int j=1;j<=len[m];j++) F[i-j]=(F[i-j]+(ll)F[i]*R[m][j])%MOD;F[i]=0; }
}
inline void qpow(int* F,int p)
{ans[0]=1;for (int i=1;i<=len[m];i++) ans[i]=0;while (p){if (p&1) mul(ans,F);mul(F,F),p>>=1;}for (int i=0;i<=len[m];i++) F[i]=ans[i];
}
inline int calc(int n)
{for (int i=0;i<=len[m];i++) cur[i]=0;cur[1]=1;qpow(cur,n);int res=0;for (int i=1;i<=len[m];i++) res=(res+(ll)S[m][i]*ans[i])%MOD;return res;
}
inline int phi(int x)
{int ans=1;for (int i=2;i*i<=x;i++)if (x%i==0){ans*=i-1,x/=i;while (x%i==0) ans*=i,x/=i;}if (x>1) ans*=x-1;return ans;
}
int main()
{int n;cin>>n>>m;int ans=0;for (int i=1;i*i<=n;i++)if (n%i==0){ans=(ans+(ll)calc(i)*phi(n/i))%MOD;if (i*i<n) ans=(ans+(ll)calc(n/i)*phi(i))%MOD;}cout<<(ll)ans*qpow(n,MOD-2)%MOD;return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/314865.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

ASP.NET Core 3.0中支持AI的生物识别安全

本文共两个部分&#xff0c;这是第一部分&#xff0c;其中介绍了 ASP.NET Core 3 中旨在将授权逻辑与基本的用户角色相分离的基于策略的授权模型。此部分提供了此授权进程的基于生物识别信息&#xff08;如人脸识别或语音识别&#xff09;的具体示例。在此示例中&#xff0c;检…

Codeforces Round #588 (Div. 2) D. Marcin and Training Camp 思维

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 有nnn个人&#xff0c;每个人都有一个能力值bib_ibi​以及他会的技能aia_iai​&#xff0c;当他会第xxx个技能的时候&#xff0c;aia_iai​的第xxx位是111。定义当xxx不会某个技能但是yyy会的时候&#xff0…

架构杂谈《五》

保证最终一致性的模式在大规模、高并发服务化系统中&#xff0c;一个功能被拆分成多个具有功能单一的子功能&#xff0c;一个流程会有多个系统的多个单一功能的服务组合实现&#xff0c;如果使用两阶段提交协议和三阶段提交协议&#xff0c;确实能解决系统间的一致性问题。其实…

【ROI 2019 Day2】课桌【贪心】【决策单调性】【分治】

题意&#xff1a;有 mmm 个班&#xff0c;每个班有 2n2n2n 个人&#xff0c;他们的身高给定。有 kkk 种双人桌&#xff0c;每张桌子有两个属性值 Li,RiL_i,R_iLi​,Ri​,一个身高为 hhh 的人坐第 iii 种桌子的不舒适度为 hhh 到区间 [Li,Ri][L_i,R_i][Li​,Ri​] 的最小距离。你…

Codeforces Round #588 (Div. 2) E. Kamil and Making a Stream 数学 + 暴力

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一颗树&#xff0c;其中根是111&#xff0c;每个点有一个点权&#xff0c;求每个点到根的所有路径的gcdgcdgcd之和。 n≤1e5n\le1e5n≤1e5 思路&#xff1a; 一看到以为是个点分治&#xff0c;让后发现…

【NOIP模拟】彩色树【树形dp】【树链剖分性质】【复杂度分析】

题意&#xff1a;一棵初始时为空的树&#xff0c;依次加入 nnn 个叶结点&#xff0c;每次加入后询问 用若干不同颜色的路径将树边染色后 每个点到根经过的颜色数 的最大值 的最小值。 n≤106n\leq 10^6n≤106 首先发现这个路径没啥用&#xff0c;其实就是个剖分方案。 然后我…

在Linux的Windows子系统上(WSL)使用Docker(Ubuntu)

背景平时开发大部人都是在提供了高效GUI的window下工作&#xff0c;但是真正部署环境普遍都是在Linux中&#xff0c;所以为了让开发环境和部署环境统一&#xff0c;我们需要在windows模拟LInux环境&#xff0c;以前我们可能通过虚拟机的方式实现&#xff0c;不过自从微软拥抱开…

NOIP2020 赛前总结

没有一眼秒的题都必须对拍&#xff0c;这里对拍包括疯狂造小数据人工检查。所以不会的话放心打暴力吧&#xff0c;反正写出来了也要对拍。码农题最后码。预估可能会失误&#xff0c;不要认为写一半就丢了很可惜&#xff0c;其他更可做的题没碰更可惜。想起了就把文件夹复制一遍…

.net持续集成cake篇之cake介绍及简单示例

cake介绍Cake 是.net平台下的一款自动化构建工具,可以完成对.net项目的编译,打包,运行单元测试,集成测试甚至发布项目等等.如果有些特征Cake没有实现,我们还可以很容易地通过扩展Cake来实现我们想要的功能.Cake有以下特点1) 使用c#语言编写,可以在Cake脚本里使用C#语言来实现我…

Educational Codeforces Round 73 (Rated for Div. 2) Make The Fence Great Again dp + 结论

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先证明一个结论&#xff1a;一个数最多被加两次。 首先假设a[i]a[i−1]a[i]a[i-1]a[i]a[i−1]或a[i]a[i1]a[i]a[i1]a[i]a[i1]&#xff0c;那么此时可以将a[i]a[i]1a[i]a[i]1a[i]a[i]1&am…

.NET Core 3.0 发布小尺寸 self-contained 单体可执行程序

.NET Core 提供的发布应用程序选项 self-contained 是共享应用程序的好方法&#xff0c;因为应用程序的发布目录包含所有组件、运行时和框架。您只需要告诉使用者应用程序的入口 exe 文件&#xff0c;就可以使程序运行起来&#xff0c;而不必担心目标计算机上是否存在.NET Core…

NOIP 2020 游记

Day -1 校内考试出阴间题&#xff0c;体验极差。 晚上写了发压位高精&#xff0c;一发 AC&#xff0c;感觉很飘&#xff08;flag&#xff09; Day 1 8:00 进考场 8:10 分发压缩包&#xff0c;随便点开看看&#xff0c;看题目名字感觉很友好。然后就在一个叫 ball 的文件夹…

Educational Codeforces Round 73 (Rated for Div. 2) E. Game With String 思维博弈 好题(2500)

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 我们将每一段...拿出来看成若干段&#xff0c;将其分成以下四种情况&#xff1a; (1)len<b(1)len<b(1)len<b (2)b≤len<a(2)b\le len<a(2)b≤len<a (3)a≤len<2∗b(3…

「Sqlserver」数据分析师有理由爱Sqlserver之三-最值得使有低投入高产出的Sqlserver功能...

数据分析师群体&#xff0c;有别于一般的传统开发群体和数据库运维群体&#xff0c;对Sqlserver的功能需求上也各不相同&#xff0c;本篇以笔者的亲身经历&#xff0c;用一种有别于一般性教程的角度来讲解Sqlserver值得我们学习&#xff0c;投入产出比高的一些功能模块。当然&a…

【WC2019】数树【子集反演】【结论】【树形dp】【生成函数】【函数求导】【多项式全家桶】

题意&#xff1a;有两棵基于同一点集的树&#xff0c;点集大小为 nnn &#xff0c;两棵树中有 opopop 棵未确定&#xff0c;可以取所有 nn−2n^{n-2}nn−2 种可能。给每个点染上 [1,y][1,y][1,y] 中的一个颜色&#xff0c;要求若 uuu 到 vvv 在两棵树上的路径完全相同&#xff…

Educational Codeforces Round 73 (Rated for Div. 2) F. Choose a Square 线段树 + 二维转一维

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你nnn个点(xi,yi)(x_i,y_i)(xi​,yi​)&#xff0c;每个点有个价值cic_ici​&#xff0c;现在你可以框一个正方形&#xff0c;要求左下角和右上角的坐标(x,y)(x,y)(x,y)必须xyxyxy&#xff0c;也就是说必须…

.NET开发框架(八)-服务器集群之网络负载平衡(视频)

【视频教程在文章底部】&#xff0c;本文讲解Windows服务器集群的网络负载平衡NLB的作用&#xff0c;以及在.NET开发框架的架构设计中&#xff0c;如何应用NLB与ARR&#xff0c;使用它们各有什么优点。视频教程目录&#xff1a;1、讲解NLB概念与演示其作用 2、安装配置负载均衡…

【CF113D】Museum【概率期望】【高斯消元】

题意&#xff1a;一张 nnn 个点的无向连通图&#xff0c;两个人开始时分别在 a,ba,ba,b。每次在 uuu 时会以 ppp 的概率原地不动&#xff0c;1−p1-p1−p 的概率等概率随机选择到一个相邻的点&#xff0c;当两人在同一点时停止。分别求在每个点相遇的概率。 n≤22n\leq 22n≤22…

人生如戏,别太入戏

这里是Z哥的个人公众号每周五早8点 按时送达当然了&#xff0c;也会时不时加个餐&#xff5e;我的第「81」篇原创敬上大家好&#xff0c;我是Z哥。对&#xff0c;就是脑袋很大的那个。大到什么程度呢&#xff1f;我给新来的小伙伴们说说。我昨天还打算配副新眼镜来着&#xff0…

Codeforces Round #588 (Div. 2) F. Konrad and Company Evaluation 图论 + 建反图 好题

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一张nnn个点mmm条边的图&#xff0c;其中每个点iii初始编号为iii&#xff0c;边是有向的&#xff0c;方向为从编号大的指向编号小的。定义一个贡献为存在某三个点a,b,ca,b,ca,b,c有两条边为a−>b,b−…