• <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>
            隨筆 - 97, 文章 - 22, 評論 - 81, 引用 - 0
            數據加載中……

            PKU 1186 方程的解數

            題目鏈接:http://poj.org/problem?id=1186
            /*
            題意:
                已知一個n元高次方程: 
            其中:x1, x2,,xn是未知數,k1,k2,,kn是系數,p1,p2,pn是指數。
            且方程中的所有數均為整數。 
            假設未知數1 <= xi <= M, i=1,,,n,求這個方程的整數解的個數。 
            1 <= n <= 6;1 <= M <= 150。   
            方程的整數解的個數小于2^31。 
            ★本題中,指數Pi(i=1,2,,n)均為正整數。 

            題解:
                DFS + HASH

            思路:
                一看到題目以為是解方程的,仔細一看,其實不然,題目有一個條件就是
            未知數的數目n很小,只有6,而且未知數也是有上限(1 <= M <= 150),所以
            當最大情況的時候可以枚舉3個未知數的值,一共150^3,然后插入到hash表中,
            再枚舉另外三個組成的所有情況,每次得到一個和為S時,只要查詢hash表中-S
            的數的數量,加到最后的答案即可。
                這樣就把本來150^6的復雜度開了個根號。枚舉數的時候因為n是不確定的
            ,所以我采用dfs來枚舉,這樣寫起來會方便許多。
            */


            #include 
            <iostream>

            using namespace std;

            #define maxn 4642307
            #define ll __int64
            int split;

            int nowPos[maxn], Case;
            int hash[maxn];
            int key[maxn];

            int n, M;

            int EXP(int a, int b) {
                
            if(b == 0)
                    
            return 1;
                
            int tmp = EXP(a*a, b/2);
                
            if(b & 1)
                    tmp 
            *= a;
                
            return tmp;
            }


            struct point {
                
            int k, p;
            }
            pt[6];

            void Insert(int val) {
                
            int s = val % maxn;
                
            if(s < 0)
                    s 
            += maxn;

                
            while(1{
                    
            if(nowPos[s] != Case) {
                        nowPos[s] 
            = Case;
                        hash[s] 
            = 1;
                        key[s] 
            = val;
                        
            return ;
                    }
            else {
                        
            if(key[s] == val) {
                            hash[s] 
            ++;
                            
            return ;
                        }

                        s 
            ++;
                        
            if(s == maxn)
                            s 
            = 0;
                    }

                }

            }


            int Query(int val) {
                
            int s = val % maxn;
                
            if(s < 0)
                    s 
            += maxn;
                
                
            while(1{
                    
            if(nowPos[s] != Case) {
                        
            return 0;
                    }
            else {
                        
            if(key[s] == val) {
                            
            return hash[s];
                        }

                        s 
            ++;
                        
            if(s == maxn)
                            s 
            = 0;
                    }

                }

            }


            void dfs0(int depth, int sum) {
                
            if(depth == split || depth == n) {
                    Insert(sum);
                    
            return ;
                }

                
            int i;
                
            for(i = 1; i <= M; i++{
                    dfs0(depth
            +1, sum + pt[depth].k * EXP(i, pt[depth].p) );
                }

            }


            void dfs1(int depth, int sum, int& ans) {
                
            if(depth == n) {
                    ans 
            += Query(- sum);
                    
            return ;
                }

                
            int i;
                
            for(i = 1; i <= M; i++{
                    dfs1(depth
            +1, sum + pt[depth].k * EXP(i, pt[depth].p), ans );
                }

            }


            int main() {
                
            int i;
                
            while(scanf("%d %d"&n, &M) != EOF) {

                    
            if(n <= 4{
                        split 
            = 2;
                    }
            else
                        split 
            = 3;

                    Case 
            ++;
                    
            for(i = 0; i < n; i++{
                        scanf(
            "%d %d"&pt[i].k, &pt[i].p);
                    }

                    dfs0(
            00);

                    
            int ans = 0;
                    
            if(n <= split) {
                        ans 
            = Query(0);
                    }
            else {
                        dfs1(split, 
            0, ans);
                    }

                    printf(
            "%d\n", ans);
                }

                
            return 0;
            }

            posted on 2011-04-09 23:55 英雄哪里出來 閱讀(1329) 評論(0)  編輯 收藏 引用 所屬分類: 數學

            波多野结衣久久一区二区| 国产一区二区精品久久凹凸| 日韩人妻无码一区二区三区久久99| 婷婷国产天堂久久综合五月| 人妻无码αv中文字幕久久琪琪布 人妻无码精品久久亚瑟影视 | 国产精品永久久久久久久久久| 久久国产视屏| 久久男人Av资源网站无码软件 | 97久久超碰成人精品网站| 精品乱码久久久久久夜夜嗨| 狠狠色婷婷久久综合频道日韩| 国产麻豆精品久久一二三| 婷婷久久综合九色综合绿巨人 | 理论片午午伦夜理片久久| 欧美黑人激情性久久| 怡红院日本一道日本久久| 99久久综合国产精品免费| 久久99精品久久久久久齐齐| 国产精品99精品久久免费| 一本色道久久88综合日韩精品| 99国内精品久久久久久久| 午夜精品久久久久久毛片| 色综合久久中文字幕综合网| 色综合久久中文色婷婷| 九九精品99久久久香蕉| 久久国产乱子伦免费精品| 久久成人小视频| 久久综合视频网| 久久综合亚洲色一区二区三区| 久久久中文字幕日本| 91亚洲国产成人久久精品| 色综合久久最新中文字幕| 色偷偷888欧美精品久久久| 久久福利青草精品资源站| 大伊人青草狠狠久久| 伊人久久免费视频| 亚洲乱亚洲乱淫久久| 欧美精品一区二区久久| 久久伊人五月丁香狠狠色| 久久久久久久久久久精品尤物| 97久久国产露脸精品国产|