如果两个数$a_{x}$和$a_{y}$,$\exists 0<i,a_{x}^{i}\equiv a_{y}(mod\ p^{k})$,就建一条$x$到$y$的有向边,再对这张图强连通分量缩点,记$s_{i}$表示第$i$个点的大小,$f_{i}$表示能到达$i$的点(初始)数量,则答案为$\sum_{i=1}^{scc}\frac{2^{s_{i}}-1}{2^{s_{i}+f_{i}}}$(对每一个点分别统计贡献)
问题就变为如何判定$\exists i,x^{i}\equiv y(mod\ p^{k})$,分为$p\mid x$和$p\nmid x$两类
对于$x\mid p$和$y\nmid p$,$x^{i}\equiv y(mod\ p^{k})$的必要条件为$x^{i}\equiv y(mod\ p)$,而$x^{i}\equiv 0\ne y(mod\ p)$(交换$x$和$y$也同理可以证明),因此两类之间没有关系
对于$x\mid p$和$y\mid p$,必然有$\forall k\le i,x^{i}\equiv 0(mod\ p^{k})$,暴力枚举$i$即可,时间复杂度$o(n\log p)$
对于$x\nmid p$和$y\nmid p$,如果求出$且ord(x)=\min_{i>0且x^{i}\equiv 1(mod\ p^{k})}i$,那么$x^{i}\equiv y(mod\ p^{k})$当且仅当$ord(y)\mid ord(x)$
用BSGS求$ord(x)$复杂度及常数较大,考虑$ord(x)\mid \varphi(p)$,枚举$\varphi(p)$的质因子并判断能否消除,复杂度为$o(n\log^{2}p)$

 1 #include<bits/stdc++.h>
2 using namespace std;
3 #define N 5005
4 #define mod 998244353
5 map<int,int>mat;
6 vector<int>v,v1[N],v2[N],v3[N];
7 int n,m,P,p,k,ans,a[N],ord[N],vis[N],dfn[N],sz[N],bl[N],dp[N];
8 int ksm(int n,int m,int p){
9 if (!m)return 1;
10 int s=ksm(n,m>>1,p);
11 s=1LL*s*s%p;
12 if (m&1)s=1LL*s*n%p;
13 return s;
14 }
15 int bsgs(int k){
16 int ans=1;
17 for(int i=0;i<v.size();i++)ans=ans*v[i];
18 for(int i=0;i<v.size();i++)
19 if (ksm(k,ans/v[i],P)==1)ans/=v[i];
20 return ans;
21 }
22 void add(int x,int y){
23 v1[x].push_back(y);
24 v2[y].push_back(x);
25 }
26 void dfs1(int k){
27 if (vis[k])return;
28 vis[k]=1;
29 for(int i=0;i<v1[k].size();i++)dfs1(v1[k][i]);
30 dfn[++dfn[0]]=k;
31 }
32 void dfs2(int k){
33 if (bl[k]){
34 if (bl[k]!=m)v3[m].push_back(bl[k]);
35 return;
36 }
37 sz[m]++;
38 bl[k]=m;
39 for(int i=0;i<v2[k].size();i++)dfs2(v2[k][i]);
40 }
41 void dfs3(int k,int x){
42 if (vis[k])return;
43 vis[k]=1;
44 dp[x]+=sz[k];
45 for(int i=0;i<v3[k].size();i++)dfs3(v3[k][i],x);
46 }
47 void init(){
48 p=P;
49 for(int i=3;i*i<=P;i++)
50 if (P%i==0){
51 p=i;
52 break;
53 }
54 for(int i=P;i>1;i/=p)k++;
55 int j=P/p*(p-1);
56 for(int i=2;i*i<=j;i++)
57 while (j%i==0){
58 j/=i;
59 v.push_back(i);
60 }
61 if (j>1)v.push_back(j);
62 }
63 void Kosaraju(){
64 memset(vis,0,sizeof(vis));
65 for(int i=1;i<=n;i++)dfs1(i);
66 for(int i=n;i;i--)
67 if (!bl[dfn[i]]){
68 m++;
69 dfs2(dfn[i]);
70 }
71 }
72 int main(){
73 scanf("%d%d",&n,&P);
74 init();
75 for(int i=1;i<=n;i++){
76 scanf("%d",&a[i]);
77 mat[a[i]]=i;
78 }
79 for(int i=1;i<=n;i++)
80 if (a[i]%p)ord[i]=bsgs(a[i]);
81 else{
82 for(int j=1,t=a[i];t>0;j++,t=1LL*t*a[i]%P)
83 if ((mat[t])&&(mat[t]!=i))add(i,mat[t]);
84 }
85 for(int i=1;i<=n;i++)
86 for(int j=1;j<=n;j++)
87 if ((i!=j)&&(ord[i])&&(ord[j])&&(ord[i]%ord[j]==0))add(i,j);
88 Kosaraju();
89 for(int i=1;i<=m;i++){
90 memset(vis,0,sizeof(vis));
91 dfs3(i,i);
92 }
93 for(int i=1;i<=n;i++)ans=(ans+(ksm(2,sz[i],mod)-1LL)*ksm(2,n-dp[i],mod))%mod;
94 printf("%d",ans);
95 return 0;
96 }

[loj3330]猜数游戏的更多相关文章

  1. 【BZOJ 1594】 [Usaco2008 Jan]猜数游戏 (二分+并查集)

    1594: [Usaco2008 Jan]猜数游戏 Description 为了提高自己低得可怜的智商,奶牛们设计了一个新的猜数游戏,来锻炼她们的逻辑推理能力. 游戏开始前,一头指定的奶牛会在牛棚后面 ...

  2. usaco 猜数游戏

    Description 为了提高智商,锻炼思维能力,奶牛设计了一个猜数游戏.游戏开始前,贝西会在牛棚后面摆上N个数字.所有数字排成一条直线,按次序从1到N编号.每个数字在1到10^9之间,没有两个数字 ...

  3. (一)Python之猜数游戏

    猜数游戏由简如深的编码学习过程: 3次机会: print('------------------我爱鱼C工作室------------------')count=0while count < 3 ...

  4. 猜数游戏-flag的运用

    package my;import java.util.Scanner;public class MyJava {        public static void main(String[] ar ...

  5. Java课程设计——猜数游戏(201521123111 陈伟泽)

    Java课程设计--猜数游戏(201521123111 陈伟泽) 1.团队课程设计博客链接 博客作业--猜数游戏 2.个人负责模块或任务说明 Answer:一些基础界面的构造,排行榜的构造,用文件录入 ...

  6. JAVA课程设计 猜数游戏 团队

    团队名称,成员介绍 名称: 猜数游戏 成员: 网络1514 201521123086 周颖强 网络1514 201521123087蒋勃超 项目git地址 git.oschina.net/jbc113 ...

  7. BZOJ 2222: [Cqoi2006]猜数游戏【神奇的做法,傻逼题,猜结论】

    2222: [Cqoi2006]猜数游戏 Time Limit: 20 Sec  Memory Limit: 259 MBSubmit: 604  Solved: 260[Submit][Status ...

  8. 从scratch到python——猜数游戏

    ` 之前讲解了从scratch到python,基于python turtle库的实现,讲解了用scratch和python turtle绘图的实现,以及让小猫动起来和当角色被单击的例子. 本节课讲继续 ...

  9. VBS猜数游戏

    VBS 猜数游戏 2018-11-09  21:19:11 by xutao msgbox "The Swami" ,,"Game" msgbox " ...

随机推荐

  1. 数据应用的变与不变,ShardingSphere 正在影响未来数字体验的建设理念

    近年来关于底层数据库的开源产品越来越多,它们也受到了许多资本的青睐. 伴随着移动互联网催生的数字化场景爆发,云计算.大数据等技术逐渐有了更加广阔的应用场景.在云计算和大数据经过十年的追赶式发展后,不只 ...

  2. Dapr-简介及环境搭建

    一.Dapr是什么? Dapr 是一个可移植的.事件驱动的运行时,它使任何开发人员能够轻松构建出弹性的.无状态和有状态的应用程序,并可运行在云平台或边缘计算中,它同时也支持多种编程语言和开发框架. 在 ...

  3. Java(5)输入和输出

    作者:季沐测试笔记 原文地址:https://www.cnblogs.com/testero/p/15201515.html 博客主页:https://www.cnblogs.com/testero ...

  4. mall笔记

    介绍 SpringBoot.SpringCloud.SpringCloudAlibaba.Nacos.Sentinel.Seata整合demo. 软件架构 JDK 1.8 Spring Boot 2. ...

  5. LeetCode:并查集

    并查集 这部分主要是学习了 labuladong 公众号中对于并查集的讲解,文章链接如下: Union-Find 并查集算法详解 Union-Find 算法怎么应用? 概述 并查集用于解决图论中「动态 ...

  6. Java:重载和重写

    Java:重载和重写 对 Java 中的 重载和重写 这个概念,做一个微不足道的小小小小结 重载 重载:编译时多态,同一个类中的同名的方法,参数列表不同,与返回值无关. 有以下几点: 方法名必须相同: ...

  7. MiniFly四轴飞行器之部分系统及电源分析

    最近硬件四轴很火,了解了很久,还是选择了MiniFly,主要还是资料多,后边可以有人讨论,不像很多就是建了个群,研究问题还是在论坛方便很多. 四轴终于拿到手,功能很强大,主要是还支持二次开发,可以研究 ...

  8. Netty:Reactor Pattern 与 Dubbo 底层传输中的 NettyServer

    首先,我们需要了解Reactor模式的三种线程模型: 1)单线程模型 Reactor 单线程模型,指的是所有的 IO 操作都在同一个 NIO 线程上面完成,NIO 线程的职责如下: 作为 NIO 服务 ...

  9. HTML+CSS基础(HTML篇)

    引言 在日常开发Android中,很多时候会遇到和WebView打交道,对CSS HTML JS不是很清楚的话是完不成一些功能的,本篇开始学习HTML,文章的主要内容是总结了慕课网中,HTML+CSS ...

  10. PyPi到底是什么?pypi有啥作用?PyPi和pip有何渊源?

    转载:https://blog.csdn.net/weixin_42139375/article/details/82711201 可能有很多刚入行不久的朋友们,每天都在用pip 命令install ...