【BZOJ2751】【codevs1853】容易题，快速幂+逆元

Time:2016.06.24
Author:xiaoyimi
转载注明出处谢谢

---

传送门1
传送门2
思路：
k=0时答案就是
$\prod^m_{i=1}\sum^n_{j=1}j$
$=(\frac{n(n+1)}2)^m$
每一种限制条件就是在乘的时候去除某个数
就比如说n=3，m=3时，答案就是
$(1+2+3)(1+2+3)(1+2+3)=216$
当k=1,限制条件为(1,1)时我们就把第一个括号里的1去掉，变成
$(2+3)(1+2+3)(1+2+3)=180$
k=2，限制条件为(1,1)(2,1)时答案就是
$(2+3)(2+3)(1+2+3)=150$
搞清楚这个原理接下来就很简单了
将k个限制条件排序，然后分开计算每一个位置里的数就好了，剩下的没被去掉的位置用快速幂
（思想如上，但代码表达上有点怪，仅供参考）

#include<cstdio>
#include<iostream>
#include<algorithm>
#define M 100004
#define LL long long
#define mo 1000000007
using namespace std;
int in()
{
    int t=0;char ch=getchar();
    while (ch>'9'||ch<'0') ch=getchar();
    while (ch>='0'&&ch<='9') t=(t<<3)+(t<<1)+ch-48,ch=getchar();
    return t;
}
LL qr(LL x,int y)
{
    LL ans=1;
    for (;y;y>>=1,x=x*x%mo)
        if (y&1) ans=ans*x%mo;
    return ans;
}
int n,m,k;
struct node{
    int x,y;
    bool operator <(node other)const{
        if(x!=other.x) return x<other.x;
        return y<other.y;
    } 
}Q[M];
LL ans,sum,inv;
main()
{
    n=in();m=in();k=in();
    for (int i=1;i<=k;i++)
        Q[i].x=in(),Q[i].y=in();
    LL nn=(LL)n*(n+1)/2%mo; 
    ans=qr(nn,m);
    inv=qr(nn,mo-2);
    sort(Q+1,Q+k+1);
    sum=Q[1].y;
    for (int i=2;i<=k;i++)
    {
        if (Q[i].x==Q[i-1].x)
            if (Q[i].y!=Q[i-1].y) sum+=Q[i].y;
            else;
        else
            ans=ans*(nn-sum)%mo*inv%mo,
            sum=Q[i].y;
    }
    ans=ans*(nn-sum)%mo*inv%mo;
    printf("%d",ans<0?ans+mo:ans);
}