題目描述
在三維空間中,平面 x = 0, y = 0, z = 0,以及平面 x + y + z = K 圍成了一個三棱錐。
整天與整數打交道的小明希望知道這個三棱錐內、上整點的數目。
他覺得數量可能很多,所以答案需要對給定的 M 取模。
整天與整數打交道的小明希望知道這個三棱錐內、上整點的數目。
他覺得數量可能很多,所以答案需要對給定的 M 取模。
輸入描述:
輸入有 1 ≤ T ≤ 105
組數據。9
每組數據中,輸入兩個整數 0 ≤ K ≤ 10
+ 7, 1 ≤ M ≤ 109
+ 7,意義如題目描述。
輸出描述:
對於每組數據,輸出一個整數,為三棱錐內、上整點的數目對 M 取模。
輸入
4 0 60 1 60 29 60 29 100007
輸出
1 4 40 4960
題意 : 詢問在三棱錐內上的整點個數。
思路分析 : 兩種方法
在計算除法取模的時候,如果除數不是素數,則不能用快速冪取求逆元解決這個問題,因此這里要借助一個除法取模的性質, (a/b)%m = (a % (b*m))/b
1 、 將三棱錐截成一個平面一個平面的三角形,然后計算每個平面內整點的數量,最后一個邊長的累加就可以計算出來,推得的公式是
代碼示例 :
int main() { //freopen("in.txt", "r", stdin); //freopen("out.txt", "w", stdout); int t; ll n, m; cin >> t; while(t--){ cin >> n >> m; ll x = ((n+1)*(n)/2)%(3*m); x = x*(n+5)%(3*m); printf("%lld\n", (x/3+(n+1)%m)%m); } return 0; }
方法二 : 隔板法
由題意可知 x + y + z <= k , 那么我們可以加一個 d ,將其變成 x + y + z + d = k,其中 x , y , z , d 均可以為 0 ,那么這個問題不久變成了很經典的隔板法問題了嗎,在 4 個桶里都先放入一個物品,所以答案是 C(n+k-1, n-1) , 即 C(n+k-1, 3)
代碼示例 :
int main() { int T; scanf("%d",&T); while(T--) { ll k,m; scanf("%lld%lld",&k,&m); ll res=(k+3)*(k+2)/2%(3*m); res=res*(k+1)%(3*m); printf("%lld\n",res/3); } return 0; }