日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

【牛客 - 181F】子序列(容斥,组合数,费马小定理)

發布時間:2023/12/10 编程问答 40 豆豆
生活随笔 收集整理的這篇文章主要介紹了 【牛客 - 181F】子序列(容斥,组合数,费马小定理) 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

?

題干:

題目描述

給出一個長度為n的序列,你需要計算出所有長度為k的子序列中,除最大最小數之外所有數的乘積相乘的結果

?

輸入描述:

第一行一個整數T,表示數據組數。 對于每組數據,第一行兩個整數N,k,含義如題所示接下來一行N個整數,表示給出的序列保證序列內的數互不相同

輸出描述:

對于每組數據,輸出一個整數表示答案,對取模 每組數據之間以換行分割

?

示例1

輸入

復制

3 4 3 5 3 1 4 5 4 3 7 5 2 1 10 3 100 1020 2050 102 12 235 4 57 32135 54354

輸出

復制

144 81000 521918013

說明

?

第一組數據解釋

所有長度為3的子序列為

最終答案為

備注:

對于的數據: 對于的數據: 對于的數據: 保證序列中的元素互不相同且,

解題報告:

首先暴力可以n^2。思路就是要算每一個數的貢獻,也就是每一個數的出現次數。

總次數是C(n-1,k-1)。我們可以枚舉???在前面選了多少個,在后面選了多少個。但是這樣復雜度就成了O(nk)了,,也高不了哪去。所以想個辦法優化,不難想到,可以求問題的反面,容斥一下,用總次數減去作為最大最小值次數。可以排個序,和前面搭配是作為最大值,和后面搭配作為最小值。

注意,以前組合數是乘數,現在組合數是指數,不要習慣性對p=1e9+7取模。

指數可以對p-1取模。因為這里p是質數,a^(p-1)=1 mod p(費馬小定理)

所以,指數減掉若干個p-1,并不影響。

當一般地,(a,p)=1而p不是質數,就歐拉定理,a^(phi(p))=1 mod p所以對phi(p)取模

AC代碼:

#include<bits/stdc++.h> #define int long long using namespace std; const int N=2005; const int mod=1e9+7; typedef long long ll; int n,q,k; ll a[N]; ll c[N][N]; ll qm(ll x,ll y) {ll ret=1;while(y) {if(y%2==1) ret=(ret*x)%mod;x=(x*x)%mod;y/=2;}if(ret<0) ret=(ret+mod)%mod;return ret%mod; } signed main() {for(int i=0; i<=2001; i++) {c[i][0]=1;for(int j=1; j<=i; j++) {c[i][j]=(c[i-1][j-1]+c[i-1][j])%(mod-1);}}scanf("%lld",&q);while(q--) {scanf("%lld%lld",&n,&k);for(int i=1; i<=n; i++)scanf("%lld",&a[i]);sort(a+1,a+n+1);ll ans=1;for(int i=1; i<=n; i++) {ll ci=(c[n-1][k-1]-c[i-1][k-1]-c[n-i][k-1]+5*(mod-1))%(mod-1);ans=(ans*qm(a[i],ci))%mod;}printf("%lld\n",ans);}return 0; }

參考博客:鏈接???????

總結

以上是生活随笔為你收集整理的【牛客 - 181F】子序列(容斥,组合数,费马小定理)的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。