• <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
            數(shù)據(jù)加載中……

            PKU 1186 方程的解數(shù)

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

            題解:
                DFS + HASH

            思路:
                一看到題目以為是解方程的,仔細(xì)一看,其實(shí)不然,題目有一個(gè)條件就是
            未知數(shù)的數(shù)目n很小,只有6,而且未知數(shù)也是有上限(1 <= M <= 150),所以
            當(dāng)最大情況的時(shí)候可以枚舉3個(gè)未知數(shù)的值,一共150^3,然后插入到hash表中,
            再枚舉另外三個(gè)組成的所有情況,每次得到一個(gè)和為S時(shí),只要查詢hash表中-S
            的數(shù)的數(shù)量,加到最后的答案即可。
                這樣就把本來150^6的復(fù)雜度開了個(gè)根號。枚舉數(shù)的時(shí)候因?yà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 英雄哪里出來 閱讀(1317) 評論(0)  編輯 收藏 引用 所屬分類: 數(shù)學(xué)

            久久成人国产精品免费软件| 中文字幕人妻色偷偷久久| 久久久久久噜噜精品免费直播| 久久亚洲AV无码西西人体| 亚洲αv久久久噜噜噜噜噜| 久久久久中文字幕| 伊人久久精品影院| 久久九九全国免费| 亚洲综合伊人久久综合| 国产69精品久久久久99尤物| 亚洲AV无码久久精品蜜桃| 精品熟女少妇aⅴ免费久久| 久久夜色精品国产欧美乱| 久久免费视频6| 久久不射电影网| 久久人人爽人人爽人人片av高请| 久久久精品人妻无码专区不卡| 国内精品人妻无码久久久影院 | 久久久久亚洲AV无码观看| 99国产精品久久| 亚洲国产精品无码成人片久久| 亚洲精品国产综合久久一线| 国产999精品久久久久久| 91精品国产色综合久久| 色8久久人人97超碰香蕉987| 久久成人小视频| 久久精品国产黑森林| 99久久成人18免费网站| 精品无码久久久久国产| 久久人人妻人人爽人人爽| 亚洲午夜久久久影院| 无码任你躁久久久久久老妇App| 精品99久久aaa一级毛片| 久久综合九色综合精品| 久久人人爽人人爽人人片av高请| 色综合久久久久综合体桃花网 | 精品久久久久久久久中文字幕| 麻豆亚洲AV永久无码精品久久| 一级女性全黄久久生活片免费| 青青久久精品国产免费看| 亚洲国产精品嫩草影院久久|