bzoj 2751(快速幂+乘法原理)

传送门
没时间写题解了。。。hfu老师要求22点前总结计划,大概就是表示出一般情况的答案(乘法原理),然后有限制/无限制分开计算,最后乘在一起。
转一篇别人的题解:http://www.cnblogs.com/Beckinsale/p/7488874.html

#include
using namespace std;
typedef long long ll;
const int MAXN=1e5+4;
const ll MOD=1e9+7;
struct Limit {
    int x,y;
    friend bool operator <(const Limit &a,const Limit &b) {
        return a.x==b.x?a.yint k;
ll n,m,s[MAXN],sum,ans=0,cnt=0;
ll fpow(ll a,ll b,ll p) {
    ll ret=1;
    while (b) {
        if (b&1) ret=ret*a%p;
        b>>=1,a=a*a%p;
    }
    return ret;
}
int main() {
    scanf("%lld%lld%d",&n,&m,&k);
    sum=(n*(n+1)>>1)%MOD;
    for (int i=1;i<=k;++i) scanf("%d%d",&a[i].x,&a[i].y);
    sort(a+1,a+k+1);
    for (int i=1;i<=k;++i) {
        if (a[i].x^a[i-1].x) s[++cnt]=sum;
        else if (a[i].y==a[i-1].y) continue;
        s[cnt]=(s[cnt]-a[i].y+MOD)%MOD;
    }
    ans+=fpow(sum,m-cnt,MOD);
    for (int i=1;i<=cnt;++i) ans=ans*s[i]%MOD;
    printf("%lld\n",ans);
    return 0;
}

你可能感兴趣的:(数论)