Problem:
ダーツ

問題

あなたは以下のルールでダーツゲームをすることになった.

あなたは,矢を的(まと)に向かって 4 本まで投げることができる.必ずしも 4 本全てを投げる必要はなく,1 本も投げなくてもかまわない.的は N 個の部分に区切られていて,各々の部分に点数 P1,..., PN が書か れている.矢が刺さった場所の点数の合計 S があなたの得点の基礎となる.S があらかじめ決められたある点数 M 以下の場合は S がそのままあなたの得点となる.しかし,S が M を超えた場合はあなたの得点は 0 点となる.

的に書かれている点数と M の値が与えられたとき,あなたが得ることのできる点数の最大値を求めるプログラムを作成せよ.

入力

入力は複数のデータセットからなる．各データセットは以下の形式で与えられる．

1 行目には,整数 N (1 ≤ N ≤ 1000) と M (1 ≤ M ≤ 200000000 = 2 × 108)がこの順に空白区切りで書かれている.2 行目以降の第 1 + i 行目 (1 ≤ i ≤ N ) には, Pi (1 ≤ Pi ≤ 100000000 = 108 ) が書かれている.

採点用データのうち, 配点の 20% 分は N ≤ 100 を満たし,配点の 50% 分はN ≤ 300 を満たす.

N, M がともに 0 のとき入力の終了を示す. データセットの数は 10 を超えない．

出力

データセットごとにあなたが得ることのできる点数の最大値を１行に出力する.

入出力例

入力例

4 50
3
14
15
9
3 21
16
11
2
0 0

出力例

48
20

1つ目の例では,15 点の部分に 3 本,3 点の部分に 1 本の矢が刺さった場合にあなたの得点は最大になり,その得点は 48 点である.

2つ目の例では,16 点の場所に 1 本,2 点の場所に 2 本の矢が刺さった場合にあなたの得点は最大になり,その得点は 20 点である.

上記問題文と自動審判に使われるデータは、情報オリンピック日本委員会が作成し公開している問題文と採点用テストデータです。

Solution in C (judge verdict: Time Limit Exceeded)

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
long p2[501501],c,all;
void add(long);
void put(long num,long place);
long search(long,long,long,long);
int cmp(const void*,const void*);
int main(){
	long p[1001],m,j,max,num,call,cp;
	int n,i;
	while(1){
		scanf("%d %ld",&n,&m);
		if(n==0) break;
		p2[0]=0;
		all=1;
		p[0]=0;
		for(i=1;i<=n;i++){
			scanf("%ld",&cp);
			p[i]=cp;
			qsort(p,i,sizeof(long),cmp);
			c=0;
			call=all;
			for(j=0;j<i;j++){
				if(p[j]+cp>m) break;
				add(p[j]+cp);
			}
		}
		max=0;
		for(j=1;j<all;j++){
			num=search(p2[j],m,0,all-1);
			max=(num>max)?num:max;
		}
		printf("%ld\n",max);
	}
	return 0;
}

void add(long p){
	for(;c<all;c++){
		if(p2[c]>p){
			put(p,c);
			all++;
			return;
		}else if(p2[c]==p){
			return;
		}
	}
	put(p,all);
	all++;
}

void put(long num,long place){
	long i;
	for(i=all;i>place;i--){
		p2[i]=p2[i-1];
	}
	p2[place]=num;
	for(i=0;i<=all;i++){
	}
}

long search(long num,long max,long start,long end){
	if(start==end) return num+p2[start];
	if(p2[(end-start+1)/2+start]>max-num){
		return search(num,max,start,(end-start+1)/2+start-1);
	}else{
		return search(num,max,(end-start+1)/2+start,end);
	}
}

int cmp(const void* a,const void* b){
	return *(long*)a-*(long*)b;
}