Problem
最近,菲菲学会了做蛋糕,她做了一个 4 D 4D 4D 的蛋糕送给牛牛,这个蛋糕的大小是: A × B × C × D A×B×C×D A×B×C×D。
蛋糕所有的表面都抹着奶油,牛牛想把蛋糕沿着与表面平行的超平面,切成 1 × 1 × 1 × 1 1×1×1×1 1×1×1×1 的小块。
现在,牛牛想知道,在这些小块蛋糕中,有 0 0 0 个、 1 1 1 个、 2 2 2 个 … … …… …… 8 8 8个表面抹着奶油的有分别有多少块。
有 T T T 组数据。
数据范围: T ≤ 10000 T\le10000 T≤10000, 1 ≤ A , B , C , D ≤ 1 0 6 1\le A,B,C,D\le10^6 1≤A,B,C,D≤106。
Solution
我们先从三维找找规律。
我们给每个方块一个坐标 ( a , b , c ) (a,b,c) (a,b,c),不难发现的是,当 a = 1 a=1 a=1 或 a = A a=A a=A 时( b , c b,c b,c 同理),这一层是涂了奶油的。
那么换到四维,很容易得到以下式子:
ans[0]=(a-2)*(b-2)*(c-2)*(d-2);
ans[1]=2*((a-2)*(b-2)*(c-2)+(a-2)*(b-2)*(d-2)+(a-2)*(c-2)*(d-2)+(b-2)*(c-2)*(d-2));
ans[2]=4*((a-2)*(b-2)+(a-2)*(c-2)+(a-2)*(d-2)+(b-2)*(c-2)+(b-2)*(d-2)+(c-2)*(d-2));
ans[3]=8*((a-2)+(b-2)+(c-2)+(d-2));
ans[4]=16;
也可以用降维的方法(代码中用的这个),不过思想都是一样的。
要注意 A = 1 A=1 A=1 等的情况。
Code
#include<cstdio>
#include<cstring>
#include<algorithm>
#define P 2148473648
#define ll long long
using namespace std;
ll a[10],ans[10];
void dfs(int x,ll val,int pos){
if(x==5){
ans[pos]=(ans[pos]+val)%P;
return;
}
if(a[x]==1){
dfs(x+1,val,pos);
return;
}
dfs(x+1,val*2%P,pos-1);
dfs(x+1,val*(a[x]-2)%P,pos-2);
}
int main(){
int T;
scanf("%d",&T);
while(T--){
memset(ans,0,sizeof(ans));
for(int i=1;i<=4;++i) scanf("%lld",&a[i]);
dfs(1,1,8);
for(int i=0;i<=8;++i) printf("%lld ",ans[i]);
puts("");
}
return 0;
}