国产探花免费观看_亚洲丰满少妇自慰呻吟_97日韩有码在线_资源在线日韩欧美_一区二区精品毛片,辰东完美世界有声小说,欢乐颂第一季,yy玄幻小说排行榜完本

首頁 > 學院 > 開發設計 > 正文

高斯消元

2019-11-08 01:40:30
字體:
來源:轉載
供稿:網友

題: bzoj4184 bzoj3168 bzoj4031 bzoj3270 bzoj3601* bzoj3143

bzoj4184

題意:

給出一堆數的插入和刪除順序,詢問每一次操作后,選出某些數異或起來的最大值a1^a2^a3^…^an (數的個數≤500000)

分析:

有插入和刪除操作不好處理,于是我們嘗試經過一波預處理將其轉化為只有插入~ 顯然一個數的出現時間是一段一段的,離散化后記錄每個數當前的個數便可以知道在哪段區間。 由于我們要求出每一個時間點的答案,便用線段樹結構記錄某個數覆蓋的區間,在最后的dfs時不斷插入到線性基中,到線段樹的葉節點時求出答案。

//核心代碼int getans(base nw){//查詢最大值 int res=0; for (int i=30;i>=0;i--) if (nw.v[i]) res^=nw.v[i]; return res;}void ins(base &nw,int x){//向線性基中插入x int ps=-1; for (int i=30;i>=0;i--)//從最高位開始貪心插入 if ((x>>i)&1){ if (!nw.v[i] && ps==-1) ps=i; else x^=nw.v[i]; } if (~ps){ nw.v[ps]=x; for (int i=30;i>ps;i--) if ((nw.v[i]>>ps)&1) nw.v[i]^=x; }}void dfs(int rt,int l,int r,base nw){ for (int i=h[rt];i;i=g[i].next) ins(nw,g[i].y); if (l==r){ ans[l]=getans(nw); return; } int mid=(l+r)>>1; dfs(2*rt,l,mid,nw); dfs(2*rt+1,mid+1,r,nw);}

bzoj3168

題意:

給出一個矩陣A和B,滿足兩矩陣大小相同。求一個字典序最小的1...n的排列a滿足將隨意一個—Ai換成Bai后矩陣A仍然滿秩。(n≤300

分析:

容易想到二分圖,但問題在于如何在O(n3)內完成構圖。 考慮一個行向量Bi,我們在A中找到最小的行向量集合滿足Bi能夠被這些行向量線性表出,那么顯然Bi僅僅能替換這些行向量

我們能夠設矩陣C滿足C?A=B,那么C=B?A?1 Ci,j≠0表示Bi的線性表出須要Aj,因此CT就是這個二分圖的鄰接矩陣

如今我們有了一個二分圖,怎樣求字典序最小的完備匹配呢?

我們能夠枚舉每一條邊,然后推斷剩余的圖是否存在一個完備匹配。可是這樣做是O(n4)

我們能夠跑兩遍匈牙利算法,第一遍求出隨意一個完備匹配。第二遍對于每一個點貪心選最小的出邊推斷是否能找到不影響前面點的交錯環

總時間復雜度O(n3)

Matrix Get_Inv(Matrix a)//矩陣求逆 { Matrix re(true);//單位矩陣 int i,j,k; for(i=1;i<=n;i++) { for(k=i;k<=n;k++) if(a[k][i]) break; for(j=1;j<=n;j++) { swap(a[i][j],a[k][j]); swap(re[i][j],re[k][j]); } long long inv=Quick_Power(a[i][i],MOD-2); for(j=1;j<=n;j++) { a[i][j]=a[i][j]*inv%MOD; re[i][j]=re[i][j]*inv%MOD; } for(k=1;k<=n;k++) if(k!=i) { long long temp=(MOD-a[k][i])%MOD; for(j=1;j<=n;j++) { (a[k][j]+=a[i][j]*temp%MOD)%=MOD; (re[k][j]+=re[i][j]*temp%MOD)%=MOD; } } } return re; }

bzoj4031

題意:

給出一個”.”“#”矩陣,”.”與”.”之間能連邊,使得任意兩點之間有且只有一條路徑相連。求方案總數。(n,m≤9)

分析:

顯然是矩陣樹定理。給矩陣每一個點標號,若兩相鄰點x,y(編號),x能**向**y連邊就將行列式a[x][x]++,a[x][y]?? 最后求一波行列式的值即可。(將行列式左下角的值全用高斯消元消成0,對角線乘積就是答案) 注意:取模高斯消元的奇怪姿勢

long long slove(){ for (int i=1;i<tot;i++) for (int j=1;j<tot;j++) a[i][j]=(a[i][j]+mod)%mod; long long ans=1,f=1; for (int i=1;i<tot;i++){ for (int j=i+1;j<tot;j++){ long long A=a[i][i],B=a[j][i]; while (B!=0){//用歐幾里德gcd求法理解 long long t=A/B;A%=B; swap(A,B); for (int k=i;k<tot;k++) a[i][k]=(a[i][k]-t*a[j][k]%mod+mod)%mod; for (int k=i;k<tot;k++) swap(a[i][k],a[j][k]); f=-f; } } if (!a[i][i]) return 0; ans=ans*a[i][i]%mod; } return (ans*f+mod)%mod;}

bzoj3270

題意:

給出一個無向圖,兩個人分別從s,t出發,在第i個點有p[i]的幾率停留,不停留則等幾率地跑到相鄰的點上,詢問兩人同時到同一點的幾率是多少。輸出所以點的幾率。(點數≤20

分析:

從最終狀態入手,也就是說當兩個人同時到達(x1,x2),顯然這個狀態是從其他狀態轉換過來的,也就是相鄰的點。 其狀態可能為(x1,x2)、(x1,PRe2)、(pre1,x2)、(pre1,pre2) 這4種狀態轉換到(x1,x2)的概率分別為 p[x1]?p[x2]、p[x1]?(1?p[x2]、(1?p[x1])?p[x2]、(1?p[x1])?(1?p[x2]) 于是每一個點都有一條表示它的答案的表達式。也就是說我們有n個點n條表達式(即方程),再用高斯消元亂搞就好了。

bzoj3143(水)

題意:

一個無向連通圖,頂點從1編號到N,邊從1編號到M。 小Z在該圖上進行隨機游走,初始時小Z在1號頂點,每一步小Z以相等的概率隨機選 擇當前頂點的某條邊,沿著這條邊走到下一個頂點,獲得等于這條邊的編號的分數。當小Z 到達N號頂點時游走結束,總分為所有獲得的分數之和。 現在,請你對這M條邊進行編號,使得小Z獲得的總分的期望值最小。

分析:

給m條邊編號,顯然期望經過次數越多的邊編號越小。 設g[i]為第i條經過的期望次數。 經過第i條邊必然是從第i條邊的一端(X)到另一端(Y)。 也就是說我們要知道經過每個點的期望次數。

f[i]為第i個經過的期望次數。 顯然:f[i]+=f[i?>next]deg[i?>next](deg[i]為第i個點的出度) 于是我們有n個未知數f[]與n條不等式f[i]=...... 那么高斯消元亂搞就好啦~~~

嗯,求出f[i]以后。現在要求的就是g[]了。 顯然:g[i]=f[x]deg[x]+f[y]deg[y](x,y為第i條邊的兩端)

void guass(){//高斯消元模板 for (int i=1;i<=n;i++){ int id=i;ld mx=-1.0; for (int j=i;j<=n;j++) if (fabs(a[j][i])+eps>mx) mx=fabs(a[j][i])+eps,id=j; if (id>n) continue; if (id!=i) swap(a[i],a[id]); for (int j=i+1;j<=n;j++) if (fabs(a[j][i])>eps){ ld t=a[j][i]/a[i][i]; for (int k=i;k<=n+1;k++) a[j][k]-=t*a[i][k]; } } for (int i=n;i;i--){ for (int j=i+1;j<=n;j++) a[i][n+1]-=a[i][j]*a[j][n+1]; a[i][n+1]/=a[i][i]; }}
發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
主站蜘蛛池模板: 焉耆| 平阴县| 千阳县| 惠来县| 尚义县| 晋宁县| 西乡县| 揭阳市| 梨树县| 会泽县| 罗平县| 蒲江县| 衡南县| 玛曲县| 宁晋县| 永善县| 屏东市| 河东区| 孟连| 塔城市| 卫辉市| 恩施市| 连城县| 九江市| 忻州市| 五寨县| 诏安县| 阿拉尔市| 大关县| 仁布县| 建德市| 道孚县| 凉山| 东丰县| 万荣县| 泊头市| 枝江市| 舒城县| 荆门市| 嵩明县| 彭山县|