hdu6021[BestCoder #93] MG loves string
这场BC实在是有趣啊,T2是个没有什么算法但是细节坑的贪心+分类讨论乱搞,T3反而码起来很顺.
然后出现了T2过的人没有T3多的现象(T2:20人,T3:30人),而且T2的AC率是惨烈的不到3%
(我T2读入写挂,思路想挂,推倒重码一次,交了7次才过QAQ)
最后狂码T3赶在结束之前调出来过了,3题最后都没有被叉也没有fst,捡了个rk9,人品好啊.
题意
给出26个字母的一个排列,对一个字符串进行一次变换表示将某个字母变成排列中对应的字母.
例如zabcdefg....wxy,1表示经过一次变换之后,所有的'a'会变成'z','b'会变成'a','c'会变成'b'...
一个长度为N的随机串经过若干次变换之后必然得到本身.那么我们可以求出得到本身的最小变换次数.
长度为N的串有\(26^N\)个,求这些串的最小变换次数之和.\(N<=10^9\)
分析
给出的置换中,对于某一个字母,它至少变换多少次之后得到自身是已知的.对于某一个随机串,它变换到
自己的最小步数只和出现过的字母种类有关,对出现过的字母变换到自身的最小步数求lcm即可.那么我们
对出现过的字母种类进行状压,用\(f[i][S]\)表示前i个字母中出现的字母集合为S的方案数.那么\(f[i-1][]\)
到\(f[i][]\)的转移是可以看作矩阵乘法形式的,于是我们可以在\(2^{26}logN\)的时间内解决暴力.
仔细分析.如果某个字母变换到自己的最小次数是3,那么它位于一个长度为3的循环中,这3种字母变换到自己
的最小步数是相同的,因此可以在状态压缩的时候视为一种.如果两个字母处在不同的循环但是长度相同,例如
"badc....",'a,'b','c','d'的最小步数均为2,也可以视为一种.那么长度为26的串里最多出现6种不同的
循环长度,我们就可以在\(O(2^6logN)\)内解决.也可以dfs枚举出所有可行的lcm,数量不会比2^6多.
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
using namespace std;
const int mod=1000000007;
int sz;
struct matrix{
int a[130][130];
matrix(){
memset(a,0,sizeof(a));
}
matrix(int x){
memset(a,0,sizeof(a));
for(int i=0;i<130;++i){
a[i][i]=x;
}
}
matrix operator *(const matrix &B)const{
matrix C;
for(int i=1;i<=sz;++i){
for(int j=1;j<=sz;++j){
for(int k=1;k<=sz;++k){
C.a[i][k]=(C.a[i][k]+a[i][j]*1ll*B.a[j][k])%mod;
}
}
}
return C;
}
};
matrix qpow(matrix A,int x){
matrix ans(1);
for(;x;x>>=1,A=A*A){
if(x&1)ans=ans*A;
}
return ans;
}
char buf[30];
int a[30];
int vis[30],T;
int w[30];
int b[30];int tot=0;
int lcm[1024];int totlcm=0;
int gcd(int a,int b){
return b==0?a:gcd(b,a%b);
}
int LCM(int a,int b){
return a/gcd(a,b)*b;
}
void dfs(int x,int tmp){
if(x==tot+1){
lcm[++totlcm]=tmp;
return;
}else{
dfs(x+1,tmp);
dfs(x+1,LCM(tmp,b[x]));
}
}
matrix A,ANS;
int main(){
int tests;scanf("%d",&tests);
while(tests--){
tot=0;
int n;scanf("%d",&n);
scanf("%s",buf+1);
for(int i=1;i<=26;++i){
a[i]=buf[i]-'a'+1;
}
++T;
for(int i=1;i<=26;++i){
if(vis[i]!=T){
int cnt=0;
for(int j=i;vis[j]!=T;j=a[j]){
cnt++;vis[j]=T;
}
b[++tot]=cnt;
w[i]=cnt;
for(int j=a[i];j!=i;j=a[j]){
w[j]=cnt;
}
}
}
sort(b+1,b+tot+1);
int sum=0;
for(int i=1;i<=tot;++i){
if(b[i]!=b[i-1]){
b[++sum]=b[i];
}
}
tot=sum;
totlcm=0;
dfs(1,1);
sort(lcm+1,lcm+totlcm+1);
sum=0;
for(int i=1;i<=totlcm;++i){
if(lcm[i]!=lcm[i-1])lcm[++sum]=lcm[i];
}
totlcm=sum;
A=matrix(0);
map<int,int> dict;
for(int i=1;i<=totlcm;++i)dict[lcm[i]]=i;
for(int i=1;i<=totlcm;++i){
for(int j=1;j<=26;++j){
A.a[i][dict[LCM(lcm[i],w[j])]]++;//printf("%d %d\n",i,dict[LCM(lcm[i],w[j])]);
}
}sz=totlcm;
ANS=qpow(A,n);
int ans=0;
for(int i=1;i<=totlcm;++i){
ans=(ans+lcm[i]*1ll*ANS.a[1][i])%mod;
}
printf("%d\n",ans);
// for(int i=1;i<=totlcm;++i)printf("%d ",lcm[i]);printf("\n");
// for(int i=1;i<=tot;++i)printf("%d ",b[i]);printf("\n");
}
return 0;
}
hdu6021[BestCoder #93] MG loves string的更多相关文章
- 【HDU 6021】 MG loves string (枚举+容斥原理)
MG loves string Accepts: 30 Submissions: 67 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: ...
- hdu 6021 MG loves string
MG loves string Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others ...
- MG loves string
MG loves string Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others ...
- hdu 6021 MG loves string (一道容斥原理神题)(转)
MG loves string Accepts: 30 Submissions: 67 Time Limit: 2000/1000 MS (Java/Others) Memory ...
- ●HDU 6021 MG loves string
题链: http://acm.hdu.edu.cn/showproblem.php?pid=6021 题解: 题意:对于一个长度为 N的由小写英文字母构成的随机字符串,当它进行一次变换,所有字符 i ...
- hdu 6020 MG loves apple 恶心模拟
题目链接:点击传送 MG loves apple Time Limit: 3000/1500 MS (Java/Others) Memory Limit: 262144/262144 K (Ja ...
- best corder MG loves gold
MG loves gold Accepts: 451 Submissions: 1382 Time Limit: 3000/1500 MS (Java/Others) Memory Limit ...
- (set)MG loves gold hdu6019
MG loves gold Time Limit: 3000/1500 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others) ...
- 【HDU 6020】 MG loves apple (乱搞?)
MG loves apple Accepts: 20 Submissions: 693 Time Limit: 3000/1500 MS (Java/Others) Memory Limit: ...
随机推荐
- WPF中Popup控件在Win7以及Win10等中的对齐点方式不一样的解决方案 - 简书
原文:WPF中Popup控件在Win7以及Win10等中的对齐点方式不一样的解决方案 - 简书 最近项目中使用弹出控件Popup,发现弹出框的对齐方式在不同的系统中存在不同(Popup在win10上是 ...
- 优步uber司机申请了为什么一直没有通过审核,帐号也显示未激活
优步uber现在是越来越火,申请注册成为优步uber司机的人数也日剧增多,申请了的车主都知道,申请后要等待审核,审核通过才可以激活帐号,快的运气好的,三五天不到一个星期就激活了,慢点的得大半个月,还有 ...
- python3工作环境部署+spyder3+jupyter notebook
1.python3安装 1)官网去下载python3.7版本,双击安装,只要注意勾选写到PATH就行,其它直接NEXT. 2)安装完成,CMD键入 python 回车,跳出python界面就是成功. ...
- power sequece
- moment.js使用方法总结
Moment.js是一个轻量级的JavaScript时间库,它方便了日常开发中对时间的操作,提高了开发效率.日常开发中,通常会对时间进行下面这几个操作:比如获取时间,设置时间,格式化时间,比较时间等等 ...
- Laya 自适应 不拉伸处理
Laya.init(640, Laya.Browser.width / 640 * 1028, WebGL); Laya.stage.scaleMode = "fixedwidth" ...
- TW实习日记:第五天
今天可以说是非常忙的一天了,要再项目中实现微信相关的功能:授权登录以及扫码登录,还有就是自建应用的发送消息.首先功能代码其实在经过了几天的学习之后并没有很难,但是最让我难受的是在项目中去加代码,首先s ...
- Visionpro学习网
重码网是一个在线机器视觉学习网站,推出了Halcon,Visionpro机器视觉学习视频教程,视频内容通俗易懂,没有编程基础的同学,照着视频练习,也同样可以学会. 学机器视觉,拿高薪,成就技术大拿.重 ...
- SAP(ABAP) ABAP内部外部数据转换常用function
文本相关CONVERSION_EXIT_CUNIT_OUTPUT 将内部单位转为单位文本CONVERSION_EXIT_ISOLA_OUTPUT 根据语言代码取文本CONVERSI ...
- 第二次c++作业(觉得渐渐入门系列)
其实说实话,我还是不敢很确定地说面向对象和面向过程这两种语言,我确实能分得开,但是我觉得倒是比以前好很多了.//(大概是谈了对象,知道了什么是面向对象编程) 1.从个人角度来说, a:面向过程就是-- ...