Codeforces 963C. Cutting Rectangle

感谢lxy教会我这题qaq

---

w和h具体是什么不重要，先将他们离散化，然后把c[i]记成c[w][h]的形式

如果有某个c[w][h]=0一定不合法，并且$c[w][1]:c[w][2]:....c[w][h]$这个比例对所有w相同，这样才存在合法方案

单独考虑一种方块c[w][h]，可以把他放到c[w][h]的约数列，令$d[w]=gcd(c[w][i])$，那么就可以把c[w][i]放成d[w]的约数列

因为比例相同，对于每一个i，所有的$c[w][i]/d[w]$是相同的，也就是在把每个w展开成d[w]列下可以找到一组合法方案，然后每个w都可以放成d[w]的约数列，同时也是d[w]的约数行，因为要求行要对齐，行数要同时满足所有d[w]的约数
于是$ans=d(gcd(d[w]))$

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
#define ld long double
using namespace std;

ll gcd(ll a,ll b){return !a?b:gcd(b%a,a);}

const int maxn = 210000;
const ld eps = 1e-7;

int n;
int wn,hn;
map<int,ll>mp[maxn];
ll wi[maxn],hi[maxn],ci[maxn];

struct node
{
    ll x; ll *i;
    friend inline bool operator <(const node a,const node b){return a.x<b.x;}
}a[maxn];
void Trans()
{
    for(int i=1;i<=n;i++) a[i].x=wi[i],a[i].i=&wi[i];
    sort(a+1,a+n+1); a[0].x=a[1].x-1;
    wn=0;
    for(int i=1;i<=n;i++)
    {
        if(a[i-1].x!=a[i].x) ++wn;
        (*a[i].i)=wn;
    }

    for(int i=1;i<=n;i++) a[i].x=hi[i],a[i].i=&hi[i];
    sort(a+1,a+n+1); a[0].x=a[1].x-1;
    hn=0;
    for(int i=1;i<=n;i++) 
    {
        if(a[i].x!=a[i-1].x) ++hn;
        (*a[i].i)=hn;
    }
}

int main()
{
    //freopen("tmp.in","r",stdin);
    //freopen("tmp.out","w",stdout);

    scanf("%d",&n); ll g=0;
    for(int i=1;i<=n;i++) 
    {
        scanf("%lld%lld%lld",&wi[i],&hi[i],&ci[i]);
        g=gcd(g,ci[i]);
    }
    Trans();

    if((ll)wn*hn!=n) return puts("0"),0;

    for(int i=1;i<=n;i++) ci[i]/=g;
    for(int i=1;i<=n;i++) mp[wi[i]][hi[i]]=ci[i];

    for(int i=1;i<wn;i++)
    {
        for(int j=2;j<=hn;j++)
        {
            ld t1=(ld)mp[i][j]/(ld)mp[i][j-1];
            ld t2=(ld)mp[i+1][j]/(ld)mp[i+1][j-1];
            if(abs(t1-t2)>eps) return puts("0"),0;
        }
    }
    ll ans=0;
    for(int i=1;(ll)i*i<=g;i++) if(g%i==0)
    {
        ans++;
        if((ll)i*i!=g) ans++;
    }
    printf("%lld\n",ans);

    return 0;
}