CF 830C(Bamboo Partition-满足\sum_{i=1}^n{ d\lceil ai/d\rceil-a_i } \leq k的d的最大值)
来源:互联网 发布:idea开发java图形界面 编辑:程序博客网 时间:2024/05/21 21:41
题意:求最大的正整数d,,使
观察发现只要
因此对每段分别求解。
分段点共
#include<bits/stdc++.h>using namespace std;#define For(i,n) for(int i=1;i<=n;i++)#define Fork(i,k,n) for(int i=k;i<=n;i++)#define Rep(i,n) for(int i=0;i<n;i++)#define ForD(i,n) for(int i=n;i;i--)#define ForkD(i,k,n) for(int i=n;i>=k;i--)#define RepD(i,n) for(int i=n;i>=0;i--)#define Forp(x) for(int p=Pre[x];p;p=Next[p])#define Forpiter(x) for(int &p=iter[x];p;p=Next[p]) #define Lson (o<<1)#define Rson ((o<<1)+1)#define MEM(a) memset(a,0,sizeof(a));#define MEMI(a) memset(a,127,sizeof(a));#define MEMi(a) memset(a,128,sizeof(a));#define INF (2139062143)#define F (1000000007)#define pb push_back#define mp make_pair #define fi first#define se second#define vi vector<int> #define pi pair<int,int>#define SI(a) ((a).size())#define Pr(kcase,ans) printf("Case #%d: %I64d\n",kcase,ans);#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;#define PRi2D(a,n,m) For(i,n) { \ For(j,m-1) cout<<a[i][j]<<' ';\ cout<<a[i][m]<<endl; \ } #pragma comment(linker, "/STACK:102400000,102400000")#define ALL(x) (x).begin(),(x).end()typedef long long ll;typedef long double ld;typedef unsigned long long ull;ll mul(ll a,ll b){return (a*b)%F;}ll add(ll a,ll b){return (a+b)%F;}ll sub(ll a,ll b){return ((a-b)%F+F)%F;}void upd(ll &a,ll b){a=(a%F+b%F)%F;}int read(){ int x=0,f=1; char ch=getchar(); while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();} while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();} return x*f;} #define MAXN (110)ll n,a[MAXN],k;ll calc(ll d) { ll p=0; For(i,n) { p+=(a[i]%d+d)%d; } return p;}bool check(ll d) { if (d<0) return 0; ll p=0; For(i,n) { p+=(a[i]%d+d)%d; }// cout<<d<<' '<<p<<' '<<endl; return p<=k;}ll work(ll l,ll r) { ll ans=-1; if (l>r) return ans; if (l==r) { if (check(l)) ans=l; return ans; } ll p=calc(l),delta=calc(l+1)-p; ll mm; if(!delta) { mm=r; } else { ll c=(k-p)/(delta); mm=min(l+c,r); } if (check(mm)) ans=mm; return ans;}vector<ll> S;int main(){// freopen("C.in","r",stdin);// freopen(".out","w",stdout); cin>>n>>k; For(i,n) a[i]=read(); sort(a+1,a+1+n); For(i,n) a[i]=-a[i]; ll ans=0,nxt; For(i,n) { ll nxt=1; for(ll pre=1;pre<=-a[i];pre=nxt+1) { nxt=a[i]/(a[i]/pre); S.pb(nxt); } } S.pb(1); sort(ALL(S)); S.erase(unique(ALL(S)),S.end()); vector<ll>::reverse_iterator it; work(35,36); ll pre=*S.rbegin(); for(it=S.rbegin();it!=S.rend();it++) { if (check(pre)) { ans=pre; break; } ll nxt=*it; ll p=work(nxt+1,pre); if (p!=-1) { ans=max(ans,p); break; } pre=nxt; } if (check(1)) ans=max(ans,1LL); ll tot=-a[n]*n; For(i,n) tot-=-a[i]; if(tot<=k) { ll p=(k-tot)/n; if(check(-a[n]+p)) ans=max(ans,-a[n]+p); } cout<<ans<<endl; return 0;}
阅读全文
0 0
- CF 830C(Bamboo Partition-满足\sum_{i=1}^n{ d\lceil ai/d\rceil-a_i } \leq k的d的最大值)
- printf("%d,%d,%d,%d\n",i++,--i,++i,i++);编译器是如何处理的。
- 关于C语言的printf("%d,%d,%d,%d\n",a,++a,a--,a--);
- 2.8.2 Wi n d o w s定义的U n i c o d e数据类型
- Σd|nφ(d)=n的证明
- n=∑ϕ(d) (d|n) 的证明
- Worm.Win32.D.i.s.k.G.e.n(磁.碟.机)的分析
- C,D,I 的C++实现
- U n i c o d e
- c++中printf("%5d%5d%5d%5d%5d", ++i, --i, i++, --i, -i--)的问题
- BZOJ 2301 Problem B(x属于[a,b],y属于[c,d]满足gcd(x,y)=k的(x,y)的有序对数)
- K D E是基于Troll Te c h n o l o g i e s公司的Q t库的
- CF 846 D. Monitor 滑动窗口,k*k子矩阵最大值
- printf("%d%d%d%d\n", a,b,c);
- 在O(n^2)内寻找矩阵A中 A(c,d) - A(a,b) 的最大值, c > a 且 d > b.
- 满足a^3 + b^3 = c^3 + d^3的abcd组合
- 2.9 成为符合A N S I和U n i c o d e的应用程序
- Windows 2000是使用U n i c o d e从头进行开发的
- IOS webview 自适应高度和加入css样式
- Linux系统安全配置基线
- CardStackView卡堆栈视图
- LeetCode Weekly Contest 40解题思路
- 代写材料的注意事项
- CF 830C(Bamboo Partition-满足\sum_{i=1}^n{ d\lceil ai/d\rceil-a_i } \leq k的d的最大值)
- 08-for循环
- Bootstrap列表组、面板和响应式嵌入组件
- tomcat
- iOS 面试题累计
- Python开发爬虫环境搭建
- 关于导出数据到记事本
- 是什么优化让 .NET Core 性能飙升?
- 第三章 活动图