8.12总结前日和今日

---

今天是组队赛的第三天，第一次打div2的比赛，表示真的有点难，还是有很长的路要走啊

---

8.11场链接

---

D题：（最小生成树）

题意：

n个点，每个点有一个权值，每两条边的权值为边的两个点的权值的gcd，求这幅图的最大生成树

解法：

首先处理出每个点权值的所有公因数，以所有的公因数为点，建一条点编号与其公因数的边，之后从大到小枚举每条边，跑一遍最大生成树即可

#include<cstdio>
#include<algorithm>
#include<vector>
#include<cmath>
using namespace std;
typedef long long ll;
const int maxn = 1e5 + 5;
vector<int>g[maxn];
int p[maxn];
int n;

int find(int x) {
	return (p[x] == x) ? x : (p[x] = find(p[x]));
}

int main() {
	freopen("dream.in", "r", stdin);
	int T; scanf("%d", &T);
	int kase = 1;
	while (T--) {
		for (int i = 0; i <= maxn; i++) {
			g[i].clear();
			p[i] = i;
		}

		scanf("%d", &n);
		for (int i = 1; i <= n; i++) {
			int x; scanf("%d", &x);
			for (int j = 1; j <= sqrt(x); j++) {
				if (x%j == 0) {
					g[j].push_back(i);
					if (x / j != j)g[x / j].push_back(i);
				}
			}
		}

		ll ans = 0;
		for (int i = maxn - 5; i >= 1; i--) {
			if (g[i].empty())continue;
			for (int j = 0; j < g[i].size()-1; j++) {
				int x = find(g[i][j]);
				int y = find(g[i][j + 1]);
				if (x != y) {
					p[y] = x;
					ans += i;
				}
			}
		}
		printf("Case %d: %d\n", kase++, ans);
	}
	return 0;
}

---

8.12场地址

J题：（二次剩余定理）

题意：有一个位数超过1e6的数，问这个数是否是完全平方数

解法：

取20个大质数，然后检查是否有解即可。

根据二次剩余，这等价于。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
using namespace std;
typedef long long ll;
const int N=1000010;
int n,i,cnt;char a[N];
bool isprime(int n){
    for(int i=2;i<=n/i;i++)if(n%i==0)return 0;
    return 1;
}
ll powmod(ll a,ll b,ll P){
    ll t=1;
    for(;b;b>>=1,a=a*a%P)if(b&1)t=t*a%P;
    return t;
}
bool work(ll n,ll p){
    if(n==0||n==1)return 1;
    if(powmod(n,p>>1,p)!=1)return 0;
    return 1;
}
bool check(int P){
    ll ret=0;
    for(int i=1;i<=n;i++)ret=(ret*10+a[i]-'0')%P;
    return work(ret,P);
}
int main(){
    scanf("%s",a+1);
    n=strlen(a+1);
    for(i=1000000000;;i++)if(isprime(i)){
        if(!check(i))return puts("No"),0;
        cnt++;
        if(cnt>=20)return puts("Yes"),0;
    }
}

---

E:

题意：（打印）

有一个长为L的母串和若干个长为n的子串，每次向前移动会增加一个单位的时间，换行也会增加一个单位的时间，问需要多个单位的时间才可以将整篇文章打印完

解法：

预处理出表示模式串第  个位置之后第一个字符  的位置，对于每一行计算需要的最大时间即可。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=1000010;
typedef long long ll;
int L,n,i,j,cur,f[100010][130];
char a[N],b[N];
ll ans;
int main(){
    gets(a);
    L=strlen(a);
    for(i=L-1;~i;i--)f[L][a[i]]=i;
    for(i=L-1;~i;i--){
        for(j=0;j<130;j++)f[i][j]=f[i+1][j];
        f[i][a[i]]=i;
    }
    while(gets(b)){
        n=strlen(b);
        int tmp=0;
        for(i=0;i<n;i++){
            char x=b[i];
            if(x==' ')continue;
            int nxt=((f[(cur+i)%L][x]-(cur+i))%L+L)%L;
            if(nxt>tmp)tmp=nxt;
        }
        tmp++;
        ans+=tmp;
        cur=(cur+tmp)%L;
    }
    printf("%lld",ans);
}