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

歡迎訪問 生活随笔!

生活随笔

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

编程问答

bzoj2751[HAOI2012]容易题(easy)

發布時間:2023/12/9 编程问答 26 豆豆
生活随笔 收集整理的這篇文章主要介紹了 bzoj2751[HAOI2012]容易题(easy) 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

bzoj2751[HAOI2012]容易題(easy)

題意:

已知一個數列A對于所有的A[i]都是1~n的自然數,一些A[i]不能取一些值,求出所有可能的數列的積的和 mod 1000000007的值。

題解:

題目中的n≤109實際上是109……首先推個方程s[l,r]=s[l,k]*s[k+1,r](s[l,r]表示l到r的所有l≤i≤r的a[i]的可能取值的和)因此s[1,n]等于所有a[i]的可能取值的和的乘積。因此我們先求出1到n的和,對每個約束條件按i排序,將這個和減掉約束條件中的不能取的數,就是這個a[i]所有可能取值的和。將這些a[i]乘起來,剩下的沒限制的a[i]用快速冪解決。

代碼:

1 #include <cstdio> 2 #include <cstring> 3 #include <algorithm> 4 #define ll long long 5 #define inc(i,j,k) for(int i=j;i<=k;i++) 6 #define mod 1000000007 7 using namespace std; 8 9 struct nd{ 10 ll a,b; 11 bool operator < (const nd &c)const{ 12 if(a!=c.a)return a<c.a;else return b<c.b; 13 } 14 }; 15 ll power(ll a,ll b){ 16 if(b==0)return 1; if(b==1)return a; ll c=power(a,b>>1)%mod; 17 if(b&1)return c*c%mod*a%mod;else return c*c%mod; 18 } 19 nd ns[200000]; 20 int main(){ 21 ll n,m,a1=0,a2,a3,a4; ll k; scanf("%lld%lld%lld",&n,&m,&k); 22 inc(i,1,k)scanf("%lld%lld",&ns[i].a,&ns[i].b); sort(ns+1,ns+k+1); 23 inc(i,1,k)if(i==1||ns[i].a!=ns[i-1].a)a1++; a2=n*(n+1)/2%mod; a3=a4=1; 24 inc(i,1,k)if(i==1||ns[i].a!=ns[i-1].a)a4=a4*a3%mod,a3=a2,a3=(a3-ns[i].b)>=0?(a3-ns[i].b)%mod:(a3-ns[i].b)+mod; 25 else if(ns[i].b!=ns[i-1].b)a3=(a3-ns[i].b)>=0?(a3-ns[i].b)%mod:(a3-ns[i].b)+mod; 26 a4=a4*a3%mod; a4=a4*power(a2,m-a1)%mod; printf("%lld",a4); 27 }

?

20160419

轉載于:https://www.cnblogs.com/YuanZiming/p/5703229.html

總結

以上是生活随笔為你收集整理的bzoj2751[HAOI2012]容易题(easy)的全部內容,希望文章能夠幫你解決所遇到的問題。

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