Description

給nn個(gè)小于pp的非負(fù)整數(shù)a1,…,na1,…,n，問(wèn)有多少對(duì)(i,j)(1≤i<j≤n)(i,j)(1≤i<j≤n)模pp在意義下滿足1ai+aj≡1ai+1aj1ai+aj≡1ai+1aj，即這兩個(gè)數(shù)的和的逆元等于這兩個(gè)數(shù)的逆元的和，注意0沒(méi)有逆元

Input

第一行一整數(shù)TT表示用例組數(shù)，每組用例首先輸入一整數(shù)nn表示序列長(zhǎng)度和一素?cái)?shù)pp表示模數(shù)，之后輸入nn個(gè)非負(fù)整數(shù)a1,…,n(1≤T≤5,1≤n≤2×105,2≤p≤1018,0≤a1,…,n<p)a1,…,n(1≤T≤5,1≤n≤2×105,2≤p≤1018,0≤a1,…,n<p)
Output

輸出滿足條件的(i,j)(1≤i<j≤n)(i,j)(1≤i<j≤n)對(duì)數(shù)

Sample Input

2
5 7
1 2 3 4 5
6 7
1 2 3 4 5 6

Sample Output

4
6


最后我明白了個(gè)道理，當(dāng)?shù)讛?shù)過(guò)大時(shí)，不能用普通乘法，更不不能用快速冪，因?yàn)槌艘槐榫捅?。于是釀成慘劇！

#include <bits/stdc++.h> using namespace std; typedef long long ll; #define _m(a, p) make_pair(a, p) map<ll, ll> mp; map<ll, ll> mm; ll mod; long long ksc(long long a, long long b, long long mod) {long long ans = 0;for (; b; b >>= 1){if (b & 1)ans = (ans + a) % mod;a = (a + a) % mod; //（計(jì)算機(jī)加法比乘法快，a+a比a*2快）}return ans; } int main() {int t, n;scanf("%d", &t);while (t--){ll cnt = 0;mp.clear();mm.clear();scanf("%d %lld", &n, &mod);for (int i = 0; i < n; i++){ll a;scanf("%lld", &a);if (!a)continue;mm[a]++;ll ans = ksc(ksc(a, a, mod), a, mod);mp[ans]++;}for (auto p : mp){ll cc = p.second;cnt += cc * (cc - 1) / 2;}if (mod != 3)for (auto m : mm){ll n = m.second;cnt -= n * (n - 1) / 2;}printf("%lld\n", cnt);}return 0; } 創(chuàng)作挑戰(zhàn)賽新人創(chuàng)作獎(jiǎng)勵(lì)來(lái)咯，堅(jiān)持創(chuàng)作打卡瓜分現(xiàn)金大獎(jiǎng)

總結(jié)

以上是生活随笔為你收集整理的数学--数论--HDU 6128 Inverse of sum （公式推导论）的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。