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

oyjpArt ACM/ICPC算法程序設(shè)計空間

// I am new in programming, welcome to my blog
I am oyjpart(alpc12, 四城)
posts - 224, comments - 694, trackbacks - 0, articles - 6

Increasing Speed Limits

Problem

You were driving along a highway when you got caught by the road police for speeding. It turns out that they've been following you, and they were amazed by the fact that you were accelerating the whole time without using the brakes! And now you desperately need an excuse to explain that.

You've decided that it would be reasonable to say "all the speed limit signs I saw were in increasing order, that's why I've been accelerating". The police officer laughs in reply, and tells you all the signs that are placed along the segment of highway you drove, and says that's unlikely that you were so lucky just to see some part of these signs that were in increasing order.

Now you need to estimate that likelihood, or, in other words, find out how many different subsequences of the given sequence are strictly increasing. The empty subsequence does not count since that would imply you didn't look at any speed limits signs at all!

For example, (1, 2, 5) is an increasing subsequence of (1, 4, 2, 3, 5, 5), and we count it twice because there are two ways to select (1, 2, 5) from the list.

Input

The first line of input gives the number of cases, N. N test cases follow. The first line of each case contains n, m, X, Y and Z each separated by a space. n will be the length of the sequence of speed limits. m will be the length of the generating array A. The next m lines will contain the m elements of A, one integer per line (from A[0] to A[m-1]).

Using A, X, Y and Z, the following pseudocode will print the speed limit sequence in order. mod indicates the remainder operation.

for i = 0 to n-1
print A[i mod m]
A[i mod m] = (X * A[i mod m] + Y * (i + 1)) mod Z

Note: The way that the input is generated has nothing to do with the intended solution and exists solely to keep the size of the input files low.

Output

For each test case you should output one line containing "Case #T: S" (quotes for clarity) where T is the number of the test case and S is the number of non-empty increasing subsequences mod 1 000 000 007.

Limits

1 ≤ N ≤ 20
1 ≤ m ≤ 100
0 ≤ X ≤ 109
0 ≤ Y ≤ 109
1 ≤ Z ≤ 109
0 ≤ A[i] < Z

Small dataset

1 ≤ mn ≤ 1000

Large dataset

1 ≤ mn ≤ 500 000

Sample


Input
 

Output
 
2
5 5 0 0 5
1
2
1
2
3
6 2 2 1000000000 6
1
2

Case #1: 15
Case #2: 13

The sequence of speed limit signs for case 2 should be 1, 2, 0, 0, 0, 4.

沒趕上Round1A 郁悶。
Round1C Solve1和2,3的large不會做,菜。Rank好像是60多,能過。

賽后學(xué)習(xí)了下,也不算太難。
本來DP方程是這樣的
for(i = 0; i < n; ++i) {
 for(j = 0; j < i; ++j) {
  if(A[j] < A[i]) {
    dp[i] += dp[j];
  }
 }
}
如果對A排序并且離散化,則變成了
for(i=0; i < n; ++i) {
  for(j = 0; j < A[i]; ++j) {
   dp[A[i]] += dp[j];
  }
}


大家注意看,內(nèi)循環(huán)其實是一個區(qū)間求和。那么對于這種求和,線段樹只可以做到NlogN的。
記得以前寫過一道題的解題報告,是類似的。
pku1769 點樹解決塊查詢點操作

下面是代碼:(solve2函數(shù)是一個n^2的DP,偶水small input用的)
// Solution by alpc12  
#include 
<stdio.h>
#include 
<cassert>
#include 
<map>
#include 
<algorithm>
using namespace std;

const int M = 100;
const int N = 500010;
const int MOD = 1000000007;

typedef 
long long LL;

int n, m, X, Y, Z;
int A[N], S[N];
int st[1048576];
int upperbound = 524288;
int dp[N];

void generate() {
    
int i;
    
for(i = 0; i < n; ++i) {
        S[i] 
= A[i%m];
        A[i
%m] = ((LL)X*A[i%m]+(LL)Y*(i+1))%Z;
    }
    
for(i = 0; i < n; ++i) {
        A[i] 
= S[i];  
    }
}

int get(int x, int y) { // 左閉右開
    x += upperbound, y += upperbound;
    
int ans = 0;
    
while(x + 1 < y) {
        
if(x&1) { // x是右子樹 
            ans = (ans + st[x]) % MOD;
            x
++;
        }
        
if(y&1) { // y是右子樹
            y--;
            ans 
= (ans + st[y]) % MOD;
        }
        x 
>>= 1;
        y 
>>= 1;
    }
    
if(x < y) 
        ans 
= (ans + st[x]) % MOD;
    
return ans;
}

void ins(int x, int a) {
    x 
+= upperbound;
    
while(x > 0) {
        st[x] 
= (st[x] + a) % MOD;
        x 
>>= 1;
    }
}

void solve() {
    memset(st, 
0sizeof(st));
    sort(S, S 
+ n);
    map
<intint> mm;
    
int i, j = 0, ans = 0;
    
for(i = 0; i < n; ++i) {
        
if(!mm.count(S[i])) {
            mm[S[i]] 
= ++j;
        }
    }
    ins(
01);
    
for(i = 0; i < n; ++i) {
        A[i] 
= mm[A[i]];
        
int sum = get(0, A[i]);
        ans 
= (ans + sum) % MOD;
        ins(A[i], sum);
    }
    printf(
"%d\n", ans);
}

void solve2() {
    
int i, j, k;
    
for(i = 0; i < n; ++i) dp[i] = 1;
    
for(i = 1; i < n; ++i) {
        
for(j = 0; j < i; ++j) {
            
if(S[j] < S[i]) {
                dp[i] 
+= dp[j];
                dp[i] 
%= MOD;
            }
        }
    }
    LL sum 
= 0;
    
for(i = 0; i < n; ++i) {
        sum 
+= dp[i];
        sum 
%= MOD;
    }
    printf(
"%I64d\n", sum);
}

int main()
{
//    freopen("C-large.in", "r", stdin);
//    freopen("C-large.txt", "w", stdout);

    
int ntc, i, j, k, tc=0;
    scanf(
"%d"&ntc);
    
while(ntc--) {
        printf(
"Case #%d: "++tc);
        scanf(
"%d%d%d%d%d"&n, &m, &X, &Y, &Z);
        
for(i = 0; i < m; ++i) scanf("%d", A+i);
        generate();
//        solve2();
        solve();
    }
    
return 0;
}

青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            亚洲免费成人av| 亚洲理伦电影| 麻豆精品在线观看| 久久久久久综合| 久久亚洲精品一区| 欧美成人精品不卡视频在线观看| 久久一区二区三区国产精品| 久久精品国产一区二区三区免费看| 亚洲深夜激情| 久久精品国产综合| 欧美精品久久99久久在免费线| 欧美日韩精品久久久| 国产精品日韩精品欧美精品| 国产亚洲二区| 亚洲日本中文字幕区| 亚洲一二三区精品| 久久久天天操| 亚洲精品欧洲| 久久精品国产成人| 欧美日本网站| 精品999在线观看| 中日韩高清电影网| 久久综合网络一区二区| 日韩午夜激情电影| 久久久久久久波多野高潮日日| 欧美精品免费视频| 国产一区二区精品久久99| 亚洲精品人人| 久久亚洲国产精品日日av夜夜| 亚洲精品欧洲精品| 狂野欧美激情性xxxx欧美| 国产精品v欧美精品∨日韩| 在线观看国产成人av片| 亚洲女同精品视频| 亚洲国产高清在线观看视频| 一片黄亚洲嫩模| 久久影院午夜论| 国产精品自拍三区| 日韩视频在线一区| 老司机精品视频网站| 99国产精品国产精品毛片| 久久全球大尺度高清视频| 国产精品久久精品日日| 亚洲欧美在线一区二区| 午夜精品久久久久久久蜜桃app| 久久成人综合网| 一本色道精品久久一区二区三区| 狂野欧美激情性xxxx| 黄色国产精品| 久久久久看片| 亚洲欧美日韩综合国产aⅴ| 欧美日韩日日骚| 一本大道av伊人久久综合| 欧美激情中文字幕一区二区| 久久精品视频一| 狠狠色丁香婷婷综合| 欧美中文字幕| 午夜精品福利在线观看| 国产精品视频一二三| 西瓜成人精品人成网站| 中文国产一区| 国产精品久线观看视频| 亚洲综合不卡| 亚洲欧美日韩国产综合精品二区 | 在线视频免费在线观看一区二区| 欧美电影专区| 欧美二区在线观看| 亚洲乱码日产精品bd| 亚洲激情在线播放| 欧美日韩视频在线一区二区观看视频| 亚洲免费黄色| aa成人免费视频| 国产精品一级久久久| 欧美一站二站| 久久人人爽人人爽| 亚洲精品在线免费| 正在播放欧美一区| 国产日韩欧美中文| 免费视频亚洲| 欧美日韩精品福利| 久久狠狠亚洲综合| 久久久久在线观看| 在线亚洲一区| 欧美一级日韩一级| 亚洲国产高清高潮精品美女| 亚洲黄色av| 国产日韩欧美亚洲| 欧美肥婆在线| 国产精品久久久久久久免费软件| 久久aⅴ国产欧美74aaa| 麻豆精品精华液| 亚洲欧美日韩另类精品一区二区三区| 欧美一区成人| 妖精视频成人观看www| 午夜精品福利视频| 日韩午夜中文字幕| 欧美一级在线亚洲天堂| 日韩视频在线一区二区| 欧美一区2区视频在线观看| 亚洲精品欧洲| 久久久久一区二区| 欧美一区二区性| 欧美激情导航| 免费视频亚洲| 宅男噜噜噜66一区二区66| 亚洲麻豆视频| 红桃视频一区| 亚洲视频在线二区| 亚洲欧洲精品一区二区三区不卡| 亚洲一级二级在线| 亚洲精品国产日韩| 欧美在线视频一区二区三区| av成人免费观看| 久久这里只有| 久久日韩粉嫩一区二区三区| 国产精品久久久久av免费| 亚洲第一黄色| 狠狠色狠狠色综合人人| 亚洲四色影视在线观看| 9久草视频在线视频精品| 久久久一本精品99久久精品66| 亚洲在线观看视频网站| 亚洲经典三级| 狼狼综合久久久久综合网| 欧美一级免费视频| 国产精品黄视频| 一区二区三区久久网| 99av国产精品欲麻豆| 欧美成人按摩| 亚洲国产一二三| 亚洲久色影视| 欧美日韩高清在线| 亚洲人成久久| 日韩午夜电影av| 欧美激情第五页| 亚洲精选一区二区| 亚洲深夜影院| 国产精品都在这里| 亚洲综合999| 欧美影院在线| 韩国福利一区| 猛干欧美女孩| 亚洲日本中文| 午夜亚洲激情| 国内成人精品视频| 久久综合九九| 最新亚洲一区| 亚洲影视在线| 国产一区二区三区高清| 久久久久国产精品一区| 欧美高清在线一区二区| 亚洲日本成人女熟在线观看| 欧美精品国产精品| 亚洲图片欧美一区| 久久久久免费视频| 亚洲高清一区二区三区| 欧美久久综合| 亚洲欧美一区二区精品久久久| 欧美中文字幕视频在线观看| 韩国av一区二区三区在线观看| 久久久免费精品| 亚洲毛片一区| 久久亚洲私人国产精品va| 亚洲国产中文字幕在线观看| 欧美日韩1区2区3区| 午夜精品久久久久影视 | 男女av一区三区二区色多| 亚洲第一中文字幕| 亚洲欧美日韩另类精品一区二区三区| 国产最新精品精品你懂的| 国产欧美一区二区白浆黑人| 亚洲国产精品一区| 亚洲自拍偷拍色片视频| 黄色小说综合网站| 欧美国产欧美亚洲国产日韩mv天天看完整 | 欧美在线网址| 亚洲免费观看在线观看| 国产日韩在线看| 欧美不卡视频| 午夜宅男欧美| 亚洲精品在线视频观看| 久久久青草婷婷精品综合日韩| 99精品热视频| 国内精品久久久久久久影视麻豆| 欧美黑人一区二区三区| 欧美在线网站| 亚洲午夜电影在线观看| 亚洲国产经典视频| 久久久国产成人精品| 亚洲图片欧美午夜| 亚洲日本aⅴ片在线观看香蕉| 国产欧美在线观看| 国产精品成人一区二区艾草| 农村妇女精品| 久久久久久久久久久久久久一区| 亚洲综合日本| 亚洲午夜精品久久久久久app| 亚洲精品视频一区二区三区| 久热精品视频在线免费观看| 久久精品91久久香蕉加勒比|