• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            算法學社
            記錄難忘的征途
            posts - 141,comments - 220,trackbacks - 0
            吐槽:
                1. 最近因為一些考試/數(shù)據(jù)結(jié)構(gòu)試驗報告,等瑣事連續(xù)N天沒有做題,也就沒有更新blog了。真是罪過罪過…………
                2. 600pt 代碼偏長(T_T) 最后時刻沒有調(diào)出來,等考完試奉獻上代碼。
                3. 和HIT的學長分到一個房間了Orz....
            250pt
            有兩個集合,每個集合有N(N<50)個錐面(表示為高H和底面半徑R,均小于10,000)。要求從集合1和集合2的笛卡爾積集合中選取最多的元素滿足如下性質(zhì):
                1. Ha / Ra > Hb / Rb
                2. Ra < Rb
            算法分析:
                裸二分匹配,貌似貪心也可
            #include<iostream>
            #include<vector>
            using namespace std;
            int G[55][55];
            int n,m,yM[55],chk[55];
            int dfs(int u){
                for(int v=0; v<m;v++) if(G[u][v] && !chk[v]){
                    chk[v] = 1;
                    if(yM[v]==-1 || dfs(yM[v])){
                        yM[v] = u; return 1;
                    }
                }
                return 0;
            }
            bool judge(int th,int tr, int bh, int br){
                return br > tr && br*th >tr*bh;
            }
            class PointyWizardHats{
                public : int getNumHats(vector <int> th, vector <int> tr, vector <int> bh, vector <int> br){
                    n = th.size(), m = bh.size();
                    for(int i=0;i<n;i++)
                        for(int j=0;j<m;j++)
                            if(judge(th[i],tr[i],bh[j],br[j]))
                                G[i][j] =1;
                    int ans=0;
                    for(int i=0;i<m;i++) yM[i] = -1;
                    for(int i=0;i<n;i++){
                        for(int j=0;j<m;j++) chk[j] = 0;
                        ans += dfs(i);
                    }
                    return ans;
                }
            };
            600pt
            A和B博弈。 在N*N(N<=13)的矩陣上,有一些 ’H‘ 格子, 一開始每個H格子里最多有一枚帶權(quán)值的硬幣。 每次小B選一個H格子,然后小A將硬幣任意排列后,將該格子的硬幣給小B。
            排列后必須滿足如下條件:
                1. 每個格子必須有最多一枚硬幣,小B選的格子可以沒有硬幣
                2. 每一行,H的個數(shù)和硬幣的個數(shù)之和必須是偶數(shù)
                3. 每一列也一樣。
            小B想讓硬幣的權(quán)值和最多,小A反之。
            小B選擇一個格子后,那個格子上就不能在轉(zhuǎn)移硬幣了。如果有就永遠有,沒有就永遠沒有。
            問兩人都用最優(yōu)策略,B得到的分數(shù)和是多少。
            算法分析:
                博弈DP: dp[mask] mask是三進制,1表示已經(jīng)選擇了那個格子,2表示那個格子上有硬幣。
                之前要預(yù)處理一遍,檢查某個狀態(tài)是否是合法的。
            #include<string>
            #include<iostream>
            #include<vector>
            #include<algorithm>
            using namespace std;
            const int N = 13;
            pair<int,int> pos[N];
            bool can[1600000],row[N],col[N];
            int count(int mask){
                int a = 0;
                while(mask) a += mask & 1, mask >>= 1;
                return a;
            }
            int dp[1600000],base[20],cont[1<<N],n,cns;
            int cal(int a,int b){
                int ans=0;
                for(int i=0;i<n;i++){
                    ans += base[i] * ((a&1) + (b&1));
                    a>>=1, b>>=1;
                }
                return ans;
            }
            const int inf = 100;
            void OP(int msk){for(int i=0;i<n;i++,msk>>=1)cout<<(msk&1);cout<<" ";}
            int dfs(int msk1,int msk2,int m){
                
                if(m==0) return 0;
                int &ans = dp[cal(msk1,msk2)];
                if(ans!=-1) return ans;
                for(int i=0;i<n;i++)if((1<<i) & msk1);else{
                    int nmask = msk1 ^ (1<<i),l=inf,r=inf;
                    if(can[cal(nmask,msk2)]) l = dfs(nmask,msk2,m-1);
                    if(can[cal(nmask,msk2 ^ (1<<i))]) r = dfs(nmask,msk2 ^ (1<<i),m-1)+1;
                    if(l > r) l = r;
                    if(l>=inf) continue;
                    if(ans < l) ans = l;
                }
                return ans;
            }
            class MagicalHats{
                public : 
                int findMaximumReward(vector <string> G, vector <int> coins, int m){
                    n = 0;
                    int flag,k, ans = 0;cns = coins.size();
                    memset(dp,-1,sizeof(dp));
                //    memset(hash,-1,sizeof(hash));
                    for(int i=0;i<G.size();i++)
                        for(int j=0;j<G[0].size();j++)
                            if(G[i][j] == 'H'){
                                pos[n++] = make_pair(i,j);
                            }
                    base[0] = 1;
                    for(int i=1;i<n;i++) base[i] = 3*base[i-1];
                    int l1=G.size(),l2=G[0].size();
                    for(int i=0;i<(1<<n);i++) cont[i] = count(i);
                    for(int msk1 = 0; msk1 < (1<<n); msk1 ++){
                        for(int msk2 = 0; msk2 < (1<<n); msk2 ++){
                            if(cont[msk2] != cns) continue;
                            if(can[k=cal(msk1,msk2 & msk1)]) continue;
                            for(int i=0;i<l1;i++) row[i] = 0;
                            for(int i=0;i<l2;i++) col[i] =0;
                            for(int i=0;i<n;i++) if((1<<i) & msk2);else row[pos[i].first]^=1, col[pos[i].second] ^=1;
                            flag = 0;
                            for(int i=0;i<G.size();i++) flag |= row[i];
                            for(int i=0;i<G[0].size();i++) flag |= col[i];
                            if(flag) continue;
                            can[k] = 1;
                        }
                        if(msk1 && !can[0]) return -1;
                    }
                //    cout<<"."<<endl;
                    if(!can[0]) return -1;
                    k = dfs(0,0,m);
                    sort(coins.begin(), coins.end());
                    for(int i=0;i < k; i++) ans += coins[i];
                    return ans;
                }
            };
            posted on 2012-07-09 21:49 西月弦 閱讀(452) 評論(0)  編輯 收藏 引用 所屬分類: 比賽感言
            欧美亚洲国产精品久久高清| 精品久久久久久成人AV| 亚洲狠狠久久综合一区77777| 青青草原综合久久| 久久久WWW成人| 久久亚洲精品国产精品| 91久久婷婷国产综合精品青草 | 国产精品99久久久精品无码| 国产成人精品综合久久久| 久久精品无码一区二区三区| 狠狠色丁香久久婷婷综合图片| 91精品国产高清久久久久久io | 久久天天躁狠狠躁夜夜不卡 | 久久强奷乱码老熟女网站| 国产精品一久久香蕉产线看| 无码任你躁久久久久久久| 久久久91精品国产一区二区三区| 综合久久精品色| 久久精品二区| 久久精品成人国产午夜| 亚洲中文精品久久久久久不卡| 麻豆久久| 久久久91人妻无码精品蜜桃HD| 国产成人久久777777| 久久久久99精品成人片试看| 国产激情久久久久久熟女老人| 免费精品国产日韩热久久| 国产精品热久久毛片| 久久综合欧美成人| 精品久久久久久中文字幕人妻最新| 久久99热这里只有精品国产| 久久热这里只有精品在线观看| 亚洲国产日韩欧美综合久久| 中文字幕无码久久久| 欧美精品一区二区久久| 久久男人AV资源网站| 日日狠狠久久偷偷色综合0| 日韩中文久久| 国产精品美女久久福利网站| 亚洲国产精品一区二区三区久久| 亚洲欧美一区二区三区久久|