2019牛客多校Monotonic Matrix
Monotonic Matrix
題意:
問(wèn)有多少個(gè)n * m的矩陣A滿足一下情況:答案mod 1e9+7
題解:
我們先看看這個(gè)式子有啥規(guī)律?
對(duì)于每一行,每一列都是非下降序列,也就是說(shuō)如果有一位是2,其后面都是2,如果有一位是1,其前面都是1
那我們考慮元素01和12的分界線
是(n,0)到(0,m)的兩條不相交(可重合的路徑)
第一個(gè)分界線以上的點(diǎn)是一種,兩條分界線之間是一種,第二個(gè)分界線以下是一種(詳細(xì)看圖)
不過(guò)現(xiàn)在這個(gè)兩個(gè)線的起終點(diǎn)一樣,我們可以將第一個(gè)線進(jìn)行偏移,變成起點(diǎn)為(n-1,-1)到(-1,m-1).(注意這兩個(gè)線不是弧形的,而是水平豎直線)
現(xiàn)在問(wèn)題就是起點(diǎn)為(n-1,-1),終點(diǎn)為(-1,m-1)和起點(diǎn)為(n,0),終點(diǎn)為(0,m),兩個(gè)不相交路徑的條數(shù)
怎么做?引入LGV定理
就有:
起點(diǎn){a1,a2a_{1},a_{2}a1?,a2?}={(n,0),(n?1,?1)(n,0),(n-1,-1)(n,0),(n?1,?1)}
終點(diǎn){b1,b2b_{1},b_{2}b1?,b2?}={(0,m),(?1,m?1)(0,m),(-1,m-1)(0,m),(?1,m?1)}
帶入公式:
ans=∣(a1,b1)(a1,b2)(a2,b1)(a2,b2)∣ans= \begin{vmatrix} (a_{1},b_{1})&(a_{1},b_{2})\\ (a_{2},b_{1})&(a_{2},b_{2})\\ \end{vmatrix} ans=∣∣∣∣?(a1?,b1?)(a2?,b1?)?(a1?,b2?)(a2?,b2?)?∣∣∣∣?
(a1,b1)(a_{1},b_{1})(a1?,b1?)表示從(n,0)到(0,m)的路徑數(shù),從(n,0)到(0,m)共有n+m步,我們選擇其中n步向上走,剩下的自然是向右走,共Cn+mnC_{n+m}^{n}Cn+mn?種方法,其他同理
最終答案ans=Cn+mn?Cn+mn?1?Cn+mm?1?Cn+mnC_{n+m}^{n}*C_{n+m}^{n-1}-C_{n+m}^{m-1}-C_{n+m}^{n}Cn+mn??Cn+mn?1??Cn+mm?1??Cn+mn?
代碼:
#include <bits/stdc++.h> #include <unordered_map> #define debug(a, b) printf("%s = %d\n", a, b); using namespace std; typedef long long ll; typedef unsigned long long ull; typedef pair<int, int> PII; clock_t startTime, endTime; //Fe~Jozky const ll INF_ll= 1e18; const int INF_int= 0x3f3f3f3f; void read(){}; template <typename _Tp, typename... _Tps> void read(_Tp& x, _Tps&... Ar) {x= 0;char c= getchar();bool flag= 0;while (c < '0' || c > '9')flag|= (c == '-'), c= getchar();while (c >= '0' && c <= '9')x= (x << 3) + (x << 1) + (c ^ 48), c= getchar();if (flag)x= -x;read(Ar...); } template <typename T> inline void write(T x) {if (x < 0) {x= ~(x - 1);putchar('-');}if (x > 9)write(x / 10);putchar(x % 10 + '0'); } void rd_test() { #ifdef ONLINE_JUDGE #elsestartTime = clock ();freopen("data.in", "r", stdin); #endif } void Time_test() { #ifdef ONLINE_JUDGE #elseendTime= clock();printf("\nRun Time:%lfs\n", (double)(endTime - startTime) / CLOCKS_PER_SEC); #endif } const int maxn=3e6+9; const int mod=1e9+7; ll inv[maxn]; ll f[maxn]; ll f0[maxn]; void init(int N){inv[1]=inv[0]=1;f[1]=f[0]=1;f0[0]=f0[1]=1;for(int i=2;i<=N;i++){f[i]=f[i-1]*i%mod;f0[i]=(mod-mod/i)*f0[mod%i]%mod;inv[i]=inv[i-1]*f0[i]%mod;} } ll C(ll a,ll b){return 1ll*f[a]*inv[b]%mod*inv[a-b]%mod; } int main() {//rd_test();init(10000);ll n,m;while(~scanf("%lld%lld",&n,&m)){cout<<(C(n+m,n)*C(n+m,n)%mod-C(n+m,n-1)*C(n+m,m-1)%mod+mod)%mod<<endl;}//Time_test(); } 創(chuàng)作挑戰(zhàn)賽新人創(chuàng)作獎(jiǎng)勵(lì)來(lái)咯,堅(jiān)持創(chuàng)作打卡瓜分現(xiàn)金大獎(jiǎng)總結(jié)
以上是生活随笔為你收集整理的2019牛客多校Monotonic Matrix的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: 关于图片预加载loading及加载失败的
- 下一篇: U盘加密软件厂商提醒:实现U盘文件防拷贝