BZOJ4513 [Sdoi2016]储能表

对i^j打个表，可以发现很多规律，具体规律可以看这里

然后我们对于一个n*m的矩形，找到小于n和m中较大值的最大的2的幂，将矩形划分为2部分或者4部分，其中1部分或者3部分都能直接计算，递归剩下那部分即可

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<algorithm>
#include<iomanip>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<bitset>
#include<set>
#include<map>
using namespace std;
#define MAXN 1010
#define MAXM 1010
#define INF 1000000000
#define eps 1e-8
#define MOD 1000000007
#define ll long long
ll n,m,k,p;
ll s(ll n,ll k){
ll re=0;
if(k<0){
k=-k;
re+=((n+1)%p)*(k%p)%p;
k=0;
}
if(k>n){
return 0;
}
n-=k;
if(n&1){
(re+=(n%p)*((n+1)/2%p)%p)%=p;
}else{
(re+=(n/2%p)*((n+1)%p)%p)%=p;
}
return re;
}
ll cal(ll n,ll m,ll k){
if(n<0||m<0){
return 0;
}
ll i;
if(n<m){
swap(n,m);
}
for(i=1;i-1<=n;i<<=1){

}
i>>=1;
i--;
if(m>=i){
return ((i+1)%p*s(i,k)%p+(m-i+n-i)%p*s(i,k-i-1)%p+cal(n-i-1,m-i-1,k))%p;
}else{
return ((m+1)%p*s(i,k)%p+cal(n-i-1,m,k-i-1))%p;
}
}
int main(){
int tmp;
scanf("%d",&tmp);
while(tmp--){
scanf("%lld%lld%lld%lld",&n,&m,&k,&p);
printf("%lld\n",cal(n-1,m-1,k));
}
return 0;
}

/*

*/