2N——数论+快速幂

Posted z8875

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了2N——数论+快速幂相关的知识,希望对你有一定的参考价值。

题目

为了使得大家高兴,小Q特意出个自认为的简单题(easy)来满足大家,这道简单题是描述如下:
有一个数列A已知对于所有的(A[ i ])都是(1~n)的自然数,并且知道对于一些(A[ i])不能取哪些值,我们定义一个数列的积为该数列所有元素的乘积,要求你求出所有可能的数列的积的和 (mod 1000000007)的值,是不是很简单呢?呵呵!

Input

第一行三个整数(n,m,k)分别表示数列元素的取值范围,数列元素个数,以及已知的限制条数。
接下来(k)行,每行两个正整数(x,y)表示(A[ x])的值不能是(y)

Output

一行一个整数表示所有可能的数列的积的和对(1000000007)取模后的结果。如果一个合法的数列都没有,答案输出(0)

Sample Input

3 4 5
1 1
1 1
2 2
2 3
4 3

Sample Output

90

样例解释

(A[ 1])不能取(1)
(A[ 2])不能去(2、3)
(A[ 4])不能取(3)
所以可能的数列有以下(12)

数列
2 1 1 1 2
2 1 1 2 4
2 1 2 1 4
2 1 2 2 8
2 1 3 1 6
2 1 3 2 12
3 1 1 1 3
3 1 1 2 6
3 1 2 1 6
3 1 2 2 12
3 1 3 1 9
3 1 3 2 18

数据范围

  • (30%)的数据(n<=4,m<=10,k<=10)

  • 另有(20%)的数据(k=0)

  • (70%)的数据(n<=1000,m<=1000,k<=1000)

  • (100%)的数据 (n<=10^9,m<=10^9,k<=10^5,1<=y<=n,1<=x<=m)

题解

解题思路

考虑如果没有限制条件,根据加法原理和乘法原理可得
((1+2+3+……+n)^m=n*(n+1)^m)
单独计算加上限制后这位数的情况

代码

#include <cstdio>
#include <map>
#define int long long
using namespace std;
const int N = 1e5+5, M = 1e9+7;
int qpow(int a, int x) {
    a %= M;
    int ans = 1;
    while (x) {
        if (x & 1) ans *= a % M, ans %= M;
        x >>= 1;
        a *= a % M, a %= M;
    }
    return ans % M;
}
map<pair<int , int>, int > v;
map<int, int> s;
int n, m, k, a[N], tot, ans = 1, sum;
signed main() {
    scanf("%lld%lld%lld", &n, &m, &k);
    sum = n * (n + 1) / 2;
    for(int i = 1; i <= k; i++) {
        int x, y;
        scanf("%lld%lld", &x, &y);
        if (!s[x]) a[++tot] = x;
        if (v[make_pair(x, y)]) continue;
        v[make_pair(x, y)] = 1;
        s[x] += y;
    }
    for(int i = 1; i <= tot; i++)
        ans *= (sum - s[a[i]]) % M, ans %= M;
    ans = ans * qpow(sum, m - tot) % M;
    printf("%lld", ans);
    return 0 ;
}

以上是关于2N——数论+快速幂的主要内容,如果未能解决你的问题,请参考以下文章

HDU 2256 Problem of Precision 数论矩阵快速幂

数论——快速幂,模运算及快速幂求逆元

数论——快速幂(C++)

数论--快速幂,矩阵快速幂

总结——数论:快速幂

数论-快速幂