SSL 2676_数学_模拟

Posted

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了SSL 2676_数学_模拟相关的知识,希望对你有一定的参考价值。

题目描述

小A得到了一个数列A,他告诉你这个数列中所有的A[i]都是1到n的自然数,并且告诉你对于一些A[i]不能取哪些值。无聊的你想要知道所有可能的数列的积的和。定义一个数列的积为这个数列所有数的乘机。由于答案太大,只要模10^9+7输出即可。


 

思路

首先有一个规律就是ans为每一项中没有被限制的数相加让后相乘

显然这样直接暴力会T,数据中的限制有 $10^5$ 个,而总数则有$10^9$ 个

所以没有限制条件的数量还是很多的,我们就可以用一个快速幂加上暴力来求解


#include <stdio.h>
#include <algorithm>
#include <iostream>
using namespace std;
#define p 1000000007
#define maxn 100001
struct arr
{
    long long x, y;
}a[maxn];
long long tot;
long long read()
{
    long long x = 0; char ch = getchar();
    while (ch < 0 || ch > 9) ch = getchar();
    while (ch >= 0 && ch <= 9) {x = (x << 1) + (x << 3) + ch -0; ch = getchar();}
    return x;
}
long long cmp(arr a, arr b)
{
    return (a.x < b.x) || (a.x == b.x && a.y < b.y);
}
long long fast(long long x)
{
    if (x == 0) return 1;
    if (x == 1) return tot;
    long long t = fast(x / 2) % p;
    if (x % 2 == 1) return (((t * t) % p) * (tot % p)) % p;
    else return (t * t) % p;

}
int main()
{
    long long n, m ,k;
    cin >> n >> m >> k;
    for (long long i = 1; i <= k; i++)
    {
        a[i].x = read();
        a[i].y = read();
    }
    sort(a + 1, a + k + 1, cmp);
    for (long long i = 1; i < k; i++)
        if (a[i].x == a[i + 1].x && a[i].y == a[i + 1].y) a[i].y = 0;
    tot = ((1 + n) * n) / 2;
    long long i = 0, ans = 1, now = 1, sum = m;
    while (i < k)
    {
        i++;
        if (a[i].x != a[i - 1].x)
        {
            ans = (ans * (now % p)) % p;
            sum--;
            now = tot;
        }
        if (a[i].y != 0) now -= a[i].y;
    }
    ans = (ans * (now % p)) % p;
    ans = (ans * fast(sum)) % p;
    cout << ans << endl;
}

 

以上是关于SSL 2676_数学_模拟的主要内容,如果未能解决你的问题,请参考以下文章

c# https请求忽略证书验证_各种编程语言忽略http的SSL证书认证

某某天的模拟赛题 IOer(组合数学)

[CSP-S模拟测试]:小P的生成树(数学+Kruskal)

[CSP-S模拟测试]:Divisors(数学)

jzoj6004. PKUWC2019模拟2019.1.17集合 (组合数学)

[数学] aw3815. 最大约数(模拟+分解质因数+CF588B)