ccpc网络赛2018

ccpc网络赛2018

体验极差，hdu老师去爬莫干山了然后OI赛制的好久。
傻逼题想多了没有出，6题滚了。

GuGu Convolution（hdu 6442）

题目描述
http://acm.hdu.edu.cn/showproblem.php?pid=6442

题解
ans= ((A+sqrt(B))^n-(A-sqrt(B))^n)/2
m模数不是质数，可能求不出逆元，不会除二，做不出来。
其实这个题这样搞出来之后就只有带根号的项了，所以我们只要算带根号的项就好了。
就是个傻逼题了。

代码

#include
#define ll long long
using namespace std;
int T,A,B,mod;ll n;
struct info{
  ll x,y,z;
  info operator*(const info &p)const{
    ll a=(x*p.x%mod+y*p.y%mod*z%mod)%mod;
    ll b=(x*p.y%mod+y*p.x%mod)%mod;
    return (info){a,b,z};
  }
}ans;
info Pow(info a,ll b)
{
  info res=a;b--;
  while(b)
  {
    if(b&1)res=res*a;
    a=a*a;b>>=1; 
  }
  return res;
}

ll cal(ll x)
{
  int res=1;
  for(int i=2;i*i<=x;i++)
  {
    int sum=0;
    while(!(x%i))x/=i,sum++,res*=i;
    if(sum&1)res/=i;
  }
  return res;
}

int main()
{
  scanf("%d",&T);
  while(T--)
  {
    scanf("%d%d%I64d%d",&A,&B,&n,&mod);
    ans=Pow((info){A,1,B},n);
    int num=cal(B);B/=num;
    ans.y=ans.y*(ll)sqrt(num)%mod;
    printf("1 %d %d\n",ans.y,B);
  }
  return 0;
} 

Search for Answer（hdu 6445）

题目描述
http://acm.hdu.edu.cn/showproblem.php?pid=6445

题解
给你一张竞赛图，对于所有的4元组，如果是四元环那么贡献是1，如果有两个点的出度是2那么贡献是-1。
我们先假设所有的四元组都是四元环，然后再把多余的扣掉。
这部分对答案的贡献是 A(n,4)。
这样，对于一个四元组，如果不是四元环的话，那么必然存在一个点，它的出度是2。
所以我们枚举每个度数>=2的点，它对答案的贡献是 A(deg[i],2)*(n-3)*4。就是对个这个点选两个点向其连边，再在其它点中随便选一个。是圆排列，所以乘了4。
所以 ans=A(n,4)-ΣA(deg[i],2)*(n-3)*4。
这东西其实就是要求 ΣC(deg[i],2)。
C(deg,2)=deg*(deg-1)/2，当前度数为n，度数每增加1，对于这个点贡献增加n。
暴力连边跑费用流就好了。

代码

#include
#define inf 1050000000
#define P 210
#define N 410
#define M 2010
using namespace std;
int Tx,n,S,T,tot,ans,d[P],du[P];char s[P][P];
int k,la[N],ff[M],w[N],pre[N],flag[N],q[N];
struct node{int a,b,c,v;}e[M];
void add(int a,int b,int c,int v)
{
  e[++k]=(node){a,b,c,v};ff[k]=la[a];la[a]=k;
  e[++k]=(node){b,a,0,-v};ff[k]=la[b];la[b]=k;
}

int next(int x){return ++x>405?1:x;}
bool spfa()
{
  for(int i=S;i<=T;i++)w[i]=inf;
  int l=1,r=2;q[1]=S;w[S]=0;
  while(l!=r)
  {
    int x=q[l];l=next(l);flag[x]=0;
    for(int a=la[x];a;a=ff[a])
      if(e[a].c&&w[e[a].b]>w[x]+e[a].v)
      {
        w[e[a].b]=w[x]+e[a].v;pre[e[a].b]=a;
        if(!flag[e[a].b])
          q[r]=e[a].b,flag[q[r]]=1,r=next(r); 
      }
  }
  return w[T]int flow=inf;
  for(int i=pre[T];i;i=pre[e[i].a])flow=min(flow,e[i].c);
  for(int i=pre[T];i;i=pre[e[i].a])
    ans+=flow*e[i].v,e[i].c-=flow,e[i^1].c+=flow;
}


int main()
{
  scanf("%d",&Tx);
  while(Tx--)
  {
    k=1;memset(la,0,sizeof(la));
    memset(ff,0,sizeof(ff));
    memset(d,0,sizeof(d));
    memset(du,0,sizeof(du));
    scanf("%d",&n);tot=n;ans=0;
    for(int i=1;i<=n;i++)
    {
      scanf(" %s",s[i]+1);
      for(int j=1;j<=n;j++)
      {
        if(s[i][j]=='1')d[i]++;
        if(s[i][j]=='2'&&i>j)
        {
          du[i]++;du[j]++;tot++;
          add(tot,i,1,0);add(tot,j,1,0);
        }
      }
    }
    S=0;T=++tot;
    for(int i=n+1;i1,0);
    for(int i=1;i<=n;i++)
    {
      ans+=d[i]*(d[i]-1)/2;
      for(int j=1;j<=du[i];j++)add(i,T,1,d[i]++);
    }
    while(spfa())add();
    printf("%d\n",n*(n-1)*(n-2)*(n-3)-ans*(n-3)*8);
  }
  return 0;
}