2017年10月30日提高组T1 数论

Description

---

聪明的0v0正在学习莫比乌斯反演。
她看到了这样的一道题：有n*m个人站成了一个n*m的方阵……
剩下的题面，聪明的0v0不记得了。但是，她通过自己高超的数论技巧，给出了一个转化后的模型：给出n和m，求

聪明的0v0当然知道怎么做了，但是她想考考你。

Input

---

一行三个正整数n，m，p。

Output

---

一行一个非负整数，设答案为x，输出x mod p。

Hint

---

30% n,m<=2000 p=998244353。
30% n*m<=10^9 n,m<=10^5 p为质数
20% n,m<=10^6 p为质数
20% n,m<=10^7 p为合数
对于所有数据，保证p<=10^9

Solution

---

不难想到当只有gcd(i,j)=1的条件时就是能在原点直接看到的点的数量，这里乘上$\min(\lfloor{\frac{n}{i}}\rfloor, \lfloor{\frac{n}{i}}\rfloor)$ 就是枚举的这条斜率覆盖了多少个点。由于每个点最多被覆盖了一次，答案就是n*m%p了

注意要开ll

Code

---

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <iostream>
#include <queue>
#include <vector>
#include <algorithm>
#define rep(i, st, ed) for (int i = st; i <= ed; i += 1)
#define drp(i, st, ed) for (int i = st; i >= ed; i -= 1)
#define erg(i, st) for (int i = ls[st]; i; i = e[i].next)
#define fill(x, t) memset(x, t, sizeof(x))
#define min(x, y) ((x)<(y)?(x):(y))
#define max(x, y) ((x)>(y)?(x):(y))
#define ld long double
#define db double
#define ll long long
#define INF 0x3f3f3f3f
#define N 1001
#define E 1001
#define L 1001
inline ll read() {
    ll x = 0, v = 1;
    char ch = getchar();
    for (; ch < '0' || ch > '9'; v *= (ch == '-')?(-1):(1), ch = getchar());
    for (; ch <= '9' && ch >= '0'; (x *= 10) += ch - '0', ch = getchar());
    return x * v;
}
int main(void) {
    ll n = read();
    ll m = read();
    ll p = read();
    std:: cout << (n % p * m % p) % p << std:: endl;
    return 0;
}