[讨论] 数组操作练习

qibbxxt

在cody上面看见这样一个问题，希望大家讨论下：
Given an array (x) of integers, the "counting" array (y) is showing the number of identical consecutive integers in x in front of the integer itself. For example, if

• x = 1
• y = [1 1],
  

because there is one "1". If then

• x = [1 1]
• y = [2 1],
  

because there are now two "1"s. Finally, a more complex example:

• x = [1 2 2 4 4 3 0 0 1]
• y = [1 1 2 2 2 4 1 3 2 0 1 1].
  

So y gets two elements for each series of identical integers in x. (I hope this problem does not exist already)

来自于：http://www.mathworks.cn/matlabcentral/cody/problems/1290-a-different-counting-method

lengyunfeng

抛砖引玉一下，供大家参考:

1. 
   
2. <P>function y=next_series(x)
   
3. if (isvector(x)~=1)                         %避开矩阵输入
   
4.     error(message('What you input is Not a vector!'));
   
5. else
   
6.     n=length(x);
   
7.     if (n==1)                               %当为标量时的处理办法
   
8.         y=[1 x];
   
9.     else
   
10.         [h,l]=size(x);                      %当为向量时的处理办法        
    
11.         if (h~=1)
    
12.             xtemp=x';
    
13.         else
    
14.             xtemp=x;
    
15.         end
    
16.         a=[ones(1,n);xtemp];
    
17.         [h,l]=size(a);
    
18.         p=1;
    
19.         for k=2:l
    
20.             if a(2,k)==a(2,p)
    
21.                 a(1,p)=a(1,p)+1;
    
22.                 a(1,k)=0;
    
23.             else
    
24.                 p=k;
    
25.             end
    
26.         end
    
27.         c=a(1,:);
    
28.         d=nonzeros(c);
    
29.         n=length(d);
    
30.         m=1;
    
31.         y=zeros(2,n);
    
32.         for k=1:l
    
33.             if a(1,k)~=0
    
34.                 y(:,m)=a(:,k);
    
35.                 m=m+1;
    
36.             end
    
37.         end
    
38.         y=y(:)';
    
39.     end
    
40. end           </P>
    
41. <P> </P>

复制代码

nwcwww

经典的run-length encoding问题:

1. function y = next_series(x)
   
2.   idx = find(diff([x,nan]));
   
3.   y0 = [diff([0 idx]); x(idx)];
   
4.   y = y0(:)';
   
5. end

复制代码

qibbxxt

感谢各位参与，下面我来大致列举一下cody上面的做法
1.利用regexprep的动态正则表达式功能来实现

1. function ans = next_series(x)
   
2. str2num(regexprep(char(x+48), '(\d)(\1*)', '${char(numel($2)+49)} $1 '));
   
3. end

复制代码

这种做法的优点在于代码简洁，缺点在于当数字不是0-9之间是，该方面失效，\d只能是0-9，但是此方法有助于练习正则表达式的使用，因为cody上面一个变量的size为2，所以大家尽量减少变量和函数的使用，这个时候就出现了ans和str2num之类，核心还是regexprep函数的使用

2. 使用2次diff

1. function ans = next_series(x)
   
2. find(diff([x inf]));
   
3. [diff([0 ans]); x(ans)];
   
4. ans(:)';
   
5. end

复制代码

diff函数在该类问题中有天然的优势，比如求解局部的极值等问题，改代码简洁高效，应用范围广，属于正解

3.  常规的做法

1. function y = next_series(x)
   
2. y=[];
   
3.   while ~isempty(x)
   
4.     val=find(x~=x(1),1);
   
5.     if isempty(val)
   
6.      val=length(x)+1;
   
7.     end
   
8. y=[y val-1 x(1)];
   
9. x=x(val:end);
   
10. end

复制代码

这个没有什么好说的，有很多不同的版本

4. 用了一次diff,用arrayfun函数操作的

1. function y = next_series(x)
   
2.   idx = find(diff(x));
   
3.   y0 = cell2mat(arrayfun(@(a, b) [ b-a+1; x(a) ],[ 1 idx + 1 ],[ idx length(x) ], 'uni', 0));
   
4.   y = y0(:)'
   
5. end

复制代码

时间有限，先不分析，这个代码，大家可以去分析分析

再次感谢大家参与，下次我还会再提供这样的一些题目，希望大家继续支持

bainhome

建议大家去cody做做题，练上100个题目估计MATLAB的理解都会上两三个层次。前一段时间做了几个正则替换表达式构造的题目，感觉就很有提高。

nwcwww

qibbxxt 发表于 2013-5-11 20:02
感谢各位参与，下面我来大致列举一下cody上面的做法
1.利用regexprep的动态正则表达式功能来实现这种做法的 ...

arrayfun的那个好像是我提交的。。后来觉得还是太丑就用RLE做了一遍