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;
}

/*

*/