Znzryb's Blog

2025“钉耙编程”中国大学生算法设计暑期联赛（1）——树上LCM

Posted on 2025-07-20 Edited on 2026-02-17

思路讲解

三个集合的容斥原理

那么不难发现规律，就是奇加偶减。这就是这段代码的由来。

1
2
3

ll cn=bi.count();
if(cn%2==1) del+=(sz+1)*sz/2;
else del-=(sz+1)*(sz)/2;

那么直接计算合法路径数量比较难，但是我们可以枚举不合法子集。如果s是101（二进制），那么就是第一个素因子要不符合要求（LCM只关注这些数当中的最大的，你可以理解为对素因子取了一个并集，那么gcd就是取了一个交集）。

// s如果这一位是1，就代表这一位的质因数必须不满足要求
	FOR(s,0,(1<<siz)-1){
		vb vis(N+5,0);
		bitset<11> bi(s);
		// 检查该数字（该点）是否满足要求
		auto check=[&](ll u)->bool{
			FOR(i,0,siz-1){
				if(!bi[i]) continue;
				// 该数不满足要求，因为我们要求这一位质因数数量缺失，或者太多
				if(fach[u][i]==cnt[i]){
					return false;
				}
			}
			return true;
		};
		FOR(u,1,N){
			if(vis[u]) continue;
			if(!valid[u]) continue;
			vis[u]=true;
			if(!check(u)) continue;
			queue<ll> q;q.push(u);
			ll sz=0;	// 维护联通块内的点的数量
			while(!q.empty()){
				ll u=q.front();q.pop();
				vis[u]=true;
				++sz;
				for(auto v:g[u]){
					if(vis[v]) continue;
					if(!valid[v]) continue;
					if(!check(v)) continue;
					q.push(v);
				}
			}
			// sz*(sz-1)/2+sz
			ll cn=bi.count();
			if(cn%2==1) del+=(sz+1)*sz/2;
			else del-=(sz+1)*(sz)/2;
			// ans+=(sz+1)*sz/2;
		}
#ifdef LOCAL
    debug(del);
#endif
	}
	ll ans=-del;
	cout<<ans<<"\n";
}

AC代码

https://acm.hdu.edu.cn/contest/view-code?cid=1172&rid=22776

源代码

// Problem: 树上LCM
// Contest: HDOJ
// URL: https://acm.hdu.edu.cn/contest/problem?cid=1172&pid=1007
// Memory Limit: 524288 KB
// Time Limit: 2000 ms
// by gospel_rock
// 
// Powered by CP Editor (https://cpeditor.org)

#include <bits/stdc++.h>
#define all(vec) vec.begin(),vec.end()
#define pb push_back
#define EB emplace_back
#define EM emplace
#define FI first
#define SE second
#define SZ(a) ((int) a.size())
#define FOR(i, a, b) for (long long i = (a); i <= (b); ++i)
#define ROF(i, a, b) for (long long i = (a); i >= (b); --i)
#define debug(var) cerr << #var <<":"<<var<<"\n";
#define lson(var) (var<<1)
#define rson(var) ((var<<1)+1)

using namespace std;

using  ll = long long;using ull = unsigned long long;
using DB=double;using LD=long double;
// using i128 = __int128;
using pdd = pair<DB,DB>;using plb = pair<ll,bool>;
using pll = pair<ll,ll>;using vll=vector<ll>;using vb=vector<bool>;
using arr3 = array<ll,3> ;using arr2 = array<ll,2>;
using tll = tuple<ll,ll,ll>;
using vpll = vector<pll>;
constexpr ll MAXN=(ll)3e5+10,INF=(ll)1e17+9; // 1e18+9也是素数
constexpr ll MAXV=(ll)1e6+10;
constexpr ll MAXM=(ll)2e6+10;
constexpr ll mod=(ll)1e9+7; // 998244353
constexpr double eps=1e-8;const double pi=acos(-1.0);
ll N,M,Q,X,K,T,lT,A[MAXN];
/*
2*3*5*7*11*13*17*19*23-1e7
213092870.0
所以说不可能超过9位（素因数）

思路总结来说就是容斥原理
*/
//  minp[i] 存储的是i的最小素因数
vector<ll> minp,primes;
void sieve(ll n){
	minp.assign(n+5,0);
	for(int i=2;i<=n;++i){
		if(minp[i]==0){		// 是素数
			minp[i]=i;
			primes.push_back(i);
		}
		for(auto p:primes){
			if(i*p>n){	// 超范围了
				break;
			}
			minp[i*p]=p;
			// 为了保证每个合数都只被其【最小的】那个质因数标记一次 break
			// i * p_next = (p * k) * p_next 最小的应该是p
			if(p==minp[i]){
				break;
			}
		}
	}
}
inline void __(){
	cin>>N>>X;
	vector<vector<ll> > g(N+5);
	FOR(i,1,N-1){
		ll u,v;
		cin>>u>>v;
		g[u].pb(v);
		g[v].pb(u);
	}
	vb valid(N+5,0);
	FOR(i,1,N){
		cin>>A[i];
		if(X%A[i]==0) valid[i]=true;
		// 有不能整除节点的路径必然不能形成lcm=x
	}
	vector<ll> fac;ll siz=0;vll cnt;
	auto primeDe=[&](ll x)->void{
		while(x>1){
			fac.push_back(minp[x]);
			x/=minp[x];
		}
		sort(all(fac));
		vll facu(all(fac));fac.resize(unique(all(fac))-fac.begin());
		siz=SZ(fac);
		for(auto &f:fac){
			ll c=upper_bound(all(facu),f)-lower_bound(all(facu),f);
			cnt.EB(c);
		}
	};
	primeDe(X);	// 返回去重的fac，以及因数计数cnt
	
	vector<vector<ll> > fach(N+5,vector<ll>(10,0));	
	FOR(i,1,N){
		if(!valid[i]) continue;
		ll a=A[i];
		while(a>1){
			ll fa=minp[a];
			a/=minp[a];
			ll idx=lower_bound(all(fac),fa)-fac.begin();
			fach[i][idx]++;
		}
	}
#ifdef LOCAL
    FOR(i,0,SZ(fac)-1){
    	debug(fac[i]);
    	debug(cnt[i]);
    }
#endif
	ll del=0;
	// s如果这一位是1，就代表这一位的质因数必须不满足要求
	FOR(s,0,(1<<siz)-1){
		vb vis(N+5,0);
		bitset<11> bi(s);
		// 检查该数字（该点）是否满足要求
		auto check=[&](ll u)->bool{
			FOR(i,0,siz-1){
				if(!bi[i]) continue;
				// 该数不满足要求，因为我们要求这一位质因数数量缺失，或者太多
				if(fach[u][i]==cnt[i]){
					return false;
				}
			}
			return true;
		};
		FOR(u,1,N){
			if(vis[u]) continue;
			if(!valid[u]) continue;
			vis[u]=true;
			if(!check(u)) continue;
			queue<ll> q;q.push(u);
			ll sz=0;	// 维护联通块内的点的数量
			while(!q.empty()){
				ll u=q.front();q.pop();
				vis[u]=true;
				++sz;
				for(auto v:g[u]){
					if(vis[v]) continue;
					if(!valid[v]) continue;
					if(!check(v)) continue;
					q.push(v);
				}
			}
			// sz*(sz-1)/2+sz
			ll cn=bi.count();
			if(cn%2==1) del+=(sz+1)*sz/2;
			else del-=(sz+1)*(sz)/2;
			// ans+=(sz+1)*sz/2;
		}
#ifdef LOCAL
    debug(del);
#endif
	}
	ll ans=-del;
	cout<<ans<<"\n";
}

int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);cout.tie(0);
	sieve((ll)1E7+10);
	cin>>lT;
	while(lT--)
		__();
	return 0;
}
/*
AC
https://acm.hdu.edu.cn/contest/view-code?cid=1172&rid=22776
*/

心路历程（WA，TLE，MLE……）

CF-1038-C. Manhattan Pairs（选择不重复的对，使得这个曼哈顿距离之和最大）

Posted on 2025-07-19 Edited on 2026-02-17

思路讲解

那么容易想到，按照x排序，可以让sumx最大的方式就是排序，然后左边匹配右边。

而且，这样子匹配，左边的谁匹配右边的谁也是不重要的，至少对于sumx不重要，因此我们就可以按照使sumy最大的方式进行匹配。

同理，计算按照y排列能得到的答案即可。

AC代码

https://codeforces.com/contest/2122/submission/329811774

源代码

#include <bits/stdc++.h>

using namespace std;

using ll = long long;

// 定义点结构，包含x坐标、y坐标和索引
struct Point {
    ll x, y;
    int id;
};

// 主函数，用于处理每个测试用例
inline void processTestCase() {
    int n;
    cin >> n;  // 读取点的数量n（偶数）

    vector<Point> points(n);  // 创建存储点的向量
    for (int i = 0; i < n; ++i) {
        cin >> points[i].x >> points[i].y;  // 读取每个点的坐标
        points[i].id = i + 1;  // 索引从1开始
    }

    // 定义按x坐标升序排序的lambda函数，如果x相同则按id升序
    auto compareByX = [](const Point& a, const Point& b) {
        return a.x < b.x || (a.x == b.x && a.id < b.id);
    };

    // 定义按y坐标升序排序的lambda函数，如果y相同则按id升序
    auto compareByY = [](const Point& a, const Point& b) {
        return a.y < b.y || (a.y == b.y && a.id < b.id);
    };

    // 定义按y坐标降序排序的lambda函数，如果y相同则按id升序
    auto compareByYDesc = [](const Point& a, const Point& b) {
        return a.y > b.y || (a.y == b.y && a.id < b.id);
    };

    // 定义按x坐标降序排序的lambda函数，如果x相同则按id升序
    auto compareByXDesc = [](const Point& a, const Point& b) {
        return a.x > b.x || (a.x == b.x && a.id < b.id);
    };

    vector<pair<int, int>> pairsOption1, pairsOption2;  // 存储两种方案的配对
    ll sumOption1 = 0, sumOption2 = 0;  // 存储两种方案的总曼哈顿距离

    // 方案1: 按x坐标分割成左右两组
    vector<Point> pointsSortedByX = points;  // 复制点集
    sort(pointsSortedByX.begin(), pointsSortedByX.end(), compareByX);  // 按x升序排序

    vector<Point> leftGroup(pointsSortedByX.begin(), pointsSortedByX.begin() + n / 2);  // 左半部分（小x）
    vector<Point> rightGroup(pointsSortedByX.begin() + n / 2, pointsSortedByX.end());  // 右半部分（大x）

    sort(leftGroup.begin(), leftGroup.end(), compareByY);  // 左组按y升序排序
    sort(rightGroup.begin(), rightGroup.end(), compareByYDesc);  // 右组按y降序排序

    for (int i = 0; i < n / 2; ++i) {
        ll dx = abs(leftGroup[i].x - rightGroup[i].x);  // 计算x差的绝对值
        ll dy = abs(leftGroup[i].y - rightGroup[i].y);  // 计算y差的绝对值
        sumOption1 += dx + dy;  // 累加到总和
        int a = leftGroup[i].id, b = rightGroup[i].id;  // 获取索引
        if (a > b) swap(a, b);  // 确保a < b
        pairsOption1.push_back({a, b});  // 添加配对
    }

    // 方案2: 按y坐标分割成上下两组
    vector<Point> pointsSortedByY = points;  // 复制点集
    sort(pointsSortedByY.begin(), pointsSortedByY.end(), compareByY);  // 按y升序排序

    vector<Point> bottomGroup(pointsSortedByY.begin(), pointsSortedByY.begin() + n / 2);  // 下半部分（小y）
    vector<Point> topGroup(pointsSortedByY.begin() + n / 2, pointsSortedByY.end());  // 上半部分（大y）

    sort(bottomGroup.begin(), bottomGroup.end(), compareByX);  // 下组按x升序排序
    sort(topGroup.begin(), topGroup.end(), compareByXDesc);  // 上组按x降序排序

    for (int i = 0; i < n / 2; ++i) {
        ll dx = abs(bottomGroup[i].x - topGroup[i].x);  // 计算x差的绝对值
        ll dy = abs(bottomGroup[i].y - topGroup[i].y);  // 计算y差的绝对值
        sumOption2 += dx + dy;  // 累加到总和
        int a = bottomGroup[i].id, b = topGroup[i].id;  // 获取索引
        if (a > b) swap(a, b);  // 确保a < b
        pairsOption2.push_back({a, b});  // 添加配对
    }

    // 选择总和更大的方案并输出配对
    if (sumOption1 >= sumOption2) {
        for (auto& pr : pairsOption1) {
            cout << pr.first << ' ' << pr.second << '\n';
        }
    } else {
        for (auto& pr : pairsOption2) {
            cout << pr.first << ' ' << pr.second << '\n';
        }
    }
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);

    int testCases;
    cin >> testCases;  // 读取测试用例数量
    while (testCases--) {
        processTestCase();  // 处理每个测试用例
    }
    return 0;
}

心路历程（WA，TLE，MLE……）

2025“钉耙编程”中国大学生算法设计暑期联赛（1）——1006景区建设（诈骗题）

Posted on 2025-07-19 Edited on 2026-02-17

思路讲解

$114∣x_1−x_2∣+5141∣y_1−y_2∣+919810|a_x - a_y|$

这个式子看着很唬人，其实还可以，不难发现

$114*100+5141*100-919810= -394310$

所以说就不用考虑这个x，也不用考虑这个y，只需要考虑高度（即a）就可以了。

找局部最高点也只需要这段代码（把旁边的四个格子扫一圈就行）。

FOR(i,1,N){
	FOR(j,1,M){
		// 已经布设了传送门
		if(i==1 && j==1) continue;
		bool isHi=true;
		FOR(k,0,3){
			ll x=i+dx[k],y=j+dy[k];
			if(x<1 || x>N) continue;
			if(y<1 || y>M) continue;
			if(A[x][y]>A[i][j]){
				isHi=false;
				break;
			}
		}
		if(isHi){
			ans+=cs;
			wp.EB(i,j);
		}
	}
}

AC代码

https://acm.hdu.edu.cn/contest/view-code?cid=1172&rid=19896

源代码

// Problem: 景区建设
// Contest: HDOJ
// URL: https://acm.hdu.edu.cn/contest/problem?cid=1172&pid=1006
// Memory Limit: 65536 MB
// Time Limit: 1000 ms
// by gospel_rock
// 
// Powered by CP Editor (https://cpeditor.org)

#include <bits/stdc++.h>
#define all(vec) vec.begin(),vec.end()
#define pb push_back
#define EB emplace_back
#define EM emplace
#define FI first
#define SE second
#define SZ(a) ((int) a.size())
#define FOR(i, a, b) for (int i = (a); i <= (b); ++i)
#define ROF(i, a, b) for (int i = (a); i >= (b); --i)
#define debug(var) cerr << #var <<":"<<var<<"\n";
#define lson(var) (var<<1)
#define rson(var) ((var<<1)+1)

using namespace std;

using  ll = long long;using ull = unsigned long long;
using DB=double;using LD=long double;
using i128 = __int128;
using pdd = pair<DB,DB>;using plb = pair<ll,bool>;
using pll = pair<ll,ll>;using vll=vector<ll>;using vb=vector<bool>;
using vpll=vector<pll>;
using arr3 = array<ll,3> ;using arr2 = array<ll,2>;
constexpr ll MAXN=(ll)1e2+10,INF=(ll)1e17+9; // 1e18+9也是素数
constexpr ll MAXV=(ll)1e6+10;
constexpr ll MAXM=(ll)2e6+10;
constexpr ll mod=(ll)1e9+7; // 998244353
constexpr double eps=1e-8;const double pi=acos(-1.0);
ll N,M,Q,X,K,T,lT;
ll A[MAXN][MAXN];
/*

*/
ll binpow(ll a,ll k){	// 迭代快速幂
	ll res=1;	// 防止a*a溢出
	while(k>0){
		if((k&1)==1) res=a*res;
		a=a*a;
		k>>=1;
	}
	return res;
}
const ll cs=binpow(2,34);
/* ->
  |  |
   <-
*/
ll dx[4]={0,1,0,-1};
ll dy[4]={1,0,-1,0};
std::istream& operator>>(std::istream& is, __int128& x) {
    std::string s;
    is >> s;
    x = 0;
    bool negative = false;

    // 处理负号
    size_t i = 0;
    if (s[0] == '-') {
        negative = true;
        i = 1;
    }

    // 将字符串转换为 __int128
    for (; i < s.length(); ++i) {
        x = x * 10 + (s[i] - '0');
    }
    if (negative) {
        x = -x;
    }
    return is;
}

std::ostream& operator<<(std::ostream& os, __int128 x) {
    if (x == 0) {
        return os << "0";
    }
    if (x < 0) {
        os << '-';
        x = -x;
    }

    // 将 __int128 转换为字符串
    std::string s;
    while (x > 0) {
        s.push_back((x % 10) + '0');
        x /= 10;
    }
    // 反转字符串以得到正确顺序
    for (int i = s.length() - 1; i >= 0; --i) {
        os << s[i];
    }
    return os;
}
bool cmp(const pll &a,const pll &b){
	return A[a.FI][a.SE]<A[b.FI][b.SE];
}
inline void __(){
	cin>>N>>M;
	FOR(i,1,N){
		FOR(j,1,M){
			cin>>A[i][j];
		}
	}
	i128 ans=0;
	vpll wp;
	FOR(i,1,N){
		FOR(j,1,M){
			// 已经布设了传送门
			if(i==1 && j==1) continue;
			bool isHi=true;
			FOR(k,0,3){
				ll x=i+dx[k],y=j+dy[k];
				if(x<1 || x>N) continue;
				if(y<1 || y>M) continue;
				if(A[x][y]>A[i][j]){
					isHi=false;
					break;
				}
			}
			if(isHi){
				ans+=cs;
				wp.EB(i,j);
			}
		}
	}
	wp.EB(1,1);
	sort(all(wp),cmp);
	FOR(i,0,SZ(wp)-2){
		auto &[x,y]=wp[i];auto &[xx,yy]=wp[i+1];
		ans+=114*abs(xx-x)+5141*abs(yy-y)+919810*abs(A[x][y]-A[xx][yy]);
	}
	cout<<ans<<"\n";
}

int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);cout.tie(0);
	
	cin>>lT;
	while(lT--)
		__();
	return 0;
}
/*
AC
https://acm.hdu.edu.cn/contest/view-code?cid=1172&rid=19896
*/

心路历程（WA，TLE，MLE……）

2025牛客暑期多校训练营2——Highway Upgrade（分治搜索，凸包）

Posted on 2025-07-19 Edited on 2026-02-17

思路讲解

START195B——Game (Hard)（分治搜索）打codechef还是有用的，或者说，非常有用。

先对询问进行排序，去重。

1 2	vll uq(all(que)); sort(all(uq));uq.resize(unique(all(uq))-uq.begin());

然后后面就是官解的思路，直接枚举边，然后可以得到一个sum→w（即sumW map）

map<ll,ll> sumW;
FOR(i,1,N){
	for(auto &[v,t,w]:g[i]){
		ll sum=revd[v]+dist[i]+t;
		if(sumW.count(sum)){
			sumW[sum]=max(sumW[sum],w);
		}else{
			sumW[sum]=w;
		}
	}
}

接着对这个sumW map进行一些处理，即使得sum越大，w也越大（sum大，w小的直接抛弃），接着装入sums数组（vector）中，方便我们二分查找。

for(auto it=sumW.begin();it!=sumW.end();){
	ll sum=it->first,w=it->second;
	if(w < mxw){  // 你的sum比别人大，你的w还比别人小，你存在的意义是什么？
		it++;
		sumW.erase(prev(it));  // 因此删除
	}else{
		mxw=w;   // 更新mxw
		it++;
	}
}
vector<pll> sums;
for(auto &[sum,w]:sumW){
	sums.emplace_back(sum,w);
}

接着利用该单调性，递归二分搜索。（这可能也是官解的思路，二分可行域）

vll cun(Q+5,0);
function<void(ll,ll,ll,ll)> dfs=[&](ll ia,ll ib,ll kc,ll kd)->void{
// kc，kd指的是询问的时候升级次数的编号，或者简单一点，就是uq的下标
	if(kc>kd) return;  
	ll k=uq[kc+kd>>1];
	ll km=kc+kd>>1;
	ll ans=INF;
	ll mid=0;
	// 在可行域中暴力计算答案
	FOR(i,ia,ib){  // ia ib 指的是sums的下标
		ll sum=sums[i].first,w=sums[i].second;
		ll lans=sum-w*k;
		if(lans<=ans){
			mid=i;
			ans=lans;
		}
	}
	
	cun[km]=ans;
	dfs(mid,ib,km+1,kd);
	dfs(ia,mid,kc,km-1);
};

AC代码

https://ac.nowcoder.com/acm/contest/view-submission?submissionId=78313063

源代码

// Problem: Highway Upgrade
// Contest: NowCoder
// URL: https://ac.nowcoder.com/acm/contest/108299/H
// Memory Limit: 1024 MB
// Time Limit: 6000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

#include <bits/stdc++.h>
#define all(vec) vec.begin(),vec.end()
#define pb push_back
#define SZ(a) ((int) a.size())
#define FOR(i, a, b) for (int i = (a); i <= (b); ++i)
#define ROF(i, a, b) for (int i = (a); i >= (b); --i)
#define debug(var) cerr << #var <<":"<<var<<"\n";
#define lson(var) (var<<1)
#define rson(var) ((var<<1)+1)
#define V vector

using namespace std;

using  ll = long long;using ull = unsigned long long;
using DB=double;using LD=long double;
// using i128 = __int128;
using pdd = pair<DB,DB>;using plb = pair<ll,bool>;
using pll = pair<ll,ll>;using vll=vector<ll>;
using arr3 = array<ll,3> ;using arr2 = array<ll,2>;
constexpr ll MAXN=static_cast<ll>(1e6)+10,INF=static_cast<ll>(1e17)+9; // 1e18+9也是素数
constexpr ll MAXM=(ll)1e6+10;constexpr ll MAXV=(ll)1e6+10;
constexpr ll mod=(ll)1e9+7; // 998244353
constexpr double eps=1e-8;const double pi=acos(-1.0);

ll N,M,Q,X,K,T,lT,A[MAXN];
/*

*/
inline void __(){
	cin>>N>>M;
	vector<vector<arr3> > g(N+5);
	FOR(i,1,M){
		ll u,v,t,w;
		cin>>u>>v>>t>>w;
		g[u].pb({v,t,w});
	}
	cin>>Q;
	vll que;
	FOR(i,1,Q){
		ll k;cin>>k;
		que.pb(k);
	}
	vll uq(all(que));
	sort(all(uq));uq.resize(unique(all(uq))-uq.begin());
	vector<ll> dist(N+5,(ll)1E18+17);
	auto dij=[&](ll s,vector<ll> &Dist,vector<vector<arr3> > & G)->void{
		using PLL=pair<ll,ll>;	// 改类型这里
		priority_queue<PLL,vector<PLL>,greater<> > q;
		vector<bool> vis(N+5,0);
		q.emplace(0,s);
		Dist[s]=0;
		while(!q.empty()){
			ll x=q.top().second;
			q.pop();
			if(vis[x]) continue;
			vis[x]=true;
			for(auto [y,w,t]:G[x]){
				if(Dist[y] > Dist[x] + w){
					Dist[y]=Dist[x]+w;
					q.emplace(Dist[y],y);
				}
			}
		}
	};
	dij(1,dist,g);
	vector<vector<arr3> > revg(N+5);
	FOR(i,1,N){
		for(auto &[v,t,w]:g[i]){
			revg[v].pb({i,t,w});
		}
	}
	vector<ll> revd(N+5,(ll)1E18+17);
	dij(N,revd,revg);
	map<ll,ll> sumW;
	FOR(i,1,N){
		for(auto &[v,t,w]:g[i]){
			ll sum=revd[v]+dist[i]+t;
			if(sumW.count(sum)){
				sumW[sum]=max(sumW[sum],w);
			}else{
				sumW[sum]=w;
			}
		}
	}
	ll mxw=0;
	for(auto it=sumW.begin();it!=sumW.end();){
		ll sum=it->first,w=it->second;
		if(w < mxw){
			it++;
			sumW.erase(prev(it));
		}else{
			mxw=w;
			it++;
		}
	}
	vector<pll> sums;
	for(auto &[sum,w]:sumW){
		sums.emplace_back(sum,w);
	}
	vll cun(Q+5,0);
	function<void(ll,ll,ll,ll)> dfs=[&](ll ia,ll ib,ll kc,ll kd)->void{
		if(kc>kd) return;
		ll k=uq[kc+kd>>1];
		ll km=kc+kd>>1;
#ifdef LOCAL
    debug(k);
#endif
		ll ans=INF;
		ll mid=0;
		FOR(i,ia,ib){
			ll sum=sums[i].first,w=sums[i].second;
			ll lans=sum-w*k;
// #ifdef LOCAL
   	// debug(sum);debug(lans);
// #endif
			if(lans<=ans){
				mid=i;
				ans=lans;
#ifdef LOCAL
    debug(km);debug(k);debug(lans);
#endif
			}
		}
		
		cun[km]=ans;
		dfs(mid,ib,km+1,kd);
		dfs(ia,mid,kc,km-1);
	};
#ifdef LOCAL
    cerr << "OK" << '\n';
#endif
	dfs(0,SZ(sums)-1,0,SZ(uq)-1);
#ifdef LOCAL
    for(auto &c:cun){
    	debug(c);
    }
    for(auto &u:uq){
    	debug(u);
    }
#endif
	FOR(i,0,Q-1){
		ll k=que[i];
		ll idx=lower_bound(all(uq),k)-uq.begin();
		cout<<cun[idx]<<"\n";
#ifdef LOCAL
    debug(idx);
#endif
	}
}

int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);cout.tie(0);
	
	cin>>lT;
	while(lT--)
		__();
	return 0;
}
/*
AC
https://ac.nowcoder.com/acm/contest/view-submission?submissionId=78313063
*/

心路历程（WA，TLE，MLE……）

vll cun(Q+5,0);

别写成N+5了。

2025“钉耙编程”中国大学生算法设计暑期联赛（1）——1009子序列（最长的两端比中间长子序列）

Posted on 2025-07-18 Edited on 2026-02-17

思路讲解

核心代码就这

	ll l=pos[0],r=pos[1];
	if(l>r) swap(l,r);
	ll ans=r-l+1;
	ll cnt=1;
	FOR(i,2,SZ(pos)-1){   // 一个点，要么在这个区间内，要么在这个区间外
		ll p=pos[i];   // 在这个区间内，必定不更新这个这个区间
		l=min(l,p),r=max(r,p);  // 在这个区间外，必定更新这个区间
		ans=max(ans,r-l+1-cnt);  // 所以说一个数更新区间的话，前面的数一定都在这个区间内
		++cnt;
	}
	cout<<ans<<"\n";
}

我们再来看看树状数组怎么做。

AC代码

https://acm.hdu.edu.cn/contest/view-code?cid=1172&rid=16219

源代码

// Problem: 子序列
// Contest: HDOJ
// URL: https://acm.hdu.edu.cn/contest/problem?cid=1172&pid=1009
// Memory Limit: 262144 MB
// Time Limit: 1000 ms
// by gospel_rock
// 
// Powered by CP Editor (https://cpeditor.org)

#include <bits/stdc++.h>
#define all(vec) vec.begin(),vec.end()
#define pb push_back
#define SZ(a) ((int) a.size())
#define FOR(i, a, b) for (int i = (a); i <= (b); ++i)
#define ROF(i, a, b) for (int i = (a); i >= (b); --i)
#define debug(var) cerr << #var <<":"<<var<<"\n";
#define lson(var) (var<<1)
#define rson(var) ((var<<1)+1)
#define V vector

using namespace std;

using  ll = long long;using ull = unsigned long long;
using DB=double;using LD=long double;
// using i128 = __int128;
using pdd = pair<DB,DB>;using plb = pair<ll,bool>;
using pll = pair<ll,ll>;using vll=vector<ll>;using vb=vector<bool>;
using arr3 = array<ll,3> ;using arr2 = array<ll,2>;
constexpr ll MAXN=(ll)2e6+10,INF=(ll)1e17+9; // 1e18+9也是素数
constexpr ll MAXV=(ll)1e6+10;
constexpr ll MAXM=(ll)2e6+10;
constexpr ll mod=(ll)1e9+7; // 998244353
constexpr double eps=1e-8;const double pi=acos(-1.0);

ll N,M,Q,X,K,T,lT,A[MAXN];
/*

*/
inline bool cmp(const ll &a,const ll &b){
	return A[a]>A[b];
}
inline void __(){
	cin>>N;
	vll pos;
	FOR(i,1,N){
		cin>>A[i];
		pos.pb(i);
	}
	if(N==1){	// 特判
		cout<<1<<"\n";
		return;
	}
	sort(all(pos),cmp);
	ll l=pos[0],r=pos[1];
	if(l>r) swap(l,r);
	ll ans=r-l+1;
	ll cnt=1;
	FOR(i,2,SZ(pos)-1){
		ll p=pos[i];
		l=min(l,p),r=max(r,p);
		ans=max(ans,r-l+1-cnt);
		++cnt;
	}
	cout<<ans<<"\n";
}

int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);cout.tie(0);
	
	cin>>lT;
	while(lT--)
		__();
	return 0;
}
/*
AC
https://acm.hdu.edu.cn/contest/view-code?cid=1172&rid=16219
*/