蓝书——AcWing 146. 序列

求n个不好求，我们可以先求2个序列，求出从小到大的前n个。

然后再与第三个序列进行求解，以此类推

这样问题就简化了

复杂度NM*logm  用优先队列维护下即可。

注意要防止   i  j+1         i+1 j，这两种情况都能递推到i+1,j+1.

我们先让i加，后让j加即可，熟悉的套路！  

 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
//typedef __int128 LL;
//typedef unsigned long long ull;
//#define F first
//#define S second
typedef long double ld;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
typedef pair<ld,ld> pdd;
const ld PI=acos(-1);
const ld eps=1e-9;
//unordered_map<int,int>mp;
#define ls (o<<1)
#define rs (o<<1|1)
#define pb push_back
//#define a(i,j) a[(i)*(m+2)+(j)]  //m是矩阵的列数
//pop_back()
const int seed=131;
const int M = 2e3+7;
/*int head[M],cnt;
struct EDGE{int to,nxt,val;}ee[M];
void add(int x,int y,int z){ee[++cnt].nxt=head[x],ee[cnt].to=y,ee[cnt].val=z,head[x]=cnt;}*/
int a[M],b[M],c[M];
struct node{
	int x,i,j,tg;
	bool operator <(const node r)const{
		return x>r.x;
	}
};
int n,m;
priority_queue<node>q;
void cal()
{
	while(!q.empty())q.pop();
//	cout<<a[1]<<"             "<<b[1]<<endl;
	q.push(node{a[1]+b[1],1,1,0});
	for(int i=1;i<=m;i++)
	{
		node tp=q.top();q.pop();
		if(tp.tg==0)q.push(node{a[tp.i+1]+b[tp.j],tp.i+1,tp.j,0});
		q.push(node{a[tp.i]+b[tp.j+1],tp.i,tp.j+1,1});
		c[i]=tp.x;
	//	cout<<tp.i<<" "<<tp.j<<" "<<tp.x<<"  -  "<<i<<endl;
	}
}
int main()
{
	ios::sync_with_stdio(false);
  	cin.tie(0);
	int t;
	cin>>t;
	while(t--)
	{	
		cin>>n>>m;
		for(int i=1;i<=m;i++)cin>>a[i];
		sort(a+1,a+1+m);
		for(int i=2;i<=n;i++)
		{
			for(int j=1;j<=m;j++)cin>>b[j];
			sort(b+1,b+1+m);
			cal();
			for(int j=1;j<=m;j++)
			a[j]=c[j];
			sort(a+1,a+1+m);
//			for(int k=1;k<=m;k++)
//			cout<<a[i]<<" ";
//			cout<<i<<endl;
		}
		for(int i=1;i<=m;i++)cout<<a[i]<<" ";
		cout<<endl;
	}
	return 0;
}