Description
为了使得大家高兴,小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
HINT
数据范围
30%的数据n<=4,m<=10,k<=10
另有20%的数据k=0
70%的数据n<=1000,m<=1000,k<=1000
100%的数据 n<=109,m<=109,k<=105,1<=y<=n,1<=x<=m
想说的话
这题哪里简单了qwq
题解
这道题我们可以优先考虑每个元素对于我们结果的贡献
每个元素的贡献值都是从1一直加到n 这个推一推就能推出来
然后我们考虑的就是针对没有限制的元素
我们发现这些元素的贡献值是一样的
对于这些元素 我们用一个快速幂取余得到结果
剩下的我们直接把限制元素减掉就是贡献、
然后直接无脑乘进去就好
然后值得注意的一点就是sample特别有良心
他给出了相同的数据
换句话说我们需要去重
但是看到恐怖的数据范围 hash是不可能的了
所以我们采用数据结构map
用bool来映射long long
这样的话就很好办了
内存也不怕 时间也可以得到优化
我们用一个结构体limit来维护我们限制的元素和限制的值
然后用一个cmp结构体排序 去重更方便
题解
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
#define MAXN 100005
#define ll long long
#define MOD 1000000007
using namespace std;
map<ll,bool>st;
inline void read(ll &x)
{
ll s=0,w=1;
char c=getchar();
while(c<'0'||c>'9')
{
if(c=='-')w=-1;
c=getchar();
}
while(c>='0'&&c<='9')
{
s=s*10+c-'0';
c=getchar();
}
x=s*w;
}
struct limit
{
ll x,y;
friend bool operator <(limit a,limit b)
{
if(a.x==b.x)return a.y<b.y;
return a.x<b.x;
}
}lim[MAXN];
ll n,m,k,ans=1,sum,tot,cnt;
ll quick_pow(ll a,ll b)
{
if(!b)return 1;
ll temp=quick_pow(a,b>>1);
temp*=temp;
temp%=MOD;
if(b&1)temp*=a;
return temp%MOD;
}
int main()
{
read(n),read(m),read(k);
sum=(n+1)*n/2;
sum%=MOD;
for(int i=1;i<=k;i++)
read(lim[i].x),read(lim[i].y);
sort(lim+1,lim+k+1);
st.clear();
for(ll i=1;i<=k;i++)
{
if(lim[i].x==lim[i+1].x)
{
if(!st[lim[i].y])
cnt+=lim[i].y,st[lim[i].y]=1;
}
else
{
tot++;
if(!st[lim[i].y])
cnt+=lim[i].y;
ans*=((sum-cnt)%MOD);
ans%=MOD;
cnt=0;
st.clear();
}
}
ans*=quick_pow(sum,m-tot);
ans%=MOD;
printf("%lld\n",(ans%MOD+MOD)%MOD);
}
写在后面
没有2333