青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

poj 2065 SETI

   題意比較糾結,搜索了把題意。
   給你一個素數P(P<=30000)和一串長為n的字符串str[]。字母'*'代表0,字母a-z分別代表1-26,這n個字符所代表的數字分別代表
f(1)、f(2)....f(n)。定義: f (k) = ∑0<=i<=n-1aiki (mod p) (1<=k<=n,0<=ai<P),求a0、a1.....an-1。題目保證肯定有唯一解。
   解題思路:高斯消元。根據上面的公式顯然可以列出有n個未知數的n個方程式:
   a0*1^0 + a1*1^1+a2*1^2+........+an-1*1^(n-1) = f(1)
   a0*2^0 + a1*2^1+a2*2^2+........+an-1*2^(n-1) = f(2)
   ..............
   a0*n^0 + a1*n^1+a2*n^2+........+an-1*n^(n-1) = f(n)
   然后采用高斯消元法來解上面的方程組即可。
   典型的高斯消元題,只是多了個modP,因此計算過程中可能需要擴展歐幾里德算法。

   說下所謂的高斯消元的思路,其實可以參看維基百科,
http://zh.wikipedia.org/wiki/%E9%AB%98%E6%96%AF%E6%B6%88%E5%8E%BB%E6%B3%95,大致過程是一直消變量。
比如剛開始,消第一個變量,消完之后只讓第一個方程含有第一個變量,然后消第二個變量,消完之后只讓第二個方程含第二個變量,以此
下去讓最后的方程含最后一個變量,而且最后一個方程中對于前N-1個變量的系數都是0,這樣就能解出這N個變量了。
   關于自由元指的是這個變量可以取任何值,得出這樣的結論是在消變量的過程中發現該變量的在第row個方程到第N方程中的系數都是0了,
所以可以取任何值。判斷無解的方式是,第row+1到第N個方程在高斯消元之后所有的系數必定是0,所以方程的值也必須是0。
   求方程的解得過程是從N個解開始逆推,第N-1個方程也就包含2個變量了,第N個變量和第N-1個變量,以此下去,就可以解出方程組了。
   具體的可以參照維基百科和代碼仔細分析。還有演算法筆記上也有高斯消元的解釋。

   代碼如下:
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
#define MAX (70 + 10)

int nMatrix[MAX][MAX];
int nAns[MAX];
void InitMatrix(char* szStr, int nN, int nP)
{
    memset(nMatrix, 0, sizeof(nMatrix));
    for (int i = 0; i < nN; ++i)
    {
        nMatrix[i][nN] = (szStr[i] == '*' ? 0 : szStr[i] - 'a' + 1);
    }
    for (int i = 0; i < nN; ++i)
    {
        int nTemp = 1;
        for (int j = 0; j < nN; ++j)
        {
            nMatrix[i][j] = nTemp;
            nTemp = (nTemp * (i + 1)) % nP;
        }
    }
}

int egcd(int nA, int nB, int& nX, int& nY)
{
    if (nA < nB)swap(nA, nB);
    if (nB == 0)
    {
        nX = 1, nY = 0;
        return nA;
    }
    int nRet = egcd(nB, nA % nB, nX, nY);
    int nT = nX;
    nX = nY;
    nY = nT - (nA / nB) * nY;
    return nRet;
}

int Gauss(int nN, int nP)
{
    int nR, nC;
    for (nR = nC = 0; nR < nN && nC < nN; ++nR, ++nC)
    {
        if (nMatrix[nR][nC] == 0)
        {
            for (int i = nR + 1; i < nN; ++i)
            {
                if (nMatrix[i][nC])
                {
                    for (int j = nC; j <= nN; ++j)
                    {
                        swap(nMatrix[nR][j], nMatrix[i][j]);
                    }
                    break;
                }
            }
        }

        if (nMatrix[nR][nC] == 0)
        {
            nR--;    //自由元
            continue;
        }
        int nA = nMatrix[nR][nC];
        for (int i = nR + 1; i < nN; ++i)
        {
            if (nMatrix[i][nC])
            {
                int nB = nMatrix[i][nC];
                for (int j = nC; j <= nN; ++j)
                {
                    nMatrix[i][j] = (nMatrix[i][j] * nA - nMatrix[nR][j] * nB) % nP;
                }
            }
        }
    }
    for (int i = nR; i < nN; ++i)
    {
        if (nMatrix[i][nN])
        {
            return -1;//無解
        }
    }
    
    int nX, nY;
    for (int i = nN - 1; i >= 0; i--)
    {
        int nSum = 0;
        for (int j = i + 1; j < nN; ++j)
        {
            nSum = (nSum + nMatrix[i][j] * nAns[j]) % nP;
        }
        
        nSum = (nMatrix[i][nN] - nSum + nP * nP) % nP;
        
        egcd(nP, (nMatrix[i][i] + nP) % nP, nX, nY);
        nY = (nY + nP) % nP;
        nAns[i] = (nY * nSum + nP) % nP;//第i個解
    }
    return 1 << (nN - nR);//返回解的個數,本題有唯一解
}

int main()
{
    int nT;

    scanf("%d", &nT);
    while (nT--)
    {
        int nP;
        int nN;
        char szStr[MAX];
        scanf("%d%s", &nP, szStr);
        nN = strlen(szStr);
        InitMatrix(szStr, nN, nP);
        Gauss(nN, nP);
        for (int i = 0; i < nN; ++i)
        {
            printf("%d%s", nAns[i], i == nN - 1 ? "\n" : " ");
        }
    }

    return 0;
}
   

posted on 2012-08-06 16:01 yx 閱讀(915) 評論(0)  編輯 收藏 引用 所屬分類: 數學題

<2012年4月>
25262728293031
1234567
891011121314
15161718192021
22232425262728
293012345

導航

統計

公告

常用鏈接

留言簿(3)

隨筆分類

隨筆檔案

me

好友

同學

網友

搜索

最新評論

閱讀排行榜

評論排行榜

青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            亚洲电影免费观看高清完整版在线| 亚洲欧美成人在线| 日韩视频第一页| 国产欧美日韩麻豆91| 欧美jizzhd精品欧美喷水 | 欧美亚日韩国产aⅴ精品中极品| 欧美有码视频| 国产欧美日韩一区| 国产精品视频网站| 欧美一区二区三区免费看| 亚洲第一狼人社区| 在线观看国产精品网站| 一本色道久久精品| 一区二区亚洲| 亚洲婷婷在线| 久久久天天操| 国产精品qvod| 国产精品嫩草久久久久| 免费观看成人鲁鲁鲁鲁鲁视频| 亚洲美女色禁图| 在线一区二区三区四区五区| 欧美人在线视频| 国产在线成人| 免费久久99精品国产自在现线| 免费亚洲婷婷| 国产精品扒开腿爽爽爽视频| 国产一区二区三区网站| 欧美有码视频| 久久国产精品久久精品国产 | 欧美欧美天天天天操| 国产精品扒开腿爽爽爽视频 | 亚洲高清不卡在线观看| 久久精品99国产精品日本| 亚洲日本久久| 亚洲国产精品成人精品| 欧美激情小视频| 欧美国产精品va在线观看| 欧美高清视频免费观看| 亚洲深夜激情| 欧美岛国激情| 欧美日韩视频第一区| 欧美电影在线观看完整版| 国产精品系列在线| 一区二区三区www| 亚洲欧美日韩专区| 亚洲经典三级| 亚洲精选久久| 一区二区三区欧美成人| 久久综合伊人77777蜜臀| 免费观看成人www动漫视频| 猛干欧美女孩| 先锋影音一区二区三区| 久久久国产精品亚洲一区 | 欧美不卡视频| 久久久精品久久久久| 国产亚洲精品一区二区| 欧美一区2区视频在线观看| 国产精品99久久久久久久女警 | 黄色一区二区三区| 久久精品免费观看| 麻豆久久婷婷| 欧美一区二区三区在线视频 | 欧美日韩一视频区二区| 亚洲美女av网站| 91久久夜色精品国产网站| 亚洲午夜视频在线观看| 久久黄金**| 国产精品免费一区豆花| 亚洲一级黄色片| 久久精品国产久精国产一老狼| 9久草视频在线视频精品| 久久成人精品无人区| 国产精品免费一区二区三区观看| 中文日韩在线视频| 久久精品一区二区| 亚洲精品国产拍免费91在线| 欧美成人69| 日韩午夜免费视频| 久久综合五月天婷婷伊人| 亚洲午夜高清视频| 欧美成人官网二区| 免费国产自线拍一欧美视频| 99精品免费| 亚洲午夜一二三区视频| 国产日韩精品一区二区三区在线| 激情视频一区二区| 欧美大片免费看| 欧美另类人妖| 午夜免费日韩视频| 久久久97精品| 亚洲一区高清| 亚洲伦理在线| 玖玖综合伊人| 激情一区二区| 亚洲欧洲日产国产网站| 国产精品国产三级国产aⅴ入口| 亚洲福利久久| 麻豆精品在线观看| 欧美国产视频日韩| 亚洲三级电影全部在线观看高清| 欧美专区亚洲专区| 免费成人网www| 在线观看日韩www视频免费| 亚洲风情亚aⅴ在线发布| 欧美视频一区二区| 亚洲一区激情| 久久久久久国产精品mv| 国产一区二区激情| 亚洲黄色在线视频| 国产美女精品视频免费观看| 午夜精彩视频在线观看不卡| 久久久国际精品| 亚洲一区图片| 免费的成人av| 久久久999国产| 欧美亚州在线观看| 亚洲国产你懂的| 在线不卡免费欧美| 亚洲欧美国产制服动漫| 99re6这里只有精品| 久久久精彩视频| 欧美在线免费视屏| 欧美无乱码久久久免费午夜一区| 欧美风情在线观看| 狠狠色狠狠色综合日日小说| 欧美激情久久久久| 嫩草国产精品入口| 欧美成人影音| 国产日韩欧美黄色| 亚洲综合不卡| 午夜精品久久久久久久久久久久| 欧美精品系列| 欧美一区二区黄| 欧美日韩日本视频| 91久久久久久久久久久久久| 亚洲黄色三级| 欧美高清不卡在线| 亚洲二区三区四区| 亚洲国产欧洲综合997久久| 久久国产精品亚洲va麻豆| 欧美中文字幕在线视频| 国产欧美在线| 亚洲在线一区| 欧美日韩一区二区精品| 亚洲乱码国产乱码精品精| 亚洲免费成人av| 欧美国产一区二区| 亚洲精品在线一区二区| 亚洲视频在线视频| 国产精品婷婷| 欧美在线观看网址综合| 亚洲日本欧美| 免费亚洲电影在线观看| 亚洲欧洲精品一区二区三区| 亚洲精选国产| 欧美性大战久久久久久久| 亚洲自拍另类| 免费成人av在线| 99精品国产福利在线观看免费| 欧美日韩亚洲不卡| 亚洲欧美一区二区三区久久 | 亚洲国产美女| 欧美精品一区二区精品网| 亚洲精品一区二区三区蜜桃久| 亚洲图片欧美日产| 国产一级揄自揄精品视频| 久久免费少妇高潮久久精品99| 亚洲午夜精品久久久久久浪潮| 欧美色精品天天在线观看视频 | 亚洲午夜免费视频| 国产亚洲欧美日韩一区二区| 久久影音先锋| 一片黄亚洲嫩模| 久久久久成人精品| 日韩视频在线你懂得| 国产精品视频一二三| 蜜月aⅴ免费一区二区三区 | 狠狠爱成人网| 欧美激情一区二区三区蜜桃视频 | 在线观看欧美成人| 欧美日韩国产色站一区二区三区| 欧美成人第一页| 亚洲欧美精品suv| 在线观看视频一区二区| 国产精品v日韩精品| 免费不卡在线观看av| 亚洲视频在线一区| 精品不卡在线| 国产精品久久久久久模特 | 欧美一区久久| 亚洲美女毛片| 女同性一区二区三区人了人一 | 国内偷自视频区视频综合| 欧美福利电影在线观看| 午夜天堂精品久久久久| 久久久久国产精品午夜一区| 国产精品卡一卡二卡三| 欧美成ee人免费视频| 午夜在线观看欧美| 日韩午夜av|