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

糯米

TI DaVinci, gstreamer, ffmpeg
隨筆 - 167, 文章 - 0, 評論 - 47, 引用 - 0
數據加載中……

POJ 3167 Cow Patterns 牛題

此題只有130人solved!也算小牛題了,剛開始看到的時候,打算不做了的,但后來想了一下,發現有一點點思路。
經過好幾個小時的奮戰,居然做出來了!那感覺非常爽!

思路:

首先此題的變態之處是,要求比較的是排名,就不是單純的字符串匹配了。
如果每次都重新求排名然后跟pattern比較,復雜度 O(NK),八成會超時。
關鍵是:
一,不能每次都重新求排名
二,不能逐個逐個的和pattern做比較
看來也只有動態規劃或者hash能滿足這種需求。

對于前 K 個數字,求出一個初始的 hash 值。
對于后 N - K 個數字,不斷的把大小為 K 的窗口向后一格格滑動,比如 [1, 5] -> [2, 6] -> [3, 7]。
hash 值隨之更新,如果和pattern的hash值一致,那么當前的區間就是一個答案了。

一開始,想到的hash函數是,
hash = rank[1] * a[1] + rank[2] * a[2] + ... rank[K] * a[K]
其中:
所有數字的類型都為32bit無符號整數
數組 a 為等差數列
rank 為著窗口內的排名,跟 pattern 的格式一樣,長度為 K。

由于 S <= 25 ,所以可以把數值相同的元素保存在一起,因它們的rank值也是一樣的,
那 hash 函數就變成了這樣子:
hash = 1 * (a[1] + a[3] ...) + 2 * (a[2] + a[7] + ...) + ...
其中 rank[1] = rank[3] = 1,rank[2] = rank[7] = 2。。。
可以維護大小為 S 的數組 pos 、rank_val:
pos[i] = { 元素 i 出現的位置對應的 hash 值,如果沒出現過,則為 0 }
pos[1] = a[1] + a[3] + ...   
pos[2] = a[2] + a[7] + ...
。。。
rank_val[i] = { 元素 i 對應的 rank }
rank_avl[1] = 1
rank_avl[2] = 2
。。。
因而:
hash = rank_val[1] * pos[1] + rank_val[2]  * pos[2] + ...  + rank_val[K] * pos[K]

窗口每次滑動一格的時候,都要先增加一個元素,然后減掉一個元素。
再維護一個大小為 S 的數組 cnt:
cnt[i] = { 元素 i 在窗口中出現的次數 }
這樣,每次增刪元素的時候,就很容易得出 rank_val 數組的內容了,就跟基數排序一個原理。

增加元素 i 的時候,首先增加 pos[i], cnt[i] 的值,然后掃描一遍 cnt 更新 rank_val 的值。
刪除元素的情況類似。

窗口向后移動了一格之后:
pos[2] = a[2] + a[7]
應該變化為
pos[2] = a[1] + a[6]
這時候,由于 a 是等差數列,假設它的公差為 d,那么只需要 pos[2] -= cnt[2] * d 就能在 O(1) 時間內完成這種“移動”操作了!

因此,在窗口移動一格的過程中,步驟為:
一,增加窗口末尾后的第一個元素
二,刪除窗口開頭元素
三,“移動”窗口中所有元素
四,與 pattern 比較 hash 值

總的時間復雜度是 O(NS)。

后來,發現有的數據過不了,hash 沖突了,又將 hash 函數變化了一下:
hash = rank_val[1]^5 * pos[1] + rank_val[1]^5 * pos[2] ... rank_val[K]^5 * pos[K]
發現求到 4 次方的時候貌似還是有沖突。。。
到 5 次方的時候就 AC 了!

代碼跑了 350 ms

代碼:
#include <stdio.h>

#define MAX_K 25032
#define MAX_N 100032
#define MAX_S 32
#define HASH_START 0xdead
#define HASH_INC 0xbeef

struct spot_node {
    unsigned 
int cnt, rank, sum;
}
;
struct seq_node {
    
struct spot_node spot[MAX_S];
    unsigned 
int hash;
}
;
struct seq_node ptn, seq;
unsigned 
int N, K, S;
unsigned 
int hash_idx[MAX_K], in[MAX_N + MAX_K], rank_pow[MAX_S];
unsigned 
int ans[MAX_N], ans_cnt;

#define always_inline __inline 

always_inline unsigned 
int rank_hash(unsigned int i)
{
    
return rank_pow[i];
}


always_inline 
void seq_insert(struct seq_node *s, unsigned int val)
{
    
struct spot_node *t;

    t 
= &s->spot[val];
    t
->cnt++;
    t
->sum += hash_idx[K + 1];
    s
->hash += hash_idx[K + 1* rank_hash(t->rank);
    
if (t->cnt != 1
        
return ;
    
for (val++; val <= S; val++{
        t 
= &s->spot[val];
        s
->hash += t->sum * (rank_hash(t->rank + 1- rank_hash(t->rank));
        t
->rank++;
    }

}


always_inline 
void seq_delete(struct seq_node *s, unsigned int val)
{
    
struct spot_node *t;

    t 
= &s->spot[val];
    t
->cnt--;
    t
->sum -= hash_idx[1];
    s
->hash -= hash_idx[1* rank_hash(t->rank);
    
if (t->cnt != 0
        
return ;
    
for (val++; val <= S; val++{
        t 
= &s->spot[val];
        s
->hash -= t->sum * (rank_hash(t->rank) - rank_hash(t->rank - 1));
        t
->rank--;
    }

}


always_inline 
void seq_advance(struct seq_node *s)
{
    unsigned 
int i;
    
struct spot_node *t;

    
for (i = 1; i <= S; i++{
        t 
= &s->spot[i];
        s
->hash -= t->cnt * rank_hash(t->rank) * HASH_INC;
        t
->sum -= t->cnt * HASH_INC;
    }

}


always_inline 
void seq_init(struct seq_node *s, unsigned int *arr)
{
    unsigned 
int i, r;
    
struct spot_node *t;

    s
->hash = 0;
    
for (i = 1; i <= K; i++{
        t 
= &s->spot[*arr++];
        t
->sum += hash_idx[i];
        t
->cnt++;
    }

    
for (i = r = 1; i <= S; i++{
        t 
= &s->spot[i];
        t
->rank = r;
        
if (t->cnt) {
            s
->hash += rank_hash(t->rank) * t->sum;
            r
++;
        }

    }

}


int main()
{
    unsigned 
int i;

    freopen(
"e:\\test\\in.txt""r", stdin);

    scanf(
"%d%d%d"&N, &K, &S);
    
for (i = 0; i < N + K; i++)
        scanf(
"%d"&in[i]);

    hash_idx[
1= HASH_START;
    
for (i = 2; i <= K + 1; i++)
        hash_idx[i] 
= hash_idx[i - 1+ HASH_INC;
    
for (i = 1; i <= S; i++)
        rank_pow[i] 
= i * i * i * i * i;

    seq_init(
&ptn, in + N);
    seq_init(
&seq, in);
    
for (i = K; ; i++{
        
if (seq.hash == ptn.hash) 
            ans[ans_cnt
++= i - K + 1;
        
if (i == N)
            
break;
        seq_insert(
&seq, in[i]);
        seq_delete(
&seq, in[i - K]);
        seq_advance(
&seq);
    }


    printf(
"%d\n", ans_cnt);
    
for (i = 0; i < ans_cnt; i++)
        printf(
"%d\n", ans[i]);

    
return 0;
}



posted on 2010-03-18 15:30 糯米 閱讀(1214) 評論(2)  編輯 收藏 引用 所屬分類: POJ

評論

# re: POJ 3167 Cow Patterns 牛題  回復  更多評論   

這個是用KMP做..
2010-04-05 17:24 | cnx

# re: POJ 3167 Cow Patterns 牛題[未登錄]  回復  更多評論   

@cnx
Disscuss里面有人說過用KMP做,但我就是想不通怎么用做。。
能否給個思路?
2010-04-05 17:26 | 糯米
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            欧美国内亚洲| 国产欧美欧美| 性欧美在线看片a免费观看| 亚洲大片精品永久免费| 欧美日韩123| 国精产品99永久一区一区| 国产一区二区黄色| 99视频一区二区| 久久精品一区二区国产| 欧美成人精品一区二区三区| 欧美日韩精品一区| 国产亚洲福利| 亚洲免费在线播放| 欧美国产精品v| 小黄鸭精品密入口导航| 美女视频网站黄色亚洲| 国产精品久久久久久久久婷婷| 国产三级欧美三级日产三级99| 久久综合久久综合九色| 国产精品久久久久91| 亚洲电影免费在线| 久久久中精品2020中文| 亚洲一区中文| 在线中文字幕一区| 免费黄网站欧美| 亚洲影音一区| 国产欧美日韩精品a在线观看| 午夜日韩激情| 久久综合久久综合久久综合| 亚洲欧美变态国产另类| 国产精品超碰97尤物18| 亚洲网站在线| 久久久久久国产精品一区| 一区二区欧美精品| 久久精品国产清高在天天线| 99亚洲视频| 亚洲综合清纯丝袜自拍| 精品不卡在线| 亚洲图片欧洲图片av| 亚洲精品免费网站| 欧美一区二区三区视频在线观看| 亚洲欧洲日韩在线| 久久久久久久久久久一区| 欧美一区二区三区四区在线观看| 久久这里有精品15一区二区三区 | 免费看精品久久片| 欧美日韩中文字幕日韩欧美| 另类av导航| 亚洲国产高清自拍| 麻豆成人在线观看| 亚洲夫妻自拍| 亚洲国产专区| 欧美理论电影在线播放| 99精品免费网| 久久久久久久一区二区| 国产一区二区精品久久91| 欧美综合第一页| 国产精品www网站| 亚洲一区二区在线| 蜜臀va亚洲va欧美va天堂| 国产一区二区三区四区老人| 欧美一区二区三区精品| 欧美经典一区二区| 99国产精品视频免费观看| 亚洲一区二区欧美日韩| 国产精品色网| 免费不卡亚洲欧美| 亚洲一区网站| 欧美成人精品激情在线观看| 亚洲最新在线视频| 国产模特精品视频久久久久| 久久动漫亚洲| 亚洲夜间福利| 亚洲第一页在线| 久久久久久香蕉网| 亚洲综合成人婷婷小说| 国产精品自拍小视频| 99视频精品全国免费| 欧美v国产在线一区二区三区| 亚洲一区二区三区高清不卡| 影音先锋日韩精品| 国产精品剧情在线亚洲| 欧美日韩国产限制| 玖玖国产精品视频| 久久精品电影| 久久精品亚洲一区二区三区浴池 | 亚洲毛片网站| 亚洲国产aⅴ天堂久久| 欧美成人自拍| 亚洲人精品午夜| 久久亚洲风情| 午夜精品国产更新| 一区二区三区日韩精品| 亚洲一区二区在线免费观看| 欧美国产一区二区在线观看 | 欧美激情成人在线视频| 欧美激情精品久久久久久久变态 | 99精品久久久| 一本一本久久a久久精品综合妖精 一本一本久久a久久精品综合麻豆 | 一本色道久久综合狠狠躁的推荐| 亚洲福利免费| 国产精品99久久久久久久女警| 亚洲午夜国产一区99re久久| 午夜精品久久久久久久蜜桃app| 亚洲视频在线观看三级| 久久久999精品免费| 久久综合一区二区三区| 亚洲精品日本| 久久久久国产精品午夜一区| 欧美精品黄色| 国产日韩欧美亚洲一区| 亚洲人成久久| 久久精品欧洲| 亚洲无线视频| 欧美大片免费观看| 国产亚洲欧美色| 亚洲一区综合| 欧美高清影院| 噜噜噜91成人网| 黄色av一区| 久久网站热最新地址| 亚洲国内在线| 一区二区激情视频| 看片网站欧美日韩| 欧美专区在线观看| 国产无一区二区| 午夜免费日韩视频| 日韩亚洲欧美中文三级| 欧美xxx在线观看| 在线精品亚洲一区二区| 久久综合图片| 久久亚洲一区二区三区四区| 欧美va亚洲va香蕉在线| 最近中文字幕日韩精品| 亚洲国产mv| 欧美日韩在线播放三区| 亚洲一区二区免费视频| 亚洲一区免费在线观看| 国产一区亚洲一区| 开心色5月久久精品| 农夫在线精品视频免费观看| 亚洲国产精品成人综合| 99国产精品久久久久久久久久 | 国产精品成人在线观看| 亚洲欧美国产制服动漫| 欧美亚洲一区| 亚洲国产欧美一区二区三区同亚洲 | 一区二区高清视频在线观看| 午夜精品久久久久久久蜜桃app| avtt综合网| 国产精品久久久久久av下载红粉| 一区二区不卡在线视频 午夜欧美不卡'| 欧美在线视频免费| 亚洲一区中文| 99re视频这里只有精品| 亚洲一区二区三区高清| 亚洲黄色高清| 欧美中文字幕在线视频| 欧美精品日韩www.p站| 香蕉久久夜色| 欧美日韩福利视频| 老司机免费视频一区二区| 国产精品国产三级国产aⅴ入口 | 亚洲免费观看高清在线观看| 亚洲一区欧美激情| 亚洲男女自偷自拍图片另类| 老司机67194精品线观看| 久久久精品动漫| 国产主播精品| 久久久91精品国产一区二区三区 | 国产精品国产亚洲精品看不卡15| 亚洲欧洲在线免费| 一级成人国产| 久久久夜精品| 妖精成人www高清在线观看| 亚洲天堂成人| 在线观看国产欧美| 欧美日韩成人在线视频| 欧美影院成年免费版| 亚洲国产精品第一区二区| 亚洲自拍偷拍色片视频| 在线观看欧美一区| 国产精品自拍小视频| 免费日韩av片| 欧美**字幕| 亚洲激情不卡| 欧美激情综合五月色丁香小说| 亚洲精品九九| 亚洲人成绝费网站色www| 在线成人激情黄色| 国产精品国产三级国产普通话蜜臀 | 校园激情久久| 这里只有精品视频在线| 你懂的视频欧美| 久久午夜羞羞影院免费观看| 欧美在现视频| 一本色道久久88综合日韩精品| 久久精品视频在线| 亚洲欧美日韩精品在线| 亚洲砖区区免费|