题目链接:https://ac.nowcoder.com/acm/contest/8563/C
题解:
首先我们要知道z^k=z+k-2*(z&k);
则只进行一次操作的代价:(r-l+1)*((z^k)+k-z)=2*(r-l+1)(k-z&k).
对于k-z&k,我们要令其尽可能的小,则z要是k的较大的子集。那么我们枚举k的子集并保存下来,对其排序。
然后我们枚举i=0->n,直到等差数列前i个数之和大于y
对于每一个i,我们计算其前i个数之和s[i],为了满足整个数列之和小于y,对于剩下的每一个数,要满足小于等于(y-s[i])/(n-i),即[0,(y-s[i])/(n-i)]
那么我们要在这个区间寻找最大的k的子集,即可对之前存下的k子集的数组进行二分查找,然后计算此时的答案。
具体见代码:
#include <bits/stdc++.h>
using namespace std;
const int maxn=1e6+10;
typedef long long ll;
ll n,a,d,k,y;
ll s[maxn];
ll si[maxn];
int cnt;
signed main() {
cin>>n>>a>>d>>k>>y;
ll x=a;
ll sub=k;
ll ans=1e18;
//枚举k的子集
do{
si[++cnt]=sub;
sub=(sub-1)&k;
} while(sub!=k);
sort(si+1,si+cnt+1);
//i=0的情况
ll minx=(y-s[0])/(n-0);
int p=upper_bound(si+1,si+1+cnt,minx)-si-1;
ans=ans=min(ans,1LL*2*(n-0)*(k-si[p]));
//1-n的情况
for(int i=1;i<=n;++i) {
if(i==1) s[i]=x;
else x+=d,s[i]=s[i-1]+x;
if(s[i]>y) break;
ll minx=(y-s[i])/(n-i);
int p=upper_bound(si+1,si+1+cnt,minx)-si-1;
ans=ans=min(ans,1LL*2*(n-i)*(k-si[p]));
}
cout<<ans<<endl;
}