LG P2155 【SDOI2008】沙拉公主的困惑解题报告

最新推荐文章于 2024-06-21 17:14:10 发布

日居月诸Rijuyuezhu

最新推荐文章于 2024-06-21 17:14:10 发布

阅读量236

点赞数

分类专栏：算法竞赛题解文章标签：算法数论欧拉函数

本文链接：https://blog.csdn.net/qq_41996523/article/details/126668531

版权

算法竞赛同时被 2 个专栏收录

5 篇文章 0 订阅

订阅专栏

题解

2 篇文章 0 订阅

订阅专栏

LG P2155 【SDOI2008】沙拉公主的困惑解题报告

题目链接

根据题意，所求即为

$ans=\sum_{i=1}^{N!}[\gcd(i,M!)=1]$

由于 $M\le N$ ，容易知道 $M!\mid N!$ 。

结合 $g c d$ 的性质，若 $\gcd(i,M!)=c$ ，则 $\gcd(i+kM!,M!)=c$ 。所以把 $M!$ 看成一种循环节，有

$ans=\dfrac{N!}{M!}\sum_{i=1}^{M!}[i\perp M!]=\dfrac{N!}{M!}\varphi(M!)$

用欧拉函数的定义展开，有

$ans=\dfrac{N!}{M!}\cdot M!\prod_{p|M!}\dfrac{p-1}{p}=N!\dfrac{\prod_{p\le M}(p-1)}{\prod_{p\le M}p}$

这样，预处理三部分，即可 $O (N)$ 预处理， $O (1)$ 计算。注意模数的分类讨论。

/*
 * @Author: rijuyuezhu
 * @Date: 2022-09-02 17:00:03
 * @Description: https://www.luogu.com.cn/problem/P2155
 * @Tag: 数论，欧拉函数
*/
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
char In[1 << 20], *ss = In, *tt = In;
#define getchar() (ss == tt && (tt = (ss = In) + fread(In, 1, 1 << 20, stdin), tt == ss) ? EOF : *ss++)
ll read() {
   ll x = 0, f = 1; char ch = getchar();
   for(; ch < '0' || ch > '9'; ch = getchar()) if(ch == '-') f = -1;
   for(; ch >= '0' && ch <= '9'; ch = getchar()) x = x * 10 + int(ch - '0');
   return x * f;
}
const int MAXN = 1e7 + 5;
int P;
int pr[MAXN], tot, fac[MAXN], pd_pd1[MAXN], ivpd_p[MAXN];
bool np[MAXN];
int qpow(int a, int n) {
   int ret = 1;
   for(; n; n >>= 1, a = 1ll * a * a % P)
      if(n & 1) ret = 1ll * ret * a % P;
   return ret;
}
void init(int n) {
   np[1] = 1;
   for(int i = 2; i <= n; i++) {
      if(!np[i])
         pr[++tot] = i;
      for(int j = 1; j <= tot && pr[j] <= n / i; j++) {
         np[i * pr[j]] = 1;
         if(i % pr[j] == 0) break;
      }
   }
   fac[0] = 1;
   for(int i = 1; i <= n; i++)
      if(i == P) fac[i] = fac[i-1];
      else fac[i] = 1ll * fac[i-1] * i % P;
   pd_pd1[0] = ivpd_p[0] = 1;
   for(int i = 1; i <= tot; i++) {
      pd_pd1[i] = 1ll * pd_pd1[i-1] * (pr[i] - 1) % P;
      if(pr[i] == P)
         ivpd_p[i] = ivpd_p[i-1];
      else 
         ivpd_p[i] = 1ll * ivpd_p[i-1] * qpow(pr[i], P-2) % P;
   }
}
int main() {
   int t = read();
   P = read();
   init(1e7);
   for(int i = 1; i <= t; i++) {
      int n = read(), m = read();
      int id = upper_bound(pr + 1, pr + 1 + tot, m) - pr - 1; // id of the first num that <= m
      if(P <= m)
         printf("%lld\n", 1ll * fac[n] * pd_pd1[id] % P * ivpd_p[id] % P); 
      else if(P < n)
         printf("0\n");
      else if(P == n) {
         if(m == n) printf("%lld\n", 1ll * fac[n] * pd_pd1[id] % P * ivpd_p[id] % P);
         else printf("0\n");
      } else //P > n
         printf("%lld\n", 1ll * fac[n] * pd_pd1[id] % P * ivpd_p[id] % P);
   }
   return 0;
}