4

サブ配列には、正と負の両方の数値が含まれています。サブ配列の長さが k 以上になるような最大合計サブ配列を見つける必要があります。

Kadane のアルゴリズムを使用した C++ のコードを次に示します。

#include <iostream>

using namespace std;

int main(){

    int n,k;
    cin >> n >> k;
    int array[n];
    int sum = 0;
    int maxsum = 0;
    int beststarts[n];

    for(int i = 0;i < n; i++){
            cin >> array[i];
    }

    for(int i = 0;i < k-1;i ++){
            sum = sum+array[i];
            beststarts[i] = 0;
    }


    for(int i =  k-1;i < n; i++){ //best end search with min length;
            sum = sum+array[i];
            int testsum = sum;
            if(i > 0){
            beststarts[i] = beststarts[i-1];
            }
            for(int j = beststarts[i] ;i-j > k-1;j ++){
                    testsum = testsum - array[j];
                    if(testsum > sum){
                            beststarts[i] = j+1;
                            sum = testsum;
                    }
            }
            if(sum > maxsum){
                    maxsum = sum;
            }
    }

    cout << maxsum;

    return 0;
}

私のコードは正常に動作していますが、非常に遅く、コードを改善する方法が思い浮かびません。私はこの質問も読みました 合計がKで割り切れる最長の部分配列を見つけますが、それは私が望むものではありません。長さはkよりも大きくなる可能性があります。

4

2 に答える 2

3

この回答に基づく解決策

ライブデモ

#include <algorithm>
#include <iterator>
#include <iostream>
#include <numeric>
#include <ostream>
#include <utility>
#include <vector>

// __________________________________________________

template<typename RandomAccessIterator> typename std::iterator_traits<RandomAccessIterator>::value_type
max_subarr_k(RandomAccessIterator first,RandomAccessIterator last,int k)
{
    using namespace std;
    typedef typename iterator_traits<RandomAccessIterator>::value_type value_type;
    if(distance(first,last) < k)
        return value_type(0);
    RandomAccessIterator tail=first;
    first+=k;
    value_type window=accumulate(tail,first,value_type(0));
    value_type max_sum=window, current_sum=window;
    while(first!=last)
    {
        window += (*first)-(*tail) ;
        current_sum = max( current_sum+(*first), window );
        max_sum = max(max_sum,current_sum);
        ++first;
        ++tail;
    }
    return max_sum;
}

// __________________________________________________

template<typename E,int N>
E *end(E (&arr)[N])
{
    return arr+N;
}

int main()
{
    using namespace std;
    int arr[]={1,2,4,-5,-4,-3,2,1,5,6,-20,1,1,1,1,1};
    cout << max_subarr_k(arr,end(arr),4) << endl;
    cout << max_subarr_k(arr,end(arr),5) << endl;
}

出力は次のとおりです。

14
11
于 2012-11-02T05:33:18.690 に答える
0 
  int w(0);
    for (int i=0; i < k; i++) w += a[i];
    int run_sum(w), max_sum(w);
    for (int i=k; i < n; i++) {
              w = a[i] + max(w, w-a[i-k]); //  window will such that it will include run_sum
              run_sum = max(run_sum + a[i], w);
              max_sum = max(run_sum, max_sum); 
    }
    return max_sum;
于 2013-07-29T11:23:28.107 に答える