1. 程式人生 > >[構造] Codeforces 618F Wunder Fund Round 2016 F. Double Knapsack

[構造] Codeforces 618F Wunder Fund Round 2016 F. Double Knapsack

給定兩個大小為n的可重集A, B ,兩個數集中的元素均為[1, n]的整數。
現在要求在兩個數集中各找出一個非空子集(子集也為可重集) ,滿足兩個集合中元素的和相等。
n<=10^6

加限制使得選的為兩個連續的段。
不妨設SA<=SB,記SA的字首和為SA(0),SA(1),SA(2),…, SB的字首和為SB(0),SB(1),SB(2),…
對於每個SA(i),找到最大的SB(j)滿足SA(i)>=SB(j),這樣可以得到n+1個SA(i)-SB(j),然而每個都是在0到n-1之間。
所以至少有兩個是一樣的。
於是SA(i)-SB(j)=SA(p)-SB(q),即SA(i)-SA(p)=SB(j)-SB(q)。

#include<cstdio>
#include<cstdlib>
#include<algorithm>
using namespace std;
typedef long long ll;
typedef pair<int,int> abcd;

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline void read
(int &x){ char c=nc(),b=1; for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1; for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b; } const int N=1000005; int n,a[N],b[N]; ll sa[N],sb[N]; int match[N]; int pos[N]; inline void print(int l,int r){ if (l>r) swap(l,r); printf
("%d\n",r-l); for (int i=l+1;i<=r;i++) printf("%d ",i); printf("\n"); } int main(){ int flag=0; freopen("t.in","r",stdin); freopen("t.out","w",stdout); read(n); for (int i=1;i<=n;i++) read(a[i]),sa[i]=sa[i-1]+a[i]; for (int i=1;i<=n;i++) read(b[i]),sb[i]=sb[i-1]+b[i]; if (sa[n]>sb[n]){ flag=1; for (int i=1;i<=n;i++) swap(a[i],b[i]),swap(sa[i],sb[i]); } sort(sb,sb+n+1); for (int i=0;i<n;i++) pos[i]=-1; for (int i=0;i<=n;i++){ int it=upper_bound(sb,sb+n+1,sa[i])-sb-1; match[i]=it; if (~pos[sa[i]-sb[it]]){ int t=pos[sa[i]-sb[it]]; if (!flag) print(t,i),print(match[t],match[i]); else print(match[t],match[i]),print(t,i); return 0; } else pos[sa[i]-sb[it]]=i; } return 0; }