bzoj 2751 快速幂

首先我们知道，对于所有种情况，我们可以将每一位可以放的

数的值加起来，所有位置的乘起来，等于的就是最后的答案，具体

为什么正确，可以根据乘法分配律来想一想。

那么对于所有不做要求的，快速幂直接算就行了，然后快排下，就知道

每个位置不放那些值，减掉后乘进去就行了。

/**************************************************************
    Problem: 2751
    User: BLADEVIL
    Language: Pascal
    Result: Accepted
    Time:344 ms
    Memory:1008 kb
****************************************************************/
 
//By BLADEVIL
const
    d39                                 =1000000007;
     
var
    n, m, k                             :longint;
    a, b                                :array[0..100010] of longint;
    ans                                 :int64;
     
procedure swap(var a,b:longint);
var
    c                                   :longint;
begin
    c:=a; a:=b; b:=c;
end;
     
procedure qs(low,high:longint);
var
    i, j, xx, yy                        :longint;
begin
    i:=low; j:=high; 
    xx:=a[(i+j) div 2]; yy:=b[(i+j) div 2];
    while i<j do
    begin
        while (a[i]<xx) or (a[i]=xx) and (b[i]<yy) do inc(i);
        while (a[j]>xx) or (a[j]=xx) and (b[j]>yy) do dec(j);
        if i<=j then
        begin
            swap(a[i],a[j]);
            swap(b[i],b[j]);
            inc(i); dec(j);
        end;
    end;
    if i<high then qs(i,high);
    if j>low then qs(low,j);
end;
     
procedure init;
var
    i                                   :longint;
begin
    read(n,m,k);
    for i:=1 to k do read(a[i],b[i]);
    qs(1,k);
end;
 
function mi(a,b:int64):int64;
var
    sum                                 :int64;
begin
    sum:=a;
    mi:=1;
    while b<>0 do
    begin
        if b mod 2=1 then mi:=mi*sum mod d39;
        sum:=sum*sum mod d39;
        b:=b div 2;
    end;
end;
 
procedure main;
var
    i                                   :longint;
    sum, x, y, z                        :int64;
begin
    sum:=m;
    x:=-1;
    for i:=1 to k do
    begin
        if a[i]<>x then
        begin
            dec(sum);
            x:=a[i];
        end;
    end;
    x:=n; y:=n+1;
    if x mod 2=0 then x:=x div 2 else y:=y div 2;
    x:=x mod d39;
    y:=y mod d39;
    x:=x*y mod d39;
    ans:=mi(x,sum);
    for i:=1 to k do if (a[i]=a[i-1]) and (b[i]=b[i-1]) then b[i-1]:=0;
    y:=-1;
    z:=-1;
    for i:=1 to k do
    begin
        if a[i]<>y then
        begin
            if i<>1 then ans:=ans*z mod d39;
            z:=x;
            y:=a[i];
            z:=((x-b[i]) mod d39+d39) mod d39;
        end else z:=((z-b[i])mod d39+d39) mod d39;
    end;
    if z<>-1 then ans:=ans*z mod d39;
    writeln(ans);
end;
     
 
begin
    init;
    main;
end.