ニウケチャレンジ45——Cフレンド問題解決

タイトルリンク：https：//ac.nowcoder.com/acm/contest/8563/C

回答：

まず、z ^ k = z + k-2 *（z＆k）;を知る必要があります。

1回の操作のコスト：（r-l + 1）*（（z ^ k）+ kz）= 2 *（r-l + 1）（kz＆k）
。kz＆kの場合、できるだけ小さくする必要があります。 、次に、zがkのより大きなサブセットである場合。次に、kのサブセットを列挙して保存し、並べ替えます。

次に、等差数列の最初のi個の数の合計がyより大きくなるまで、i = 0-> nを列挙します。

各iについて、前のi数の合計s [i]を計算します。シーケンス全体の合計がy未満になるようにするには、残りの各数について、（ys [ i]）/（ni）、すなわち[0、（ys [i]）/（ni）]

次に、この間隔で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;
}