2019 멀티 - 대학 교육 대회 3
http://acm.hdu.edu.cn/contests/contest_show.php?cid=850
책 (1004) 배포
현재의 이분법이다 착수, 이분법 적 대답을 고려 (X 축 \)를 \ .
세트 (\ F_i \) 순방향 나타내는 \ (나는 \) 번째 세그먼트의 최대 개수로 분할 될 수 있고, 각 세그먼트는 것을 \ (\ leqslant X \) .
전송 분명히하는 열거 \ (J를 \) 경우 \ (s_i는-s_j \ leqslant X \)이 전사된다 (\ S_I \) 와 프리픽스이다.
어떤 수식이 될 것입니다 \ (S_I \는 leqslant X 축 + S_j \) 다음 각, \ (f_i \) 던져 값 \ (X + s_j \) 이 위치에, 단순히 각각의 질문에, 어떤 유지 보수를 살펴 그래서, 복잡성 \ (O (TN \ 로그 2N ^) \) .
#include<bits/stdc++.h>
using namespace std;
#define int long long
void read(int &x) {
x=0;int f=1;char ch=getchar();
for(;!isdigit(ch);ch=getchar()) if(ch=='-') f=-f;
for(;isdigit(ch);ch=getchar()) x=x*10+ch-'0';x*=f;
}
void print(int x) {
if(x<0) putchar('-'),x=-x;
if(!x) return ;print(x/10),putchar(x%10+48);
}
void write(int x) {if(!x) putchar('0');else print(x);putchar('\n');}
#define lf double
#define pii pair<int,int >
#define vec vector<int >
#define pb push_back
#define mp make_pair
#define fr first
#define sc second
#define FOR(i,l,r) for(int i=l,i##_r=r;i<=i##_r;i++)
const int maxn = 4e5+10;
const int inf = 1e9;
const lf eps = 1e-8;
const int mod = 1e9+7;
int n,a[maxn],f[maxn],s[maxn],k;
#define ls p<<1
#define rs p<<1|1
struct Segment_Tree {
int t[maxn*8];
void modify(int p,int l,int r,int x,int v) {
t[p]=max(t[p],v);
if(l==r) return ;int mid=(l+r)>>1;
if(x<=mid) modify(ls,l,mid,x,v);
else modify(rs,mid+1,r,x,v);
}
int query(int p,int l,int r,int x,int y) {
if(x<=l&&r<=y) return t[p];
int mid=(l+r)>>1,ans=0;
if(x<=mid) ans=max(ans,query(ls,l,mid,x,y));
if(y>mid) ans=max(ans,query(rs,mid+1,r,x,y));
return ans;
}
void build(int p,int l,int r) {
t[p]=0;if(l==r) return ;int mid=(l+r)>>1;
build(ls,l,mid),build(rs,mid+1,r);
}
}T;
int r[maxn],tot;
int check(int x) {
tot=0;
FOR(i,1,n) f[i]=0,r[++tot]=s[i],r[++tot]=s[i]+x;
sort(r+1,r+tot+1);int m=unique(r+1,r+tot)-r-1;
T.build(1,1,m);
for(int i=1;i<=n;i++) {
int pps=lower_bound(r+1,r+m+1,s[i])-r;
int ppsw=lower_bound(r+1,r+m+1,s[i]+x)-r;
int a=T.query(1,1,m,pps,m);
if(!a) f[i]=s[i]<=x?1:0;
else f[i]=a+1;
if(f[i]>=k) return 1;
T.modify(1,1,m,ppsw,f[i]);
}return 0;
}
void solve() {
read(n),read(k);
for(int i=1;i<=n;i++) read(a[i]);
for(int i=1;i<=n;i++) s[i]=s[i-1]+a[i];
int l=-2e14,r=1e9,ans=r;
while(l<=r) {
int mid=(l+r)>>1;
if(check(mid)) r=mid-1,ans=mid;
else l=mid+1;
}write(ans);
}
signed main() {
int t;read(t);while(t--) solve();
return 0;
}
1006 팬 블로그
아주 지루한 주제. . . 당신이 소수의 밀도와 윌슨의 정리되지 않습니다 테스트 (이름은 내가 문제에 대한 해결책을 알지 못했다 본 것입니다). .
어쨌든, 그건 프라임에 대한 말을하는 것입니다 \ (ρ-의 \) , 만날 수밖에 없다 \ ((는 P! - 2) \ p 형 PMOD \ 당량 1 \) .
덜 직접적인 폭력처럼, 난 그냥 Kuai 내 이전했다 (\ RM 밀러 \ _robon \ \를 ) 보드.
#include<bits/stdc++.h>
using namespace std;
#define int long long
void read(int &x) {
x=0;int f=1;char ch=getchar();
for(;!isdigit(ch);ch=getchar()) if(ch=='-') f=-f;
for(;isdigit(ch);ch=getchar()) x=x*10+ch-'0';x*=f;
}
void print(int x) {
if(x<0) putchar('-'),x=-x;
if(!x) return ;print(x/10),putchar(x%10+48);
}
void write(int x) {if(!x) putchar('0');else print(x);putchar('\n');}
#define lf double
#define pii pair<int,int >
#define vec vector<int >
#define pb push_back
#define mp make_pair
#define fr first
#define sc second
#define FOR(i,l,r) for(int i=l,i##_r=r;i<=i##_r;i++)
const int maxn = 1e6+10;
const int inf = 1e9;
const lf eps = 1e-8;
const int mod = 1e9+7;
const int pri[] = {0,2,3,5,7,11,13,23,31};
int gcd(int a,int b) {return !b?a:gcd(b,a%b);}
int mul(int x,int y,int p) {
int t=(long double)x*y/p;
int res=x*y-t*p;res%=p;if(res<0) res+=p;
return res;
}
int qpow(int a,int x,int p) {
int res=1;
for(;x;x>>=1,a=mul(a,a,p)) if(x&1) res=mul(res,a,p);
return res;
}
bool MR(int n) {
if(n<2) return 1;
int e=(n-1)>>__builtin_ctz(n-1);
for(int i=1;i<=8;i++) {
if(pri[i]>=n) break;
for(int d=e,lst=qpow(pri[i],d,n),now;d!=n-1;d<<=1,lst=now) {
now=mul(lst,lst,n);
if(lst!=1&&lst!=n-1&&now==1) return 0;
}if(qpow(pri[i],n-1,n)!=1) return 0;
}return 1;
}
int p,q;
void solve() {
read(p);int ans=p-1;
for(int i=p-1;;i--) {
if(MR(i)) {break;}
ans=mul(ans,qpow(i,p-2,p),p);
}
write(ans);
}
signed main() {
int t;read(t);while(t--) solve();
return 0;
}
1011 Squrirrel
물론 비교 바람직 세트는 뿌리 직경에서 찾을 수있다.
적은 직접 \ (DP \) , 먹거리 아웃의 직경 \ (F_ {I, 0 /. 1} \) 대표 \ (나는 \) 하위 트리가 최대 매직 트리 사용하지 않은 \을 (DP \) 종료 체인 거의에서 \ (DP \) 다시 같은, 아주 열심히 또한 코드를 작성합니다.
복잡성 \ (O (N-) \) .
#include<bits/stdc++.h>
using namespace std;
void read(int &x) {
x=0;int f=1;char ch=getchar();
for(;!isdigit(ch);ch=getchar()) if(ch=='-') f=-f;
for(;isdigit(ch);ch=getchar()) x=x*10+ch-'0';x*=f;
}
void print(int x) {
if(x<0) putchar('-'),x=-x;
if(!x) return ;print(x/10),putchar(x%10+48);
}
void write(int x) {if(!x) putchar('0');else print(x);putchar('\n');}
#define lf double
#define pii pair<int,int >
#define vec vector<int >
#define pb push_back
#define mp make_pair
#define fr first
#define sc second
#define FOR(i,l,r) for(int i=l,i##_r=r;i<=i##_r;i++)
const int maxn = 2e5+10;
const int inf = 1e9;
const lf eps = 1e-8;
const int mod = 1e9+7;
int dep[maxn],head[maxn],tot,n,nxt[maxn],id[maxn],vis[maxn],f[maxn][2],g[maxn][2],h[maxn][2],val[maxn],rt;
struct edge{int to,nxt,w;}e[maxn<<1];
void add(int u,int v,int w) {e[++tot]=(edge){v,head[u],w},head[u]=tot;}
void ins(int u,int v,int w) {add(u,v,w),add(v,u,w);}
void dfs(int x,int fa) {
if(dep[x]>dep[rt]) rt=x;nxt[x]=fa;
for(int i=head[x];i;i=e[i].nxt)
if(e[i].to!=fa) dep[e[i].to]=dep[x]+e[i].w,dfs(e[i].to,x);
}
void dp(int x,int fa) {
int mx1=0,mx2=0,ww,p;
for(int i=head[x];i;i=e[i].nxt)
if(e[i].to!=fa&&!vis[e[i].to]) {
dp(e[i].to,x);int t=f[e[i].to][0]+e[i].w;
f[x][0]=max(f[e[i].to][0]+e[i].w,f[x][0]);
if(t>=mx1) mx2=mx1,mx1=t,ww=e[i].w,p=e[i].to;
else if(t>mx2) mx2=t;
}
if(mx1==mx2) f[x][1]=f[x][0];
else f[x][1]=min(max(mx2,mx1-ww),max(mx2,f[p][1]+ww));
}
void solve() {
read(n);
for(int i=1,x,y,z;i<n;i++) read(x),read(y),read(z),ins(x,y,z);
dep[1]=0;dfs(1,0);int sz=0;dep[rt]=0;dfs(rt,0);
for(int u=rt;u;u=nxt[u]) id[++sz]=u,vis[u]=1;
for(int i=1;i<=sz;i++) dp(id[i],0);
g[1][0]=f[id[1]][0],g[1][1]=f[id[1]][1];
for(int i=2;i<=sz;i++) {
int ww;
for(int j=head[id[i-1]];j;j=e[j].nxt) if(e[j].to==id[i]) ww=e[j].w;//,printf("ww :: %d %d %d\n",id[i-1],id[i],ww);
g[i][0]=max(g[i-1][0]+ww,f[id[i]][0]);
g[i][1]=min(max(g[i-1][0],f[id[i]][0]),max(g[i-1][1]+ww,f[id[i]][0]));
g[i][1]=min(g[i][1],max(g[i-1][0]+ww,f[id[i]][1]));
}
h[sz][0]=f[id[sz]][0],h[sz][1]=f[id[sz]][1];int ans=1e9,p;
for(int i=sz-1;i;i--) {
int ww;
for(int j=head[id[i+1]];j;j=e[j].nxt) if(e[j].to==id[i]) ww=e[j].w;
h[i][0]=max(h[i+1][0]+ww,f[id[i]][0]);
h[i][1]=min(max(h[i+1][0],f[id[i]][0]),max(h[i+1][1]+ww,f[id[i]][0]));
h[i][1]=min(h[i][1],max(h[i+1][0]+ww,f[id[i]][1]));
}
for(int i=1;i<=sz;i++) {
//printf("%d ",id[i]);
// printf("%d %d %d\n",id[i],h[i][0],h[i][1]);
int tmp=min(max(g[i][0],h[i][1]),max(g[i][1],h[i][0]));
if(tmp<ans) ans=tmp,p=id[i];
}//puts("");
printf("%d %d\n",p,ans);
// write(ans);
}
void clear() {
for(int i=0;i<=n;i++) head[i]=g[i][0]=g[i][1]=h[i][0]=h[i][1]=f[i][0]=f[i][1]=vis[i]=0;
rt=tot=0;
}
int main() {
int t;read(t);while(t--) solve(),clear();
return 0;
}